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Abstract 

We present a review of the history and the present state of the fractal approach to the large-scale distribution of 
galaxies. The roots of the modern idea of cosmic fractality go to the beginning of the 20th century, when hierarchical 
world models were proposed by Fournier and Charlier. Fundamental aspects of hierarchical matter distribution was 
discussed by Einstein and Selety in the correspondence concerning inhomogeneous cosmological models. 

The Great Debate on the nature of spiral nebulae went over to a struggle on the structure and extension of 
the galaxy clustering. Already during the epoch of galaxy angular catalogues astronomers detected superclusters of 
galaxies on the sky. However, early counts of bright galaxies were close to the 0.6m-law, pointing to homogeneity. 
Also the debated variable extinction seemed to be a reasonable explanation of apparent galaxy clustering. Angular 
correlation function analysis gave small a homogeneity scale Rhom ~ 10 Mpc within which the correlation exponent 
corresponds to a fractal dimension D ~ 1.2 and outside which the galaxy distribution is homogeneous. 

It was realized later that a normalization condition for the reduced correlation function estimator results in 
distorted values for both -Rhom and D. Moreover, according to a theorem on projections of fractals, galaxy angular 
catalogues can not be used for detecting a structure with the fractal dimension D > 2. For this 3-d maps are 
required, and indeed modern extensive redshift-based 3-d maps have revealed the "hidden" fractal dimension of 
about 2, and have confirmed superclustering at scales even up to 500 Mpc (the Sloan Great Wall). On scales, where 
the fractal analysis is possible in completely embedded spheres, a power-law density field has been found. 

Two new fundamental cosmic numbers have appeared, the fractal dimension D and the crossover scale to 
homogeneity Rhom- Their values have been debated, and D = 1.2 indirectly deduced from angular catalogues 
has been replaced by D = 2.2 ± 0.2 directly obtained from 3-d maps and Rhom has expanded from 10 Mpc to 
scales approaching 100 Mpc. In concordance with the 3-d map results, modern all sky galaxy counts in the interval 
10 m -r 15 m give a 0.44m-law which corresponds to D = 2.2 within a radius of 100/i^qq Mpc. The narrow cones of the 
existing deep galaxy surveys and poorly known peculiar velocities at small scales are still the main limiting factors 
hampering precise estimates of D and Rhom- We emphasize that the fractal mass-radius law of galaxy clustering 
has become a key phenomenon in observational cosmology. It creates novel challenges for theoretical understanding 
of the origin and evolution of the galaxy distribution, including the role of dark matter and dark energy. 
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1 Introduction 

In this review we discuss the historical roots, methodological problems and modern results that branch of observational 
cosmology, which has given rise to the "Great Debate on Fractality" . This sharp and sometimes dramatic debate has 
been in the limelight almost the whole 20th century and still is. It has involved such persons as Einstein, Hubble, 
Sandage, Peebles, Charlier, Selety, Lundmark, de Vaucouleurs, Mandelbrot, Pietronero and many others who appear 
in the story. 

The debate on the fractality of the large scale galaxy distribution has been going on around two new fundamental 
empirical cosmic numbers, — the fractal dimension D and the bordering scale where fractality transforms into homo- 
geneity i?hom- The discussion of galaxy clustering started from scales 1 — 10 Mpc, then observations of the large scale 
structure (LSS) have shifted to the scales of 10 -j- 100 Mpc, and now we are entering gigantic scales of 100 -j- 1000 Mpc. 

There are several recent papers and books devoted to the LSS analysis, which are close to the subject of our 
review (Sylos Labini et al. 1998; Martinez & Saar 2002; Gabrielli et al. 2004; Jones et al. 2004). However, what has 
been lacking is a picture of historical steps in the studies of the large scale structure with a focus on methodological 
problems and challenges of current research. Our review addresses these questions. 

A fundamental task of practical cosmology is to study how matter is distributed in space and how it has evolved 
in cosmic time. The discovery of the strongly inhomogcneous spatial distribution of galaxies, at scales from galaxies 
to superclusters, i.e. over four orders of magnitude in scale, was of profound cosmological significance. Only faintly 
anticipated from photographic surveys, the surprisingly rich texture of galaxies became visible thanks to a large 
progress in measuring distances by rcdshifts for thousands of galaxies. The observed clustering is not just random 
clumping, but obeys a universal law. The two major theoretical tools, the correlation function and the conditional 
density analyses, both have revealed a power law dependence of the galaxy number density on the length scale, on 
scales 0.1-j- 10 Mpc. This signature of the scale invariance or fractality of the galaxy distribution has opened a new 
application of fractal geometry, widely used in modern statistical physics. 

In this review we present the growing evidence for fractality of the large scale galaxy clustering. In sect. 2 we give 
a brief introduction to the concept of fractal density field. In sect. 3 basic statistical methods for analysis of galaxy 
distribution are described. Sect. 4 is devoted to the main arguments for and against galaxy clustering, debated during 



the epoch of angular galaxy catalogues. Such 2-d data do not contain information on distance, and what is even 
more important, they are not sensitive to a distribution with the fractal dimension D > 2. In sect. 5 we summarize 
results obtained in the modern 3-d epoch from extensive redshift surveys which have revealed the "hidden" fractal 
dimension of about two, and have detected structures at scales up to 500 Mpc. One has seen how the value of D = 1.2 
derived from the angular catalogues has been replaced by D = 2.2 ± 0.2 in 3-d maps and the maximum observed 
scale of fractality Rhom has increased from 10 Mpc to scales approaching 100 Mpc. In sect. 6 we consider cosmological 
importance of the observed large scale fractality. 

2 The fractal view of the large scale structure of the Universe 

The last two decades have seen the first extensive surveys of galaxy redshifts, which have permitted us to move from 
the study of the angular distribution of galaxies on the celestial sphere to the analysis of their three-dimensional 
distribution in space. 

Already the first redshift surveys revealed a rich variety of structures in the galaxy universe. These have been 
characterized by astronomers using terms such as binaries, triplets, groups, rich, regular, and irregular clusters, walls, 
supcrclustcrs, voids, filaments, cells, soap bubbles, sponges, great attractors, clumps, concentrations, associations ... Of 
course each form of structure deserves separate studies, but they can be also viewed as natural appearances of one 
global master entity which is called a fractal. By using this self-similar structure one can describe an inhomogeneous 
galaxy distribution by means of one basic parameter, the fractal dimension D, which determines the global mass-radius 
behaviour of the Universe. 

However, it took a lot of time to recognize in the clumpy galaxy distribution on the celestial vault the signature 
of fractal structures. It is interesting and useful to see how the idea of self-similarity emerged in ancient times and 
developed into first protofractal models. 

2.1 The idea of a self-similar universe 
2.1.1 Protofractals 

The ancestors of fractals appeared in old cosmological thinking. In the presocratic Greece the philosopher Anaxagoras 
of Clazomenae (about 500 - 428 B.C.), who spent a large part of his life teaching at Athens, put forward his theory 
of seeds. Only fragments of his text have been preserved, but these give the impression that the theory was based on 
ideas reminiscent of self-similarity (e.g. Grujic 2001, 2002). Contrary to atomists, he regarded the matter divisible 
without a limit, and his famous words are: "In all things there is a portion of everything." 

What has been called by Mandelbrot a protofractal, appeared as simple hierarchies in the 18th century works by 
Swedenborg, Kant and Lambert. Emanuel Swedenborg (1688-1772) was a Swedish scientist and visionary who was 
appointed, at the age of 28, assessor extraordinar in the Swedish College of Mines. A very productive thinker and 
writer, he discussed practically all fields of science of his time. In 1734, in his Principia, Swedenborg put forward the 
remarkable view of self- similarity and cosmic hierarchy, which was related to his general opinion that everything in 
the world is constructed according to a common plan. 

The next steps were made by Immanuel Kant (1724-1804) and Johann Lambert (1728-1777) with their hierarchic 
systems. For example, Kant wrote in his book available in translation as Kant (1755): We see the first members of a 



progressive relationship of worlds and systems; and the first part of this infinite progression enables us already to recognize what 
must be conjectured of the whole. There is no end but an abyss ■ ■ ■ without bound. 

But there was an essential difference between these world models. Kant, and perhaps Swedenborg, imagined that 
the hierarchy continues without end towards larger levels of celestial systems: it was an infinite hierarchy. Lambert 
thought that after a large (he cites the figure 1000 as an example), but finite number of steps, the hierarchy ends. He 
thought that the stellar systems are kept together by the gravity of a dark mass: ... in the end you arrive at the middle 
point of the whole world structure and there I find my ultimate mass which governs the whole creation. 

John Herschel (1792-1871) was, among many others, intrigued by the dark sky or Olbers's Paradox (Harrison 
1987). And he outlined a solution in a private letter: ... it is easy to imagine a constitution of a universe literally infinite 
which would allow of any amount of such directions of penetration as not to encounter a star. Granting that it consists of 
systems subdivided according to the law that every higher order of bodies in it should be immensely more distant from the centre 
than those of the next inferior order - this would happen. 

Herschel had in mind some kind of a hierarchical system. In another text he gives as examples the satellites of the 
planets of the solar system and the large distances between stars and asserts that "the principle of subordinate grouping" 
assumes "the character and importance of a cosmical law". 

2.1.2 Fournier d'Albe and Carl Charlier 

The idea of a hierarchic structure of the universe was blown into life again by Edmund Fournier d'Albe (1868-1933). 
In 1907 he published the remarkable book Two new worlds where one finds a mathematical description of a possible 
hierarchical distribution of stars. In Fournier's world (Fig^) the stars are distributed in a hierarchy of spherical 
clusters in an infinite space, so that the mass inside each sphere increases directly proportionally to its radius: 



Of course, this makes it dramatically different from the mass-radius behaviour in a homogeneous universe where 
M(R) (x i? 3 . This was Fournier's idea of how to avoid cosmological paradoxes appearing in Newton's universe. 

The cosmic hierarchy was further studied by the Swedish astronomer Carl Charlier (1908) in his article How an 
infinite world may be built up. After enthusiastic inspection of Fournier's book, he developed more general models of 
stellar distributions which also solve the Olbers's paradox and the riddle of infinite gravitational potential. 

Charlier found a criterion which the hierarchy must fulfil so that it will solve the mentioned paradoxes. The decisive 
factor is how fast the density decreases from one level (i) to the next one (i+1), and this depends on the ratio of the 
sizes of the successive elements and on the number JVj+i of the lower elements forming the upper element. Denoting 
the sizes (radia) with Ri and Ri+i, Charlier's first criterion may be written as 



or the size of upper level element divided by size of lower level element is larger or equal to the number of lower elements 
forming the upper elements. For example, Fournier's illustration fulfils Charlier's first criterion: there iVj = Ni + i = 5, 
and the size ratio is always about 7. 

In Charlier (1922), after a note by Selety (1922), a second criterion was derived: 



M(R) cx R. 



(1) 



Ri+i/Ri > 



(2) 



Ri+i/Ri > v 7 ^ 



(3) 




Figure 1: Founder's (1907) regular hierarchic world model. The same elements are repeated on different levels in a self-similar way. The 
average density p(R) within the scale R decreases with increasing scale R, hence p(R) — » for R — > oo. 

In terms of a continuous mass-radius behaviour for identical particles with mass m the first criterion corresponds to 
M(R) = mN(R) oc R 1 . The second criterion, M{R) oc R 2 , and it is sufficient to cope with Olbers' paradox and the 
infinite gravity force. The original Fournier's (and the first Charlier's) criterion is a stronger condition, and allows 
also a finite gravitational potential and finite stellar velocities. 

2.2 Genuine fractal structures 

Even though Fournier's and Charlier's hierarchies were overly simple for the real world, as we now know, they contained 
the seeds of the modern concept of genuine fractal. Regular hierarchical models (like in Fig^) have a number of 
preferred scales, corresponding to the sizes of clusters at the level "«" . A great advantage of stochastic fractals is that 
they can be used to model scale- invariant galaxy clustering without preferred scales (Fig|5J). 

Though fractal geometry emerged just a few decades ago, some elements of it can be found already in the works 
of Poincare and Hausdorff about a century ago. Mandelbrot (1975) introduced the name "Fractals" and gave the 
following definition: A fractal is a set for which the Hausdorff dimension strictly exceeds the topological dimension. 

Mandelbrot realized that fractal geometry is a powerful tool to characterize intrinsically irregular system. Nature 
is full of strongly irregular structures; trees, clouds, mountains and lightnings are familiar objects, which have in 
common the property that if one magnifies a small portion of them, a complexity comparable to that of the entire 
structure is revealed. This is geometric self-similarity. 

Fractals are simple but subtle, as Luciano Pietronero likes to say, and becoming familiar with them requires some 
training of intuition. Here we briefly describe the essential properties of fractals, which are used in studies of the galaxy 
distribution. We recommend Mandelbrot (1982) for an original presentation by the father of fractals and Falconer 
(1990) for a special mathematical treatment. Here we emphasize observational consequences of self-similarity so that 
if this important property is actually present in galaxy data, one will be able to detect it correctly. 



Figure 2: An example of a genuine stochastic fractal structure with the fractal dimension D = 1.26 in two-dimensional space (B. 
Mandelbrot 1989). 



2.2.1 Self-similarity and power law. 

The difference between a self-similar distribution and a distribution with an intrinsic characteristic scale was clearly 
discussed by Coleman & Pietronero (1992). From a mathematical point of view self-similarity implies that the rescaling 
of the length r by a factor b 

r — » r = br (4) 

leaves the considered property, presented by an arbitrary function f(r), unchanged apart from a renormalization that 
depends on b but not on the variable r. This leads to the functional relation 

f(r') = f(b-r) = A(b)-f(r) (5) 

This is satisfied by a power law with any exponent. In fact for 

f(r) = / r Q (6) 

we have 

f(r') = fo(br) a = (b) a f(r). (7) 

Here the exponent a defines the behaviour of the function everywhere. There is no preferred scale. It is true that the 
condition on the amplitude /(ro) = 1 implies a certain length ro 

ro = fo 1/a , (8) 

and one might be tempted to call r a characteristic length, but this is quite misleading for self-similar structures! 
The power law refers to a fractal structure that was at the beginning constructed as self-similar and therefore cannot 
posses a characteristic length. 

In eqlHlthe value of ro is just related to the amplitude of the power law, and the amplitude has nothing to do with 
the scaling property. Instead of the value 1 for the amplitude, one could have used any other number in the relation 



/( r o) = 1 to obtain other lengths. This is a subtle point of self-similarity; there is no reference value (like the average 
density) with respect to which one can define what is big or small. 

This behaviour of the power-law is in contrast with the case of the exponential decay function where there is an 
intrinsic characteristic parameter ro 

g(r) = 9o e- r/ro (9) 

which determines a preferred length scale for the behaviour of the function. In the power-law the constant, dimen- 
sionless exponent a is not related to length scales at all. 

2.2.2 Fractal dimension 

The basic characteristics of a fractal structure is its dimension. The fractal dimension is a measure of the "strength 
of singularity" around the structure points. If there is a zero-level of structure elements, as in physical fractals where 
there is no mathematical singularity, the rate of growth of density with a decreasing spatial scale still defines the 
fractal dimension. 

Let us consider the simplest example of a regular fractal structure to illustrate how the fractal dimension can be 
determined. Starting from a point occupied by an object we count how many objects are present within a sphere of 
radius r in order to establish a number-radius relation from which one can define the fractal dimension D. Suppose 
that in the structure of Fig^we can find Nq objects with mass mo in a volume of size ro. If we consider a larger 
volume of size r*i = k r ■ ro we will find Ni — fejy • No objects. In a self-similar structure the parameters k r and /cat will 
be the same also for other changes of scale. So, in general in a structure of size r n = fe™ • ro we will have N n = k 7 ^ ■ No 
objects. We can then write the relation between the number TV (or mass M = moN) and length r in the form 

N(r) =B -r D , or M(r) = m B ■ r D (10) 

where the fractal dimension is the exponent D of the power law, i.e. 

D = ^ (11) 
log k r 

depends on the rescaling factors k r and k^. The prefactor B is related to the zero-level parameters Nq and ro, 



B=^. (12) 

We note that EqO corresponds to a smooth continuous presentation of a strongly fluctuating function as is evident 
from Figin The smooth average power-law for a fractal structure is always accompanied by large fluctuations and 
clustering at all scales. 

Simple algorithms for constructing stochastic fractal structures with given fractal dimension one may find e.g. in 
Gabrielli et al. (2004). With these algorithms one may build artificial galaxy catalogues useful for testing different 
methods of structure analysis. 

2.3 Concepts of density field 

The definition of density field is the basis of the large scale structure analysis. There is an essential conceptual 
difference between ordinary and fractal density fields. The former kind of model is usually utilized for the description 
of gas or fluid dynamics having short range of correlations, while the latter one emerges in physical systems with 
strong long range scale-invariant fluctuations. 



2.3.1 Ordinary fluid-like density fields 

The concept of the density of a continuous medium (approximating fluid or gas), as it is normally used in hydrody- 
namics, contains the assumption that there exists a value of density independently of the size of the volume clement 
dV. Then one may define the density p[x) at a point x and regard it as a usual continuous function of position in 
space: 

_ ,. M(x,V) 

P fluid = P( X ) = 1™ y ( 13 ) 

where M is the mass of the fluid inside the volume V around the point x. For ordinary continuous media the limit 
exists and does not depend on the volume V , because at sufficiently small scales homogeneity is reached, i.e. M = pV . 

When one studies fluctuations of an ordinary fluid, one can regard q(x) as one realization of a stochastic process 
for which the usual statistical moments are defined - average, dispersion etc. 

Models for ordinary density fields go under the common name fluid-like correlated distributions. These have a 
uniform average background with superimposed fluctuations that are correlated. An example of such a fluctuation 
may be seen in Fig|3 top. Main properties of such distributions, in comparison with fractal distributions, are discussed 
by Gabrielli et al. (2004). 

An ordinary stationary stochastic density field p{x) may be represented as a sum of density fluctuations Sp(x) and 
the mean constant density po = (p(x)), so that 

p{x)=p + Sp, (14) 
or in terms of the dimensionless relative density fluctuation: 

* (f) = p_p^ = 5 A*1 . (15) 

Po Po 

Note that the relative fluctuation 8{x) in Eg 1151 has positive and negative values in the range — 1 < S(x) < oo, 
while the density field is always positive p(x) > for positive masses of particles. 

One usually considers 8p{x) as a realization of a Gaussian stochastic process, for which the phases of fluctuations 
are uncorrelatcd. Here a fundamental role is played by the average density p — const > which should exist and 
be well defined and positive for each realization of the random process p(x). "Well defined" means that the average 
density does not depend on the size and location of a fair test volume. 

2.3.2 Ordinary stochastic discrete processes 

Ordinary density field may be also presented by a discrete stochastic process, called a stochastic point process or a 
point-particle distribution. Here discreteness introduces some new aspects related to point-like singularity of particles. 

An important example of a homogeneous stochastic discrete density field is the Poisson process, giving rise to a 
number density of particles n(x) = S^ 1 (5(x — Xi). According to the Poisson law the probability P to find N particles 
in a volume V is 

. (N) N exp(-(A0) , , 

P(N,V)= K ' £ { (16) 

where (N) = uqV is the average number of particles in the volume V(r). The only parameter of the Poisson distribution 
is the constant number density no = (n(x)), or intensity of the Poisson process. It determines a characteristic scale 
Ao for this process 

A W % 1/3 « Rsep , (17) 



which is about the average distance between the particles also denoted as R sep - 

The normalized number (mass) variance a 2 (r) in a sphere with radius r is defined as 



a 2 (r) 



(N(r) 2 ) - (N{r)) 2 



(18) 



(iV(r)) 2 

which is an important quantity to characterize a stochastic process at both small and large spatial scales r. 

Discreteness noise. Finite distances between point-particles of a discrete stochastic process make an essential 
noise of discreteness or shot noise appear when considered scales are less that average distance between particles 
r < RfL- This noise increases with decreasing r as 

-3/2 



1 



a(r) « — « ( - 



N 



(19) 



and becomes infinitely large at very small scales. 

Large spatial scales. On large spatial scales r the distribution becomes homogeneous with its normalized dispersion 
fea ll9|l approaching zero as 1/ \/~N when the number of particles increases indefinitely. The scale i?hom ~ Ao is the 
homogeneity scale of the Poisson process, it is defined from the condition that the normalized number variance 

C 2 (i?hom) = 1- 

The Poisson distribution is a homogeneous stochastic discrete density field without correlations, so that its correla- 
tion function is zero (neglecting the so-called " diagonal" term corresponding to the point singularity at r = 0) . This 
means that the positions of points are independent of each other, so that there are no genuine structures. Though the 
eye may see apparent structures, these are just random fluctuations in an outcome (realization) of the process. 

Gabrielli et al. (2004) considered another important example, the superhomogeneous discrete process like particles 
in a lattice with small correlated shifts of the particles around the regular lattice knots. Such a process is used for 
generating initial conditions for cosmological N-body simulations. 
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Figure 3: Top: an ordinary (fluid-like) discrete density field: a fluctuation on a Poisson background; bottom: a stochastic fractal density 
field (from Sylos Labini et al. 1998). 



2.3.3 Fractal density fields 

In the large scale structure analysis one applies a concept of discrete fractal point processes for which one may consider 
the density field in the form of a spatial distribution of N particles in positions x a with masses m a in a volume V, 
then 

N 

p(x) =^2m a S(x - x a ) . (20) 

a=l 

In the case of identical particles with m a — rn one may simply use the concept of the number density field 

N 

n{x)=^5{x-x a ). (21) 

a=l 

In order to describe the continuous hierarchy of clustering, which is a new characteristic of a fractal stochastic 
process, we first take into use a new independent variable, the radius r of the spherical volume V(r) in which the 
particles are counted. Then the fractal mass density may be defined as a function of two variables, the position x a of 
a particle and the radius r of the volume. 

M{Xg,V) 

Pftacw = Pi x a, r) = , (22) 

where M{x a , V) is the mass inside the volume V(r) located around the particle position x a . 

For a mathematical infinite ideal fractal the number of points of the structure in a finite volume is infinite, so the 
mass is also infinite. In physics this problem may be avoided by some natural lower limit in sizes of elements, making 
the zero-level of the hierarchy. In this case one may speak about basic structure elements — point mass particles. 
Then the number of particles and the mass within the volume V(r) is finite, i.e. in this case Ea l22l defincs a physically 
measurable density. Though the quantity p{x a , r) is highly fluctuating from one particle position to another, as we 
shall see below it is possible to consider statistical average which is more stable characteristics of a fractal structure. 

In Fig|3|an ordinary (fluid-like) density fluctuation on a Poisson background is compared with a stochastic fractal 
density fluctuations. In case of fractal structures the ordinary concept of the mass density of a continuous medium 
is not applicable. This is because the mass density can be defined only if both the position x and the volume V are 
considered. In every volume V containing a part of the structure there is a hierarchy of clusters and the value of the 
mass density strongly depends on the size of the volume. This is totally different from the usual calculus. Now there 
is no limit for the mass-volume ratio fEa ll 3|l and the density increases indefinitely when the volume tends towards 
zero: 

-°o, fOT V^O} . (23) 
If the basic zero level elements exist then it determines the maximum value of the fractal density for the structure. 

2.4 Exclusive properties of fractal density fields 
2.4.1 Power-law density-radius relation 

A characteristic feature of the fractal density field is a power-law behaviour of the density with increasing of the radius 
of the spherical volume V(r) centered at a structure point. We can illustrate this property for the case of regular 
fractal structure (Fig.l), where the number of subelements within an element of the higher level is given by Ea llUI So 



in continuous representation, using Ea l22l one can define the fractal number density, which related to the elements of 
radius r, as 

»-> - m - < 24 > 

where D is the fractal dimension of the structure. 

Hence, in order to calculate the fractal density one can start from a basic element which belongs to the structure 
and count the number of objects within the sphere V(r) and divide it by the volume of the sphere V(r) centered at the 
point. Apart from some fluctuation, depending on the actual position of the point within the structure the power-law 
En Ell will be obtained. 

Note that the result is unexpected for our usual intuition — the density decreases from each point of the fractal 
structure. It seems like each point were the centre of the structure from which the density decreases outwards following 
the law in Eq|^| This property is essential for cosmological implications of the fractal structures (sec. 6. 5). 

2.4.2 Massive, zero-density universes 

An important consequence from Eg 1241 is that in infinite space the fractal density field differs from an ordinary fluid-like 
density field at the limit of large volumes V where 

{P M -0, for V^oc} (25) 

This property is due to a growing dilution of the hierarchy with increasing scales, so that a fractal structure is 
asymptotically dominated by voids. 

Hence an infinite fractal universe can contain an infinitely large number of objects (hence an infinite mass) simul- 
taneously with the zero density of the whole Universe. This unusual property of a hierarchical structure was exploited 
in old cosmological models to avoid gravitational and photometric paradoxes of the Newtonian infinite universe. This 
follows from the relations <p oc M/R and F cx M/R 2 . 

2.4.3 The role of lower and upper cutoffs 

In the realm of physics real structures usually have a lower scale i? m i n and an upper scale i? max between which the 
physical system follows fractal self-similar behaviour. These scales are called lower and upper cutoffs. 

In studies of large-scale galaxy distribution the lower cutoff i? m i n is assumed to be equal to the size of a galaxy, 
while galaxies play a role of point-like particles. For different cosmological problems there could be different choices 
of the lower cutoff: dark matter clumps of 10 6 ^ s Mq, stars, comet-size objects, atoms, elementary particles. So the 
lower cutoff is usually a well-defined quantity. 

The upper cutoff presents a much more complicated problem in studies of the galaxy distribution. In principle, one 
should apply such methods of the large scale structure analysis which allow one to determine directly from a galaxy 
survey the scale R max = Rhom where the galaxy distribution becomes homogeneous. However, such methods need a 
large survey volume whose size should correspond to several times the scale. 

Up to now the largest galaxy redshift surveys cover a small part of the sky which hinders a firm estimation of the 
size i?hom- Is there an upper cutoff for the large-scale galaxy distribution and what is its value? These are the primary 
questions around which the most acute discussion is going on. 



2.4.4 Lacunarity 

Two fractal structures with the same fractal dimension may look very different. In particular, when one makes fractal 
models of the galaxy distribution, it is quite essential how large a relative volume is occupied by voids, on a given 
scale. 

This property was termed lacunarity by Mandelbrot (1982), from the word "lacuna" meaning hole or gap in Latin. 
Quantitatively lacunarity may be characterized by the constant of proportionality F in the relation 

N v (\ > A) =FA-° (26) 

where N v is the number of voids with size A > A within a fixed volume inside the structure. 

Another definition was introduced by Blumenfeld & Mandelbrot (1997). They used the variation of the prefactor 
B(x a ) in the density law computed for each structure point x a within a ball with a fixed radius R, N(r,x a , R) = 
B(x a , R) r D (Ea 110(1 . Then lacunarity is defined as a normalized dispersion of the distribution of the prefactor 

..fflL^k (27, 

Concrete examples of structures inside fixed sample volumes, having different lacunarities according to this definition 
may be found in Martinez & Saar (2002). 

We note that the high lacunarity of the Rayleigh- Levy flight fractal was the reason why it was rejected as a model 
for the real distribution of galaxies (Peebles 1980). However, later Mandelbrot (1998) demonstrated that fractal 
structures with a small lacunarity resemble more closely the arrangement of galaxies (see FigEJ - 

2.4.5 Projection and intersection 

The properties of orthogonal projections and intersections of a fractal structure play an important role in the analysis 
of galaxy samples with different geometries, both from angular 2-d and spatial 3-d catalogues. 

Orthogonal projection. Let an object (structure) with a fractal dimension D, embedded in an Euclidean space 
of dimension d = 3, be orthogonally projected onto an Euclidean plane with dl — 2. Then according to a general 
theorem of fractal projections (see Mandelbrot 1982; Falconer 1990), the projection as a fractal object receives the 
fractal dimension D pr so that 

D. pr = D if D < 2 (28) 

and 

D pr = 2 if D > 2 . (29) 

This means that in 3-d space a cloud having the fractal dimension D k 2.5 satisfies eqlJHl and hence gives rise to 
a homogeneous shadow (D pr = 2) on the ground. Consequently, the orthogonal projection hides from view fractal 
structures with D > 2. This has an important implication for the apparent distribution of galaxies on the sky (sec. 
4.4). 

Intersection of a fractal. If an object with a fractal dimension D, embedded in a d = 3 Euclidean space, intersects 
an object with the dimension D', then according to the law of co-dimension additivity (see Mandelbrot 1982; Falconer 
1990), the dimension of the intersection D- lnt becomes 

A„t =D + D'-d. (30) 



For example, if a fractal structure with D = 2 in 3-d space is intersected by a plane with D' = 2, then we obtain 
for the fractal dimension of the thin intersection D m t = 2 + 2 — 3 = 1. This property of intersections explains why a 
fractal structure with D w 2 may look as a fractal with D w 1 when inspected on large scales from a sample coming 
from a thin slice-like galaxy survey. 

2.4.6 Multifractal structures 

In fractal models of the galaxy distribution one usually utilizes only spatial positions for a sample of galaxies. This 
allows one to describe the distribution by means of only one parameter — the exponent in the power-law or the fractal 
dimension of the structure. 

Real observational data contain also other important astrophysical information for each galaxy, such as luminosity, 
morphology, spectral properties, stellar contents e.t.c. In this case the scaling properties can be different for different 
types of galaxies. To take into account the dependence of the distribution on these parameters one has to introduce 
a more general model, called multifractal structures, which is characterized by a continuous set of fractal dimensions. 
Such approach was firstly suggested for galaxy distributions by Pietronero (1987). For recent discussions of this subject 
see Gabrielli et al. (2004), Martinez & Saar (2002) and Jones et al. (2004). 

We note that multifractality may be viewed differently thanks to the complexity of the problem (even for fractals 
there is no unique definition). Multifractals are in contrast with homogeneity exactly like fractals are. In fact, the 
multifractal picture is a refinement and generalization of fractal properties (Paladin & Vulpiani 1987; Benzi et al. 
1984). 

2.5 Modern redshift and photometric distance surveys 

For many years, astronomers could make only indirect conclusions about the distribution of galaxies on the basis of 
their two-dimensional projected locations on the celestial sphere. Such studies of projections are well reviewed in 
Peebles (1980). 

In recent years the situation was completely changed when it became possible to measure the 3-dimensional 
distribution of galaxies using data from massive surveys of galaxy redshifts. At the present time there are several 
approaches for investigating the space distribution of matter (luminous and dark): photometric distance measurements, 
extensive redshift surveys of galaxies and quasars, the analysis of counts of galaxies, and the study of image distortion 
effects produced by weak gravitational lensing. All these observational studies have shown that clustering is a common 
phenomenon in the realm of nebulae. 

2.5.1 Redshift surveys 

Nature has given the astronomer, in the form of the linear redshift- distance law, an accurate way to measure 
extragalactic distances, which is generally more precise than photometric methods. For example, for a velocity 
dispersion of 50 km/s, typical for field galaxies, one can measure distances with an accuracy of about 1 Mpc. In order 
to get in this way a deep, 3-dimensional map of the surrounding galaxy universe, it is necessary to make deep surveys 
of redshifts, complete up to sufficiently faint magnitude limits. 

Over 2700 galaxies had their redshift listed in the Second Reference Catalogue by de Vaucouleurs, de Vaucouleurs 
& Corwin (1976). This was the breakthrough which made it possible to use redshifts for mapping the structures made 



by galaxies. 

Giovanelli & Haynes (1991) emphasized that "In the last fifteen years, advances in detector and spectrometer 
technology at both optical and radio wavelengths have spurred a tremendous explosion in the galaxy redshift tally." 
Indeed, this explosion has continued with an exponential rate up to present. Currently more than one million redshifts 
are known, almost all from optical spectra. This "redshift industry" continuously produces new points for the 3-d 
maps of spatial galaxy distribution. Special telescopes are dedicated to measurements of redshifts. 

Table 1: Some recent galaxy surveys. The columns give: name of survey, solid angle O covered by survey, apparent magnitude limit, 
total number of galaxies N, the method of distance determination. 
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Many extensive redshift surveys have already been completed, among them what are known by the abbreviations 
CfA, SSRS, LCRS, ESP. Their relevant parameters are presented in Tab l2.5.ll For a more detailed review the reader 
may consult Sylos Labini et al. (1998). Based on these surveys several 3-d maps of galaxy distribution have become 
available: both wide angle such as CfAl, CfA2, SSRS1, SSRS2, Perseus-Pisces, LEDA, and narrow angle such as 
LCRS, ESP. 

The last decade saw the appearance of essentially larger surveys with hundreds of thousand of redshifts: the two- 
degree field galaxy redshift survey 2dF and the Sloan Digital Sky Survey SDSS. The depth of these galaxy catalogues 
allows one to detect and analyze structures with sizes up to 100 Mpc (see FigQJ. 

Limits of a survey. A redshift survey is basically restricted by two limits: the apparent magnitude limit of the survey and 
the distance modulus limit which is different for different absolute magnitudes. In addition, it is important for the structure 



Figure 4: The large scale structure from the CfA and the first release of SDSS galaxy rcdshift survey (from Courtois et al. 2004a). One 
may easily recognize different kinds of structure on scales of about 100 Mpc. 

analysis to have a sufficiently large sky coverage. 

The number of galaxies per square degree increases steeply with magnitude. There is roughly one galaxy/deg 2 up to 
m w 15.5 and about 85 up to 19 mag. On the other hand, this may be coped with the modern MOS (multiple object 
spectrograph) techniques which allow a simultaneous measurement of several tens of spectra of galaxies in the field of the 
telescope. Even if the sample is complete up to a certain apparent magnitude limit, its completeness in volume depends on the 
absolute magnitude: for galaxies with absolute magnitude M, the spatial completeness limit of the survey is at the distance 
modulus 

Mlim = miim - M (31) 

For example, with mii m = 15.5 and M — —20, /i = 35.5 which corresponds to the distance r\ = 125 Mpc (corresponding to 
9300/i 75 km/s). On the other hand, galaxies with M = — 18 will be completely sampled only up to riim = 50 Mpc (3700/i 75 
km/s). This shows that redshift surveys probe the distribution of only the most luminous galaxies at large distances, which may 
cause a biased picture. In order to cope with this problem, the only way is to push the surveys to fainter and fainter magnitude 
limits, because one cannot tell beforehand whether a galaxy taken is intrinsically bright or faint (or whether it is distant or 
relatively near). 

2.5.2 Galaxy catalogues with photometric distances. 

For the study of large-scale structure, redshift catalogues are generally superior to those based on photometric distances: 
1) Redshift catalogues are much less time-consuming to make, 2) generally photometric distances are less accurate 
than redshift distances, and 3) redshifts may be measured for all Hubble types of galaxies. 

However, there are questions which necessarily require large samples of galaxies for which both redshifts and 
independent photometric distances have been measured (e.g. the Hubble constant, deviations from the Hubble law, 
peculiar velocities and large scale streams). On the other hand, even photometric distances as such are valuable for 
investigating structures. They can be used for mapping the environment of the Local Group and for deriving the 
average number density law around us. 



As radial velocities of galaxies, especially of those within groups and clusters, give only an approximate estimate of 
distances, Igor Karachentsev started a vast program of distance measurements to nearby galaxies, using the luminosity 
of their brightest blue and red stars (Karachentsev & Tikhonov, 1994; Karachentsev et al., 1997) and the luminosity 
of the tip of the red giant branch stars (Karachentsev et al. 2003). Over the last 10-15 years many nearby galaxies 
have been resolved into stars for the first time. This labour-consuming program requiring a lot of observing time with 
the largest ground-based telescopes, as well as the Hubble Space Telescope, is not yet complete. So far, the distances 
have been measured for about 150 galaxies situated within 9h^Q Mpc (Karachentsev et al. 2003). 

During the last years, a special effort has been undertaken to increase the Local Volume sample. "Blind" surveys of 
the sky in the 21 cm line (Kilborn et al. 2002), infrared and radio surveys of the Zone of Avoidance (Kraan-Korteweg & 
Lahav, 2000) and searches for new dwarf galaxies of very low surface brightness based on the POSS-II and ESO/SERC 
plates (Karachentseva & Karachentsev, 1998, 2000) have made the total number of the Local Volume galaxies increase 
more than two times. 

A formidable example of a galaxy sample with photometric distances up to w 100 Mpc is the KLUN galaxy sample 
and its growing version KLUN+ (see KLUN+ home page at http://klun.obs-nancay.fr/KLUN-l-). These contain 
thousands of spiral galaxies for which both photometric magnitudes and the width of the neutral hydrogen 21cm 
line have been measured. Using the Tully-Fisher relation these two quantities (properly corrected and reduced to 
homogenized systems) allow one to make an estimate of the distance of a galaxy. Originally planned for measuring the 
value of the Hubble constant, the KLUN with its 5500 galaxies was used to study the radial number density distribution 
of galaxies around us (Teerikorpi et al. 1998). The new KLUN+ project is based on a Cosmological Key Project at 
the refurbished Nancay radio telescope in France. The HI survey plans to build a uniquely large magnitude-limited 
sample of 20 000 spiral galaxies distributed on the 80 percent of the sky visible from Nancay (S > —40°). Theureau 
et al. (2004) give information and first results on the new HI measurements. The photometry comes from the DENIS 
(Near Infrared Survey) and 2MASS (2 Micron All Sky Survey). The aim is to build a sample complete to well defined 
magnitude limits in five photometric bands B, I, J, H and K (the earlier KLUN had only B magnitudes and diameters). 

All these data give us the possibility to study the large scale galaxy distribution by means of different statistical 
methods which we describe below. 

2.5.3 How to discover fractal structures 

From the 3-d map shown in fig|U (SDSS) one may recognize by eye structures with different sizes up to 100 Mpc. 
However, a quantitative analysis of the observed inhomogeneities in the spatial distribution of galaxies is a hard 
problem around which there is going on the debate on fractality of observed structures. 

In order to understand the meaning of conclusions from different analyses of the observational data, it is highly 
advisable to investigate idiosyncracies and limitations of the used methods. The next section addresses this issue 
and discusses main mathematical tools used to discover fractal structures. We restrict our consideration to the 
methodological questions arising in practical applications of the methods to analysis of galaxy samples. 

3 Statistical methods to detect fractality in galaxy distribution 

Several statistical methods have been used for the analysis of both 2-d (angular distributions) and 3-d (spatial distri- 
butions) galaxy catalogues. A comprehensive review of mathematical approaches for the description of the large scale 



galaxy distribution may be found in the book by Martinez & Saar (2002). 

Here we confine our discussion to the fractal approach. Detailed descriptions of such methods for the analysis of 
the distribution of galaxies are given by Sylos Labini et al. (1998) and Gabrielli et al. (2004). 

One of our goals is to compare the standard method of the £ correlation function (Peebles, 1980; 1993) with 
the fractal approach (Pietronero 1987; Gabrielli et al. 2004). We demonstrate that it is essential for the study of 
the distribution of galaxies to utilize mathematical instruments which are adequate for the existing structures. It is 
especially important to use undistorted estimators for measuring the fractal dimension and the range of fractality. 

3.1 Definitions for correlation functions 

The theory of stochastic processes introduces and studies different functions intended for the correlation analysis (see 
e.g. sect. 2 of Gabrielli at al. 2004). 

3.1.1 Complete and reduced correlation functions 

The complete two-point correlation function (we call it simply the complete correlation function) of a stationary 
isotropic process fj,(f) is defined as 

R^(r) = Ou(riMr 2 )) (32) 

where r =| r \ — \ rj — r% | is the mutual distance between considered points, and (•) is the ensemble average over all 
realizations of the stochastic process. 

Taking into account the truly constant mean value of the process 

(j.0 = (M^)) = cons t (33) 

one may define the reduced two-point correlation function Ci for the fluctuations around /io (we call it simply the 
reduced correlation function) as 

CM = ((Kn) - Mo)(/x(f 2 ) - M o)) = R^(r) - Mo • (34) 

For r = it expresses the squared dispersion of the process as cr^ t = 6*2(0). 

We emphasize an important difference between the complete and reduced correlation functions i?^ M (r) and C^r). 
For a stochastic process with long range power-law correlations the complete correlation function R^ (r) has a power- 
law form, while the reduced correlation function 62(7"), according to its definition feci I34JI . cannot be a power law in 
this case. 

3.1.2 Mass variation in spheres and characteristic scales 

In the applications that we discuss the stochastic process (x{r) will describe the density field p(r). 

Conditional and unconditional functions. It is important to distinct between conditional and unconditional 
functions (or statistics). For instance, when one considers such statistics which are defined with the condition that 
there is a fixed point-particle relative to which other particles of a process are considered, then one speaks about a 
conditional function. We will see below that the two major tools of the LSS analysis, the T and £ functions, are both 
conditional correlation functions. 



As an example of unconditional statistics we consider mass (number) fluctuations inside a sphere of radius R. Let 
us define in addition to p(r) a new stochastic variable M(R) as 

M{R) = { p(r)d 3 r (35) 

JV(K) 

For a given radius R fluctuations of this mass calculated in different positions in space can be characterized by the 
normalized mass variance o~ 2 M {R): 

, t ^_{M{Rf)-{M{R)f 



where 



um{r) = — ww 2 — 5 (36) 
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(M(R)) = -jpoR 3 (37) 



and for volume V = V(R) 

(M(R) 2 ) = [ d 3 n [ d 3 r 2 (p(n)p(r 2 )) (38) 
Jv Jv 

Here there is no condition on the locations of the centre of the sphere, which may be put anywhere in the space within 
the sample regardless of the positions of the particles, also "between" them. 

The scale of homogeneity. Our intuitive vision of uniformity may be formalized by means of the variable M(R). 
E.g. the homogeneity scale i?hom may be defined as the scale at which afj (i?hom) = 1 (or some other threshold value). 
This means that it is possible to regard a distribution of particles as approaching homogeneity if the average mass 
fluctuation within spheres of radius i?hom is about the average mass (Af(i?hom))- This scale is well defined when 
a 2 M (R) —> for R > R h 

om ■ 

The correlation length. The second scale is the correlation length R COI , which does not depend on the amplitude of 
the correlation function and just characterizes the rate of decrease of the correlation function. The correlation length 
may be infinite, as it is for a power law correlation i? MAI (r) oc r -7 , or finite, as for an exponential correlation function 
R^(r) oc e-'V-R- 

3.2 The method of the ^-correlation function 

A widely used classical approach to the analysis of the large scale structure is the method of correlation function. It 
was first introduced to the galaxy analysis by Totsuji & Kihara (1969) who adopted this method from the statistical 
physics of ordinary gas density fluctuations (e.g. Landau & Lifshitz 1958). It was further developed and extensively 
applied to galaxy data by Peebles (1980, 1993) and others (for recent reviews see Martinez & Saar 2002; Jones et al. 
2004). 

3.2.1 Peebles' ^-correlation function 

According to Peebles (1980) the two-point correlation function £(r) is defined as the dimensionless reduced correlation 
function of the density fluctuations Sp(r) = p(r) — po around the average density po 

_ ($P(ri)5p{r 2 )) _ (p(fi)p(r 2 )) - p% ^ 
Po Po 

In fact the ^-function is simply the reduced correlation function fea !34|l divided by the squared mean value of the 
process (eqEHJl, i.e. 

Of) = C 2 (r)/p (40) 



In the case of a distribution of identical particles (with masses m = mo) one uses a number density n(r) = p(r)/mo, 
whose average is (n(r)) — no- Then 

((r)= (n(rWi + r1)_ L 

This dimensionless function measures correlations of fluctuations relative to a constant average number density n . 

In the theory of stochastic processes one usually considers a normalized correlation function which is defined as 
Kfj,fj,(r) = C2(r)/(Jx — — xfy/a^.. Then there is the normalization condition K^ifi) = 1. The definition for 

the ^-function feci 139(1 implies the condition £(0) = o^/po- 

Definition via Poisson process. The correlation function may also be defined as a measure of the deflection 
of a distribution of particles from a Poisson (uniform) distribution (Peebles 1980). In this case one considers two 
infinitesimal spheres at the points r% and f*2 with volumes dV\ and dV 2 and with the mutual distance r\i. Then the 
joint probability to find one particle in the volume dV\ and another particle in the volume dV2 is proportional to the 
number of pairs dN\ 2 

dN 12 = nldVxdVzil + £(r 12 )] , (42) 

where no is the ensemble average number density of particles and £(^12) measures the deflection from the Poisson 
distribution. This definition implies that £(^12) = automatically for a Poisson process. 

For a statistically isotropic distribution the function £(^12) = £(r) depends on the separation r only. For the case 
when an object is chosen at random from the sample, the probability of finding that it has a neighbour at a distance 
r in dV (e.g. Amr 2 dr) is proportional to the expected number of neighbours dN 

dN = n dV[l + £(r)] . (43) 

Here £(r) is considered as the same two-point correlation function as defined by eq^J (see sect. 33 of Peebles 1980). 
It is a measure of finding an excess number of particles relative to the Poisson distribution, at the distance r > 
provided that there is a particle at r — 0. 

3.2.2 ^-function estimators 

In the theory of stochastic processes it is important to make a distinction between functions (e.g. £(r)) defined by 
ensemble averages and their estimators e.g. applied to a finite galaxy sample. 

To estimate the two-point correlation function from an available data sample of N s objects within a volume V s , 
one generally uses a method based on artificially generated Poisson process, which fills the same volume V s of the 
sample. Then the ^-function for a given scale r is estimated as the ratio of the number of pairs with such mutual 
distance in the sample to the number of such pairs in the artificial Poisson distribution. There are several different 
pairwise estimators (Kerscher et al. 2000; Martinez & Saar 2002; Gabrielli et al. 2004), and the difference between 
them lies mainly in their method of edge correction. 

For instance, the Davis-Peebles estimator weights the points according to the part of the spherical shell volume 
which is contained in the volume of the sample. It has the form 

Here N p (r) is the number of data-data pairs of observed objects in the catalogue having their mutual distance in the 
interval (r, r + dr). N Pj( i r (r) conies from the joint catalogue of data and artificial random distributions in the same 



volume V s . It is the number of data-random pairs with the distance r in the joint catalogue. N s and N r d are the total 

numbers of objects in the real sample and the random distribution, respectively. 

An essential assumption of the correlation function method is the hypothesis of homogeneity according to which 

the true average of objects no is estimated from the observed sample as 

N a , . 

n Q =n = — (45) 

with a high formal accuracy <7„ m l/s/N s , where N s is the total number of objects in the volume V s of the "fair" 
sample, which is assumed to be representative of the homogeneous distribution of galaxies in the whole Universe. 

3.2.3 The normalization condition for £ estimators 

A significant point related to ^-function estimators was emphasized by Pietronero (1987) and Calzetti et al. (1988), 
and in more detail by Gabrielli et al. (2004). Namely, the definition of the correlation function as a deflection from the 
Poisson distribution fea !43|l implies an integral condition for the £ function estimated from a finite sample of galaxies. 

This comes from the fact that for any sample with a finite number of galaxies N s in a volume V s the estimation of 
the average number density is n = N s /V s . Integrating the left side of ea 1431 over the sample volume we get 

dN = N s - 1 (46) 

v 3 

where N s — I is the number of neighbours, i.e. the total number of particles in the volume V s without the one whose 
neighbours are counted. Then the integration of the right side of ea 1431 over the sample volume gives 



N s — 1 = / ndV + n / £(r)dV. (47) 

Jv s JVs 

The first term on the right side is J v ndV = N s , the total number of the particles in the sample. Hence the second 
term will satisfy the condition 



/ 



£(r)dV = -- (48) 
iv s n 

or in dimensionless form: 

In the case of fluid-like correlated distributions the effective number density of particles may be arbitrarily large and 
hence the condition of ea 1151 becomes 

/ i(r)dV = (50) 
JVs 

These restrictions lie behind some controversial results obtained by the £ function method of the large scale structure 
analysis. In particular, we have in mind the inevitable non-power law behaviour of the £ estimator. From eas !48l and 
1501 follows that there is a distance r z where £(r z ) — 0. Here the correlation function estimator changes its sign from 
positive to negative values, which is impossible for a power-law function. 

3.2.4 A systematic distortion of the true power-law correlation due to the ^-estimator 

As was shown above, if the complete correlation function is a power-law then neither £(r) nor £(r) can be a power-law 
function. Nevertheless, in practice £(r) is usually presented in the form 

|(r) = , n < r < r 2 , (51) 



valid for some range of scales r\ < r < r?,- Here the parameter ro defines the amplitude. 

From such a power-law presentation one usually derives two numbers: the unit scale ro and the correlation exponent 
7. We emphasize that due to the normalization condition fea l49fl both numbers give systematically distorted values 
for the homogeneity scale and the power-law exponent of the true complete correlation function describing the 

density field. 

The unit scale ro (which is often called, somewhat misleadingly, as correlation length) is defined from the relation 

£>o) = 1 , (52) 

which characterizes the amplitude of density fluctuations at the scale ro. In fact, it is a distorted value of a true 
homogeneity scale of the process if the true value is equal to or larger than the size of a maximum sphere embedded 
completely in the sample volume V s . 

The correlation exponent 7 in the power-law representation of £(r) (eal51[l describes correctly only a restricted 
interval of scales r± < r < T2 . On scales r > ri this does not represent the true value of the exponent, because there 
the estimated value is distorted as the normalization condition fea l49J) makes £(r) deflect from the inherent power-law 
and to cross zero level. For example, below it will be shown that for the exponent 7 is derived twice the true correlation 
exponent at the unit scale ro (sect. 3.4.1). 

On scales r < T\ the true value of the exponent is distorted due to the noise of discreteness, which behaves as 
a oc r~ 3 / 2 fea ll9J) . The error will essentially increase for scales smaller than the average distance between particles in 
a sample (e.g. r < Ao). We will see later examples of how this has happened in data analysis. 

3.2.5 Redshift-space and the peculiar velocity field 

From a galaxy redshift survey one obtains a redshift-space map, i.e. (a , 5, z bs) coordinates, where (a , S) give the 
position on the sky and z b s gives the observed redshift of a galaxy in the sample. The value of z bs in principle 
contains all possible contributions from different physical causes according to the relation 

1 + Z ohs = (1 + Z cos )(l + Z v )(l + Z<j,)(l + Zncw) (53) 

Here z cos is the cosmological redshift which determines the true distance to the galaxy r ga i(z) calculated from the 
empirical distance-redshift relation (i.e. from the Hubble law, which gives r ga i = cz cos /Ho) or from an adopted 
cosmological model for large redshifts. The z v is the redshift component caused by the peculiar velocity of the galaxy, 
Z0 is the gravitational part of the redshift caused by the local gravitational potential of the galaxy (e.g. in a cluster), 
and z n ew is a possible component due to unknown cosmological physics. 

Distance error due to peculiar velocity. Let us consider the influence of the peculiar velocity field on the distance 
estimation. For peculiar velocities v << c the Doppler part of the observed redshift is determined by the radial 
component v r of the velocity as 

z v sa v r /c . (54) 

So for small radial peculiar velocities v r we obtain 

Zobs = z cos H (1 + ^cos) (55) 

c 

Hence the true spatial galaxy distribution will be distorted by a peculiar velocity field in the line of sight direction 

Hlobs = r || + Ar„ (56) 



by the value 

_ u r (1 + z cos ) _ _ x v r (l + z cos ) 

Ar "- Ho - jM P cft ioo 500km/s t 57 ) 

Note that the factor (1 + z cos ) leads to an increasing influence of v pcc on the distance distortion for deep redshift 

surveys. 

Real-space and redshift- space £ functions. A directly observed ^-function is called the redshift- space correlation 
function £ z (s)- In order to obtain the real-space correlation function £ rca i = £(r) one should extract and delete all 
non-cosmological contributions to z b s - This is a hard problem because it requires a priori knowledge of the peculiar 
velocity field, the total mass around the galaxy, and restrictions on new physics. 

The shape of the observed £ z (s) is determined by the character of the peculiar velocity field. In virialized clusters 
the velocity dispersion leads to the so-called "fmgers-of-God" effect, i.e. an elongated shape along the line of sight 
direction ru . The mean tendency of galaxies at larger scales to approach each other due to the gravity of large-scale 
structures will appear as a compression of £ 2 (s) in the direction ru. As these two effects are related to different spatial 
scales, they do not compensate each other. 

Peebles (1980; sec. 76) and Davis & Peebles (1983) suggested a procedure for the restoration of both the true £(r) 
and the relative peculiar velocity distribution f{v) from the observed correlation function £ 2 (rx, r|| obs ) where rj_ an 
r ||obs are the observed perpendicular and parallel to the line of sight components of the separation s = yV^ + ^jj obs - 

Derivation of the real-space ^-function. The method is based on the calculation of the projected correlation 
function w(r±) which does not depend on the peculiar velocity field, if the distribution of the radial peculiar velocities 
is symmetrical around each galaxy of the sample. Then integrating along the line of sight we obtain: 



w(rj_) = 2 £z(r±,ri\ obs )dri\ obs , (58) 
Jo 

where in practice the interval of integration is restricted by chosen radial velocity limits. Then the wanted inverse is 
the Abel integral 

1 f°° w'dr ± , . 

where w' = dw(r±)/drj_. Eg 1591 gives the solution for the problem of restoration of the real-space £(r) correlation 
function. 

For the power law £ = (r /r) 7 the integral ifB"^ gives 

w(r ± )=Ar 1 ± 1 , (60) 

where 

7 r e (i/2)r e (( 7 -i)/2) 
A ~ r ° r„( 7 /2) (61) 

and r e (ir) is the Euler gamma function. 

Limitations of the projection method. It is clear from ea !6UI that such a solution for the real-space £ correlation 
function is valid only if the exponent 7 > 1. For example, 7=1 gives w(r±) = constant, which demonstrates that a 
uniform background galaxy distribution may be confused with the projection of real-space non-uniform distribution. 

According to the theorem on fractal projections (sec. 2. 4. 5) such a method inevitably leads to elimination of in- 
formation on structures with the fractal dimension D > 2. Therefore to take into account the peculiar velocity field 
within fractal structures with D > 2 (such a structure is actually observed, see sec. 5), it is necessary to use another 



method of restoration for the correlation function, which is free from the above limitation. Also a more careful study 
of distance errors in both r\\ and rj_ components is required. 

Estimation of the relative velocity dispersion. In the case, when density and velocity fields are weakly coupled, 
the observed correlation function £ 2 (s) can be modelled as a convolution of the real space correlation function £(r) 
with the galaxy pairwise velocity distribution f(v). According to Peebles (1980, sec. 76) and Davis & Peebles (1983) 
this equation may be presented in the form 



1 + e»(r±,r|| obB ) = Ho J [1 + t(y/f± + rjpj f(v)dr {l (62) 

where 

v = H r l]obs - Hor^ + v 12 (r) (63) 

and U12O") is the mean radial pairwise velocity of galaxies at separation r, which is represented by a model. Davis & 
Peebles (1983) adopted the model 

*»w = Twh (64) 

and an exponential form for /(«): 

f(v) = B exp (-2 1 / 2 I v I /<ti 3 ) (65) 
As a result of this approach one obtains the pairwise velocity dispersion 012. 

3.3 The method of conditional density 
3.3.1 Definitions 

The method of conditional density has been successfully used for analysis of fractal structures in modern statistical 
physics. This method was proposed for extragalactic astronomy by Pietronero (1987) and has been applied to 3-d 
galaxy catalogues by many authors (see reviews in Coleman, Pietronero 1992; Sylos Labini, Montuori, Pietronero 
1998; Gabrielli et al. 2004). Conditional density method has the advantage that it gives an undistorted estimation of 
the true power law correlation and the true fractal dimension. It also may be used for finding an undistorted value of 
the homogeneity scale of a galaxy sample. 

Continuous stochastic processes. The conditional density T(r) may be defined by means of the complete correlation 
function fea l32|) in the following form: 

r(r) = ^ w(r) = MfM?l±?0) ( 66 ) 
Mo Po 

Here p(f) is the stochastic density field and po is the ensemble average density. The r-function has the physical 
dimension of density [g/cm 3 ], and it is a measure of correlation in the total density field without subtraction of the 
average density. The physical dimension of the T function agrees with the common interpretation of T(r) as an average 
density law around each point of the structure. This makes its estimator a natural detector of fractality. 

As we shall see below this definition allows one to construct such an estimator which has no additional restrictions 
like the normalization in eg 1481 and hence is able to give an undistorted value of the exponent for true power-law 
correlations. 



Discrete stochastic fractal processes. Let us consider a discrete stochastic process, one realization of which is a 
set of identical particles at randomly selected positions {x a }, a = X, ■■■,N, so that the number density n(x) is given 
by the expression 



N 

n(x) = y^ S(x- x a ) . (67) 

a=l 

If the stochastic process generates a fractal, then it is natural to define the number density as a function of two 
variables: n = n(x, r). The first variable describes the position x a of a structure particle, and the second variable gives 
the radius r of a ball inside which one calculates the number of particles of the structure. The variable r serves for 
constructing a statistics which can measure the strength of the singularity around a particle of the fractal structure 
where the number of particles grows as a power-law N(r) oc r D . 

Denote by Nv(x a ,r) the number of particles in a sphere of radius r, centered at the particle a with the coordinates 
x a , belonging to the structure: 



Nv(x a ,r) — / n(x~a + x)Attx dx , (68) 
Jo 

and Ns(x a , t) is the number of particles in the spherical shell r, r + Ar, with the centre at x a : 



/•r+Ar 

Ns(x a ,r) = / n(xT a + x)4:TTX 2 dx . (69) 

J r 

From one realization to another these quantities fluctuate, but after averaging over many realizations the stable 
power-law dependence on the scale r emerges. In the case of ergodic processes averaging over many realizations may 
be replaced by many points in one realization. Following the work by Pietronero (1987) we define the conditional 
(number) density of a stochastic fractal process in the form: 

and the conditional volume density as 



r*(r) = ( 7;7°'V ) =^r-^- D \ (71) 



Ny(x a ,r) \ = 3B r _ {3 _ D 
(47r/3)r 3 / Sg 4tt 

where (-)s a means averaging over all points x a in one realization with the condition that the centres of the spheres lie 
at the particles of a realization (this explains the word "conditional" ) . The last equalities in l|7U|l and (|T1|) relate to 
ideal fractal structures, for which 

r*(r) = |r(r) (72) 

3.3.2 T-function estimator. 

Consider a stochastic fractal process where the number of particles Ny(x a , r) in a sphere of radius r, centered at the 
point x a and the number Ns{x a , r) of particles in the shell (r, r + Ar) are given by eas !68l and l69l Taking into account 
definitions of conditional densities (eqs. I7UI and I7T)) one can use following two statistics for their estimation from one 
realization (a finite galaxy sample): 



I N 2. rr+Ar 

f (r) = N ^ 4^ / n(s « + ^ xHx (?3) 
for the shell conditional density T, and 



N 

(74) 



1 3 /"' 

a=l J ° 



for the volume conditional density T*. So the use of the conditional density method is in principle quite simple, just 
counting the number of particles inside the spherical volume V(r) or inside the shell S(r)Ar. This is done for each 
structure point and then the average is calculated. For the T-function estimation one need not generate artificial 
Poisson distributions, which was necessary for the ^-function method. 

Fractal dimension and co-dimension. For a fractal structure both the T function fea l7L)fl and the estimator (ea l73|l 
have a power-law form 

f(r) = T r-^ (75) 

This very important property of the T-estimator allows one to obtain an undistorted value of the fractal dimension in 
a galaxy sample. 

The exponent that defines the decay of the conditional density 

7 = D - 3 (76) 

is called the co- dimension, where D is the fractal dimension (or the correlation dimension D 2 ). The amplitude r 
of the estimator T(r) does not change when the sample volume V s is increased, only the range of available scales r 
increases. This corresponds to the meaning of T(r) as characterizing the number density behaviour. 

Homogeneity scale. For a fractal structure which has an upper cutoff at a homogeneity scale -Rhom, after which 
the distribution becomes uniform, the fractal dimension D — 3, and the estimator of the T-function is 

T(r) = constant, for r > Rhom ■ (77) 

Thus the method of conditional density is a powerful instrument when one searches for the crossover from the regime 
of fractal clustering to the realm of homogeneity. 

So, for processes with a finite upper cutoff scale of fractality, beyond which the distribution of particles turns into 
homogeneity, the statistics l|73|) and (|74|l give constant values corresponding to the fractal dimension D = 3 as it 
indeed should be for uniform structures. 

3.3.3 Redshift-space T z (s) and v pec 

As we discussed in sec. 3. 2. 5 an estimation of the spatial distribution of galaxies from redshift catalogues is based on 
the redshift-space (a , 5, z t, s ). Using the same notations as in sec. 3. 2. 5 for the £ function (s, r±, r||) we can right for 
the redshift-space conditional density 

r z (s) =r 2 (r ± ,r,| obs ) . (78) 
The relation between the real-space r(r) and redshift-space T z (s) conditional density is 

r,(r ± ,r|| obs ) = f g(f, w)T(r)d 3 w , (79) 



where w = V\ — V2 is the relative peculiar velocity of a galaxy pair at separation r, and g(r, w) is the relative peculiar 
velocity distribution. Here the components of the relative distance r are given by the following formulae: r\ = 
^3 = r|| obs - w 3 /H Q , and r 2 = r\ + (r|| obs - w 3 /H ) 2 . 

In order to restore the real-space conditional density from the directly observed redshift-space conditional density 
it is necessary to make computer simulations of artificial fractal structures with known peculiar velocity fields and 
then compare the modelled redshift-space r mo( j(s) with the observed T z (s). 

3.3.4 T-function for 2-d intersections 

If in 3-d space a fractal structure is intersected by a plane then the expected value of the fractal dimension for the 
intersection is given by ea 1301 

Ant = D - 1 (80) 

To make the T-function analysis for the sample which presented the 2-d intersection we shall use the 2-d coordinate 
system y — (yi, j/a) for which we can calculate the T- function for the intersection ri nt (y). 

Such a situation may occur in a slice-like galaxy survey for scales r larger than the thickness of the survey. E.g. 
for the true fractal dimension D = 2 the fractal dimension of the intersection will be D lnt = 1. Hence we expect to 
obtain a power-law behaviour of the corresponding ri nt (y) cx y~ l ■ We will see below that the intersection theorem 
helps one to understand the behaviour of the power-spectrum derived from slice-like galaxy surveys. 

3.4 Comparison of correlation function and conditional density 
3.4.1 The relation between T and £. 

From the definitions of the conditional density V (ea 166(1 and correlation function (ea 141(1 we have the relation: 

T(r) 

f(r = - 1 (81) 
n 

if the average number density no of a considered stochastic process exists. Both T and £ functions are conditional 
characteristics of a stochastic process, i.e. they are defined on the condition that the centres of counting spheres are 
set to structure particles. However, there is still a deep difference between them. The T(r) represents the complete 
correlation function, while the £(r) represents the reduced correlation function of the stochastic process. This fact 
makes the properties of the corresponding £ and T estimators very different. Finally, it results in conflicting values for 
the estimated correlation exponent and homogeneity scale of a galaxy sample. 

From eqOHl follows a similar relation between the estimators for a finite sample of galaxies: 

e>) = ^-i (82) 

n 

Here £fs is called the "full shell" estimator because it is defined through the T estimator which is calculated using 
full shells completely embedded in the sample volume. 

The estimator T (ea 173(1 is always a positive function and has a power-law form for fractal structures. On the 
contrary, the estimator £ fca l44|) inevitably changes its sign and hence cannot be presented as a power law even for 
scale invariant structures. All estimators of the ^-function, which are based on counting of pairs relative to an artificial 
Poisson distribution, have a common drawback. They give essentially distorted values for the true correlation exponent 



of the complete correlation function of long range power-law correlated processes and for fractal distributions. But 
the T-function estimator is specially constructed in order to give undistorted values of the correlation exponent and 
fractal dimension. 

The T-function estimator relates to intrinsic properties of the sample, while the ^-function estimator depends on 
both intrinsic and external properties of the sample. In fact, T measures the behaviour of the total density inside 
spheres within a sample, while £ measures density fluctuations relative to the average density, which is assumed to 
be valid for all space outside a finite sample. This can be illustrated also by the following reasoning. Let us consider 
counts around a fixed point. The expected number of points in a shell with radius r and volume dV is n(r)dV , where 
n(r) is the conditional density T(r) describing the density-radius law. On the other hand, the same expected number 
may be calculated with the correlation function £(r) as no(l + £(r))dV. So 

r(r) = (l + £(r))no (83) 

It is important to note that the right-hand side of eq|H3l becomes defined only after the mean density no is calculated 
for the whole sample, while T(r) always exists locally. Remember that £(r) was defined for fluctuations around the 
mean hq. 

Fractal density field. For the case of a scale-invariant stochastic fractal density field, the complete correlation 
function T(r) has the power-law form 

r(r) = ™ r-f , (84) 
while for the same fractal structure the reduced correlation function £(r) will be 

£(r) = J?" r~ 7 " 1 (85) 

which is not a power-law. This difference between complete and reduced correlation functions was pointed out by 
Pietronero & Kuper (1986). 

Thus the ^-function may be approximated by a power-law only for such r when £(r) 3> 1, which corresponds to 
small scales r << ro- However, on small scales the noise of discreteness is essential, which also leads to distortion of a 
true power-law. Hence a ^-function estimation gives a distorted values of the correlation exponent not only on large 
scales (normalization), but also on small scales (discreteness). 

It is instructive to calculate the exponent 7* of the correlation function (eg 185( 1 on scales close to the unit scale r . 
Taking the logarithmic derivative of ea!85l one obtains (Joyce, Montuori, Sylos Labini 1999): 

(r)= _feM = Wren (86) 
751 ' dlogr 2(r/ro)-7-l 1 ' 

Therefore at the unit scale r — ro we get the remarkable result: 

7«M = 2 7 (87) 



For example, for a true density power-law with 7=1, one would infer an apparent slope 7^ = 2 for the correlation 
£- function, if measured at scales close to the "correlation length" ro! 



3.4.2 The dependence of ro on sample depth. 

The sample depth R s is an important global parameter of an observed galaxy distribution. For a fractal structure 
sampled inside a spherical volume V s — ^-R 3 with h = N s /V s = 3BR®~ 3 /Air, the eal551 yields 

<M-?(ir->. 

and the unit scale ro comes to depend on the sample parameters. 

Inserting r = ro into ea 18 II and taking into account that £(r = ro) = 1 one gets 



l/ 7 



where n is the average number density of objects in the sample. Hence one obtains: 

r(ro) = IF = m = S^" (90) 

From this follows a simple relation between the correlation length ro and the depth of the sample R s : 

r ° = {^ 1 ) lh Rs (91) 

We note again that this is true for a spherical sample (i.e. R s — R^ax) an< ^ a nxe d luminosity (i.e. L — const) of 
sample galaxies. 

3.4.3 Geometry of a survey and characteristic scales. 

A strong restriction for the practical application of the T-function method is the requirement that there should be 
room for the whole sphere in the volume of the considered sample. For example, for galaxy surveys with slice-like 
geometry, this makes it impossible to measure the conditional density for scales larger than the thickness of the slice, 
i.e. the diameter of the maximum sphere fully contained in the survey volume. 

The above dependencies between ro and various sample parameters were derived for the ideal case of a simple 
fractal in a sufficiently large spherical volume inside which the correlation function is estimated. For a non-spherical 
survey geometry (like a slice or cone) these relations are valid only for such scales r for which the survey galaxies 
are contained completely within the sphere with the radius r. For non-spherical geometries these dependencies are 
expected to differ from the above predictions. 

For a galaxy sample under consideration one should always control the following characteristic distance scales: 

Rsep, Rmax> ■ (92) 

The separation distance R sep between galaxies in a sample may be roughly estimated as n -1 / 3 or calculated from 
the nearest neighbours distribution. We may define R sep — aR sep - av , where a rs 0.1 and R sep - av is the average 
distance between nearest neighbours in a sample. For r < R sep the discreteness noise is important. 

The radius of the maximum sphere i!^ is related to the complete embedded sphere in the sample and play a 
crucial role in estimation of correlation properties of the galaxies from the sample. The depth R s of a survey is related 
to the largest distances of galaxies in the sample, it essentially differs from the radius of the maximal sphere in the 
case of a slice-like surveys. 



As a general rule for a correlation analysis of slice-like samples one should consider separately two regions of scales. 
First, 

R sep < r < RZ h ax , (93) 
where it is possible to use the T function method and estimate the value of the true fractal dimension, and secondly, 

JCi < r < R s , (94) 

where the T function cannot be applied. In this region other methods should be developed such as the T function for 
fractal intersections and the two-point conditional column density. 
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Figure 5: Probability distributions for finding the nearest neighbour at a distance r for two cases: (top) the homogeneous Poisson 
distribution (D = 3) and (bottom) a fractal distribution with D = 2 (a random Cantor set). The distributions are in 3-d space inside a 
sphere of diameter = 1, the total number of points in each case is N = 2.5 ■ 10 3 . The dotted curves correspond to theoretical predictions 
according to eqs. 99, 100 (from Vasiliev 2004). 



3.5 Nearest neighbours distribution. 

The distribution of distances to the nearest neighbour point is a useful property of fractals, which may be used to 
make a distinction between fractal and ordinary distributions. It is discussed in detail by Gabrielli et al. (2004). 

Poisson distribution. For example, for a Poisson process in 3-d space the probability density distribution uu{r) for 
finding the nearest neighbour at a distance r is 



oj{r) = iTTnor 2 exp 



(95) 



— 1/ 3 

For Poisson process the average distance between point-particles is R se p ~ n D ' 

Fractal distribution. There is no exact formulae for the general case of a fractal structure. However, Gabrielli et 
al. (2004) derived a useful approximation for the probability density distribution for the nearest neighbour in a fractal 
structure: 

w(r) =4^Cr I5 - 1 exp[ , (96) 

where C = DB/4ir. 

Figure [3 presents to(r) for two cases: a homogeneous Poisson distribution (D = 3) and a fractal distribution with 
D = 2. The dotted curves correspond to theoretical predictions from eqEHJand eq|^Hl It is seen that the lower is the 
value of the fractal dimension D, the closer are the neighbours within the fractal structure. This example shows that 
the nearest neighbour probability function may be used as an additional method for estimating the fractal dimension 
of a galaxy sample. 

3.6 Two-point conditional column density 
3.6.1 Definitions 

The conditional densities of stochastic fractal processes above discussed were one-point, as the centre of the sphere 
in which particles are counted lies at one fixed point {a} with the coordinates {x a } and the counts are made around 
each point {a} of the sample. However, in some cosmological studies, e.g. related to gravitational lensing (Baryshev 
& Ezova 1997), it becomes necessary to use two-point conditional densities, whereby one simultaneously fixes two 
particles {a, b} with the coordinates {x a , Xb} and counts galaxies within a thin cylinder between these points. 

In order to define the distribution of particles along such a cylinder, whose axis connects two structure points 
{a, b} C {xi , i = 1, ...,N}, Baryshev & Bukhmastova (2004) introduced the concept of 2-point conditional column 
density r) a b{f) for a stochastic fractal process. In applications to galaxies, according to the cosmological principle of 
Mandelbrot the particles a and b are statistically equivalent, hence for each of them the 1-point conditional density is 
given by the expression (|70|l . which is proportional to the probability of finding particles at the distance r from fixed 
structure points. 

Let us now take two independent points of the structure at the distance r a b = \x a — Xb\ from each other. Denote 
by C the event, when particles appear at the distance r a from a and independently at the distance rb from b. Then 
C is given by the union C = A (J B of the two events related to a and b. So the 2-point conditional column density, 
which is proportional to the probability of finding particles around the particles a and b, may be presented as a sum of 
1-point conditional densities. The assumption of independent events is a first approximation. Numerical simulations 
have shown that this approximation is sufficient for the analysis of usual fractal structures (Baryshev & Bukhmastova 
2004). Then the 2-point conditional column density can be expressed as 

Vab(r) = \ [T a (r) + r b (r ab - r)} 



DB _ 7 

r ab 



r \ l r 
— + 1 

Tab J \ Tab 



(97) 



8tt 

where 7 = D — 3. The distance r is measured along the segment of the line connecting the particles a and b, and 
at the same time it defines the radius r of the sphere having its centre at the first point and the radius r a b — r for 



another sphere having its centre at the second point. The constant B defines the normalization. We note that in this 
formula the volume elements are taken along the line connecting the two points - in this sense the coordinate r is 
not a spherical radial coordinate, but a Cartesian coordinate labelling volume elements ("tablets") with thickness dr 
along this cylinder. 

i— • j 




Fi gure 6: (left) The observed distribution of the galaxy two-point conditional column density fj a b{ r ) of the LEDA sample along cylinders 
with I < 100 Mpc. (right) The same distribution in logarithmic scale and the fitted theoretical curve »7 a (,(r), with the fractal dimension 
D = 2.02 (from Baryshev and Bukhmastova 2004). 



3.6.2 Estimation 

In order to estimate T]ab(r) one may use the statistic 

' N c (x a ,x b , r, h, Ar) 
7r/i 2 Ar 



Vab(r) 

{a.b} 



^ N a b ^ nr+Ar rh 

= ^£^/ l<^ hdhd r, (98) 

{a.b} 

where N c is the number of particles in a volume element of the cylinder with a diameter h and height Ar and whose 
axis connects the structure particles a and b. The volume element lies at the distance r from a, which corresponds to 
the distance r a b — r from b. Averaging is performed for every pair of particles with connecting cylinders having the 
length in the interval (I, I + Al). The parameters h and Ar are chosen to be a fraction of the mean separation of the 
particles in the sample. Practical modelling of artificial fractals show that the estimated fractal dimension is robust 
to reasonable variations of the tablet size (Baryshev & Bukhmastova 2004; Vasiliev 2004). 

A more general situation, which may be encountered in practice, is that both a fractal structure and a homogeneous 
background exist in the sample. Then the practical calculation of the fractal dimension D by means of the estimator 
r] a b(r) may be done by fitting through the observations the dimensionless probability distribution with three free 
parameters 7, R\, R 2 



where N(y) is the observed number of points which are found in each tablet, i.e. within the small intervals (y, y + Ay) 
along the cylinder with a length I. The variable y is the relative distance measured along the line connecting the two 
points (y = r/r a b = r/l). N is the total number of points within cylinders of length /. In the case of genuine self- 
similar structures it is possible to calculate these numbers for all cylinders of different length simultaneously. R 2 shows 
the contribution of a possible Poisson background. 

It is useful to introduce a measure of the relative strength j3 of the fractal component as was done by Vasiliev 
(2004): 

fl-i^s. doc, 

When (3 = 1, the contribution from a fractal structure is equal to the Poisson background contribution (R 2 = 0.5). 

The main advantage of the cylinder method is that it can work for slice-like surveys allowing the study of scales 
comparable to the depth of the survey. 

3.7 Fourier analysis of the galaxy distribution 
3.7.1 Ordinary density fields 

If the spatial distribution of objects is given by a stochastic density field, then fluctuations of this field may be 
represented by the Fourier integral as a superposition of plane spatial waves with the wave number k — \k\ — 2tt/X: 

m = ew-e = _1_ j me ^ d 3 k (101) 

Here the Fourier transform S(k) of the density fluctuations S(x) is a complex quantity and thus may be given in the 
form 

6(k) = \6(k)\exp(i<p(k)). (102) 

Here we see that for a complete description of the spatial distribution of objects the analysis should include both the 
amplitude spectrum \S(k)\ and the phase spectrum <p{k). 

In the case of a Gaussian random process the phases of plane waves are distributed uniformly in the interval [0, 2tt] 
and to describe the density field it is sufficient to consider only the power spectrum 

P(k) = (\~5(k)f) . (103) 
When the distribution is isotropic, the power spectrum and the correlation function arc connected by the relation 



(104) 



Thus for a power law correlation function £(r) oc r~ 7 the power spectrum has also a power law form P{k) oc fc 7 ~ 3 for 
restricted range of scales k\ < k < k 2 . For small k (large scales A > R s ) there is a limit due to the size of a survey 
R S} which give P(k) — > for A — > R s . 

The phase spectrum 4>(k) carries important information, complementary to P(k). Possible non-Gaussianity of the 
clustering properties i.e. non-uniformity of the phase spectrum was considered in Chiang & Coles 2000 and Coles & 
Chang 2000. 



The main problem of the power spectrum analysis is the same as for the correlation function: the average density 
p should be well defined within the sample volume. 

However an advantage of the P(fc)-method is that the correlation exponent 7 may be estimated without distortion 
from P{k) for the scales A < R^x- 

3.7.2 Fractal density fields 

For stochastic fractal processes one may perform the Fourier-analysis of fractal density field ( Sylos-Labini & Amendola 
1996, Sylos Labini et al. 1998). They introduced the concept of generalized power spectrum II(fc) of a fractal process, 
which is defined by the Fourier transformation of the conditional density T(r): 

n(fc) = 4tt / T(r)^^-r 2 dr . (105) 
J kr 

Because of a power-law form T(r) oc r~( 3 ~-°) for the conditional density of a stochastic fractal process, we also get 
for the generalized power spectrum a power law H(k) oc k~ D . Here D is the fractal dimension, which hence may be 
directly measured from the power spectrum. The generalized power spectrum is a suitable tool for describing highly 
irregular systems with fractal properties. 

3.7.3 Role of geometry of a sample 

In the case of a slice-like survey for scales larger than maximum sphere completely embedded in the sample, the 
geometry of considered structure become effectively be close to the case of an intersection of a 3-d fractal structure by 
a plane. Hence according to the theorem on intersections of fractals (sec. 2. 4. 5) it is expected that for a structure with 
the true fractal dimension D the measured fractal dimension of the sample will be Di nt = D — 1 . It means that there 
are three characteristic intervals of scales for behaviour of power spectra P(k) or TL(k): 1) for scales A < R%£ x we get 
P{k) oc k~ D with exponent equal the true fractal dimension; 2) for scales R^x < A < i? s we get P(k) oc fc _ ^ D_1 ' 
with exponent equal the fractal dimension of intersection; and 3) for scales A > R s we get P(k) — > due to a finite 
size of the sample. 

3.8 Multifractals and luminosity function 
3.8.1 Spectrum of fractal dimensions 

There are several different definitions for the fractal dimension: Hausdorf dimension Dh, box or capacity dimension 
Dt, correlation dimension D 2l mass-length dimension D m (for brief descriptions for astronomers see e.g. Gabrielli 
et al. 2004; Martinez & Saar 2002). For self-similar distributions all these dimensions are usually equal (and are 
< d, the dimension of the embedding space). In such simple cases we shall use the symbol D for this common fractal 
dimension. For multifractal analysis which operates with a whole spectrum of fractal dimensions all dimensions should 
be considered separately. 

We have already encountered the mass-length fractal dimension D m in the Fournier world model which was 
constructed so that M(r) oc r and hence D m = 1. It was the first mathematical model of a regular fractal structure 
applied to the whole universe. 

The box dimension or capacity of a set S C W is defined as 

D 6 = lim^M (106) 
e^o log(l/e) 



where N(e) is the minimum number of p-dimensional boxes of size e needed to cover completely the set S. 

In the case of spatial galaxy distribution p = 3 and one should extract from a sample the cube (3-d box) within 
which to perform the calculation of the box dimension £V The slope of the plot of logiV(e) versus log(l/e) gives an 
estimation of the box dimension. 

The correlation dimension D2 is defined via the correlation exponent 7 in the complete correlation function or 
conditional density T. Our discussion in the previous section was intended to distinguish between homogeneity and 
scale invariant properties and, for this purpose, it remains perfectly appropriate even if the galaxy distribution were 
multifractal. In this case the correlation functions we have considered would correspond to a single exponent of the 
multifractal spectrum, but the issue of homogeneity versus scale invariance (fractal or multifractal) is exactly the 
same. 

We emphasize also that a discussion in the literature around the multifractal properties of galaxy distribution had 
its origin in the difference between the values of the fractal dimension obtained by two methods: box-counting and 
correlation function (Jones et al. 1988; Martinez & Jones 1991; Martinez & Saar 2002). However, as we have seen 
above, the difference was caused by the use of the ^-function estimator which gave a distorted value for the correlation 
dimension Z?2- If one uses the appropriate T-function estimator then the unique fractal dimension D w 2 well fits 
observational data and hence eliminates the claimed need for that kind of multifractality. 

3.8.2 Luminosity function 

Above we considered fractal structures made of identical particles with unit mass. In connection with galaxies, Sylos 
Labini et al. (1998) and Gabrielli et al. (2004) use the term multifractality in the case where the so-called fractal 
support (positions of particles) has a unique dimension. In this case a stochastic process labels the fractal support 
particles by different values of some random quantity n (for example luminosity L or mass ra). 

Real galaxies have different masses and different luminosities which may be characterized by a luminosity function 
(f){L). This is usually represented by Schechter's law 

<j>(L) dL = (p* ■ (L/L*) a ■ exp(-L/L*) dL , (107) 

which gives the fraction of galaxies in the unit volume in the luminosity interval (L, L + dL). The parameters a and 
L* are to be determined from observations and (j>* is the normalizing coefficient from the condition <p(L) dL = 1, 
so that 

^r(a+^) (108) 

Here T(n, x) is the incomplete gamma function and (3 is a parameter which determines the luminosity truncation in 
the faint tail. 



3.8.3 Space-luminosity correlation 

Let us consider a realization of a stochastic process for which one may define the luminosity (or mass) density fi in 
the form 



N 

fi(x) = Mi S(x-Xj). (109) 
i=i 

where ^ is the luminosity (or mass) of i-th particle. 

Pietronero (1987), Coleman & Pietronero (1992) and Sylos-Labini & Pietronero (1996) regard luminosity or mass 
density functions as multifractal measures on the set of realizations. Multifractals are characterized by the spectrum 
of fractal dimensions D(L) which determines the dependence of the fractal dimension on luminosity or mass. This 
means that luminous and faint galaxies may have different spatial distributions. One prediction of the multifractal 
model is that the fractal dimension decreases for increasingly luminous objects. 

Let NL,s(x a , X, r) be the number of galaxies having the luminosity in the interval (X, X + AX) in a spherical 
shell S(r) = Airr 2 Ar with its centre at the point x a belonging to the structure. In order to generalize the concept of 
conditional density to the case of particles of different luminosity or mass, we define the conditional luminosity (mass) 
density as 

KAr) = / ^#^ ) = ("») 



4irr 2 Ar AX . 

-, N , pL + AL pr+Ar 

-ffEssasTEjL / «*>*'"■■ (1U) 

As shown by Sylos-Labini & Pietronero (1996)for a wide class of multifractal stochastic processes the conditional 
luminosity density has the general form 

KX,r) = r (X)-r L (r)= (112) 
L ^ Q .exp(X/X*(r)).^^. r D ^- 3 . (113) 



X*(r) / v ' v " 4?r 

This means that the Schechter law for the luminosity function is an observable consequence of multifractality and not 
just a convenient analytical form! 

Special features of the expression l|113|) are the dependence of the "knee" of the LF, X*(r), on the radius of the 
volume and also the dependence of the fractal dimension on luminosity D(L). These properties may be used for 
testing multifractality. However, this needs large samples of galaxies, because fractal analysis must be made for each 
luminosity interval separately. If the distribution is multifractal then the brightest luminosity X max in a sample is 
related to its spherical depth R s by the relation (Coleman & Pietronero 1992) 

X max cx flf , (114) 



where the exponent /3 depends on multifractal spectrum. 



4 The epoch of angular-position galaxy catalogues 

The history of main events illuminating the fractal debate in the 20th century is presented in the Tab0]"The debate 
on large scale fractality" . This subject has always been at the centre of cosmological thinking, even sharp debates, 
because of its direct link to cosmological principles. Of course, inhomogeneities in our neighbourhood are evident 
and the border to the assumed uniformity must lie somewhere at a larger distance. During the history the border to 
uniformity has gradually shifted outwards from the crystal sphere of fixed stars to Newton's evenly scattered stars, 
and then to Hubble's uniform galaxy distribution. At our times the observed superclusters of galaxies have put the 
border of uniformity to a scale of at least 100 Mpc. 

Together with the border to homogeneity, the centre of the universe has moved from the earth to the sun, and to 
the Milky Way. One of the most famous events in 20th-century astronomy was the Great Debate between Harlow 
Shaplcy and Hcbcr D. Curtis in 1920 about the scale and structure of the Universe (Smith 1983). In fact, this debate 
heralded the change of the Milky Way into an ordinary galaxy, whereby the centre of the universe finally disappeared 
into the realm of galaxies. 

Meanwhile, a new debate around the nature of galaxy clustering and on the border of the large scale inhomogeneity 
emerged. For almost the whole 20th century this struggle was going on, involving such figures as Charlier, Einstein, 
Selety, Hubble, Lundmark, de Vaucouleurs, Sandage, Peebles and others. As always in astronomy and in particular 
in this new extragalactic field, the scarcity of available observations at any historical moment leads to uncertain 
interpretations of the data. The modern phase of this fractal debate concerns observational tests with specially 
dedicated galaxy surveys to find the spatial scale where homogeneity becomes reality. 

New aspects and instruments for tackling the inhomogeneity problem entered the scene, when the concept of 
fractal was introduced by Benoit Mandelbrot (1975, 1977, 1982, 1988). Fractals, self-similar structures with long- 
range correlations, had been disclosed in physics and were then extended to astronomical scales, from the solar system 
and interstellar clouds to clusters of galaxies. Fractals on the largest scales became observable entities thanks to new 
astronomical techniques for measuring redshifts for thousands of galaxies, together with theoretical recipes for powerful 
data analysis. The primary open questions are: 1) Where is the border of transition from fractality to homogeneity?, 
and 2) what is the value of the fractal dimension of the galaxy distribution? 

4.1 The birth of the debate 

4.1.1 Einstein Selety correspondence 

In his first paper on cosmology Einstein (1917) gave arguments for homogeneity. He notes that any finite stellar 
system will evaporate into infinity due to internal gravitational interactions. He also emphasized that stellar velocities 
are small, which spoke against large potential differences and supported uniform large scale mass distribution. For 
Einstein a strong theoretical argument was that the Poisson equation modified by a cosmological constant term A 

Acj) ~ \(f) = 4nGp (115) 

has a solution for the density p = constant: 



(116) 



Table 2: The history of the debate on large scale fractality. 



Years 


Authors 


Subject 


1900- 


Fournier d'Albe 


regular hierarchical models 


1920s 


Charlier, Selety, Lundmark 


criteria for infinite world 


paradoxes 


Einstein, Selety 


Mach, stability, middle point 




Shapley, Zwicky, Abell 


strong galaxy clustering 


1930s - 


Carpenter, Kiang, Karachentsev 


superclusters up to 100 Mpc 


1970s 


Neyman, Scott 


2-level hierarchical model 




de Vaucouleurs 


cosmic density-radius law 


clusters 




p{r) oc r~ 7 with 7 = 1.7 


and 












uniformity 


Hubble 


uniformity from galaxy counts 






log N(m) — 0.6m+ constant 




Ambartsumian, Holmberg, Fesenko 


variable extinction in MW 




Sandage, Tammann, Hardy 


linear Hubble law at < 30 Mpc 




Webster, Longair 


isotropy of radio sources 




Wertz, Bonnor, Wesson, Alfven 


physics of hierarchy 


1970s - 


Haggerty, Severne, Prigogine 


N-body dynamics in hierarchies 


1980s 


Totsuji, Kihara, Peebles 


w oc 6>~ ' 8 => £ oc r" 1 - 8 




Mandelbrot 


fractals, multifractals 


£ 


Baryshev, Perdang 


first evidence for D = 2: 


and 

r 


Lerner, Schulman, Seiden 


M(r), z(r), stability, percolation 


function 


Davis, Peebles 


7 = 1.8, r = 5 Mpc, 2000 galaxies 




Einasto, Klypin, Kopylov, Bahcall 


r depends on depth and type 




Pietronero 


the method of V function 




Pietronero, Sanders, Coleman 


the first fractal analysis of CfA 




Pietronero, Ruffini, Calzctti et al. 


explanation of ro(R s ) 




Jones, Martinez, Saar, Einasto 


D 2 = 1.2, D = 2, multifractal? 




Sylos Labini, Montuori, Pietronero 


fl«2 from all 3-d catalogs 


1990s - 


Wu, Lahav, Rees 


fractality at r < 30 Mpc 


2000s 


Teerikorpi, Hanski, Theureau et al. 


TF < 200 Mpc D = 2.2, KLUN 




Paturel, Teerikorpi, Courtois 


LEDA counts < 15 m : 0.44m, D = 2.2 


fractals 


Baryshev, Bukhmastova 


2-point column density: D = 2.1 


in 


Zchavi (SDSS team) 


£(s) oc s" 12 , 29 300 gal. 


3-d maps 


Hawkins (2dFGRS team) 


£(s) oc s-°- 75 , 200 000 gal., D = 2.25 






distortion by peculiar velocities? 




Gott et al. (SDSS team) 


500 Mpc Sloan Great Wall 



so that a static, homogeneous matter distribution is possible, which also explains the small stellar velocities. The 
homogeneity scale was supposed to be about the mean distance between stars. 

In his paper Einstein then extended this result to general relativity, obtaining the spherical world model with 
a finite radius of curvature. At the end of his article Einstein's mentions that he docs not ponder if his model is 
compatible with available astronomical observations. Later this crucial hypothesis of homogeneity came to be called 
Einstein's cosmological principle of homogeneity. 

The Austrian scientist Franz Selety was aware of both Einstein's homogeneous and Charlier's hierarchical models. 
In an article in Annalen der Physik, Selety (1922) argued that it is possible to construct hierarchical worlds which 
fulfil simultaneously the following conditions: 

• infinite space 

• infinite total mass 

• mass filling space so that locally there is everywhere a finite density 

• zero average density of the mass in the whole world 

• non-existence of a unique middle point or middle region of the world 

In fact, Selety was one of the first to realize that the cosmological principle of "no centre" may also be valid for 
hierarchical systems. He expressed the essence of hierarchic models in that the universe appears for an observer in a 
"molecular-hierarchic" system everywhere basically similar. He also raised the question of Mach's principle in such 
universes and argued that it can be fulfilled. He pointed out that in such models a zero average density for the whole 
universe exists simultaneously with its infinite total mass. 

Einstein (1922) quickly replied to Selety. He expressed his opinion that Mach's principle is not fulfilled in a zero- 
density universe. Selety (1923) did not agree with Einstein and once more discussed the crucial points of his model. 
Summarizing the arguments which were raised by Einstein and Selety, we see as main objections to hierarchical models 
in the 1920's: 

• Mach 's principle is not valid for a hierarchic world model with zero global density. 

• Large potential differences in a strongly inhomogeneous universe lead to too high a velocity dispersion for stars, 
which is not observed. 

• A hierarchic stellar system will evaporate and stars will fill the voids, leading to a homogeneous distribution. 

• A hierarchic world contains a preferred middle point. 

4.1.2 A retrospective view on the posed questions. 

In retrospect, one may say that intriguingly, all these arguments arc still actual in modern cosmology. 

First, Mach's principle links the inertial mass of a body to the large scale mass distribution in the universe. Only 
relative to those distant masses one can define the acceleration of a test particle. In fact, the nature of inertial mass 
is still a challenge for modern theoretical physics, including general relativity, where the rest mass of a particle is 



regarded as a relativistic invariant, independent of the cosmologically distributed mass around the particle. Thus 
Mach's principle cannot any more be considered as a reason to reject the hierarchic models. 

Second, the small velocity dispersion of stars is due to their motion in our Galaxy and not related to the universe as 
a whole. In hierarchic models with D = 1 there is a constant velocity dispersion for each hierarchic level (i> 2 oc M/r oc 
constant). In fact, the problem of velocity dispersion has moved from stars to galaxies. In recent years, there has been 
special concern why the velocity dispersion around the local Hubble flow is so small inside the highly inhomogeneous 
galaxy distribution (Sandagc, Tammann & Hardy 1972; Chernin 2001; Baryshcv, Chernin & Teerikorpi 2001). One 
possibility is that the A-term, which Einstein introduced in his 1917 paper, or its modern generalization dark energy, 
is responsible for the smooth Hubble flow. 

Third, the question of stability of hierarchical (fractal) structures of gravitating particles is one of the open modern 
topics of gravithermodynamics. Perdang (1990) and de Vega, Sanches & Combes (1996) concluded that a statistical 
equilibrium may be possible for fractal structures with D w 2. 

Fourth, the question of the middle point is interestingly related to the cosmological principle. A stochastic fractal 
structure does not contain a privileged centre. Fractals preserve important properties of the old hierarchical systems 
and are more realistic models of the real galaxy distribution (Mandelbrot 1989; Pietronero, Montuori & Sylos Labini 
1997). 

The arguments of Einstein continue to inspire physical questions. Now they are not reasons for rejecting inhomo- 
geneous world models, but define important directions to study fractals. 

4.2 Early arguments for galaxy clustering. 

In the meanwhile, Charlier (1922) prepared a sky map for the distribution of 11475 nebulae, taken from Dreyer's New 
General Catalogue and his two Index catalogues. This led him to conclude that: A remarkable property of the image is 
that the nebulae seem to be piled up in clouds. In the debate on the nature of nebulae, Charlier considered: it appropriate 
to regard the spiral nebulae as foreign Galaxies similar to our own. Thus he related the observed clustering to the global 
matter distribution in the universe, i.e. he saw in it evidence for his hierarchical world models. 

Charlier's successor in the professor's chair at Lund University, Knut Lundmark, in his doctoral dissertation written 
in 1919, had expressed his conviction that: We can certainly expect to find a very complicated structure in the doubtlessly 
gigantic universal system, which is formed by the spiral nebulae. Later one of his main occupation was to build the Lund 
General Catalogue of galaxies in order to study the real distribution of galaxies. However, this work on thousands of 
galaxies described on separate cards was never completed. 

4.2.1 Observations disclose clusters of galaxies. 

After the discovery of the galaxy universe it soon became clear that in addition to field galaxies there are pairs, groups 
and clusters of galaxies. For example, in the 1930's clusters of galaxies were already routinely used for extending the 
redshift-distance law to larger distances by Humason (1931) and Hubble & Humason (1931). 
Several observations were used as evidence for the large scale clustering: 

• Shapley's metagalactic clouds. 

• De Vaucouleurs's Local Supercluster and density-radius relation. 



• Abell's rich clusters and their super clusters. 

• Shane-Wirtanen clouds of galaxies in the Lick counts. 

Shaplcy initiated wide photographic surveys of galaxies. The Shapley-Ames catalogue of 1249 bright galaxies from 
the year 1932 formed the basis for de Vaucouleurs's Reference Catalogue in the 1960's. Inspecting the distribution 
of galaxy clusters, Shapley came to the conclusion that there are "metagalactic clouds" (today's superclusters) , for 
example in the constellations of Coma, Centaurus and Hercules. The Centaurus cloud is nowadays called Shapley's 
supercluster. It is interesting to mention that Clyde Tombaugh, the discoverer of the planet Pluto, noted as a by- 
product of his extensive planet searches the Perseus-Pisces supercluster. He counted 1800 galaxies in this elongated 
cloud which is now a much studied agglomeration of clusters of galaxies at a distance of about 100 Mpc. 

Harlow Shapley's book The Inner Metagalaxy (Shapley 1957) is an interesting and illustrative outcome of the 
2-d epoch, where he summarizes the work on the clustering of galaxies performed at the Harvard observatory. In 
a section "Introduction on depth surveys" (p. 77) Shapley writes: The distribution of galaxies on the surface of the sky is 
easily examined on any uniform collection of long-exposure photographs. An effective study, however, of the distribution in the 
line of sight requires much greater labour. It is complicated by the difficulties of nebular photometry as well as by uncertainties 
introduced through the considerable dispersion in the intrinsic luminosities of galaxies. Systems side by side in space can differ 
by five or more magnitudes in apparent brightness, as for example the Andromeda nebula and its companions; and a pair with 
equal apparent brightness may differ in distance by a factor of ten. [. . .] In the study of the radial distribution of population 
it is necessary to use photometric methods for estimating distances, relative or absolute. A modern reader in the middle 
of large redshift catalogues is stricken by the fact that nowhere the redshift is mentioned as a possible indicator of 
distances. But at that time the local universe was not yet a subject of redshift surveys, the precious telescope time 
went to extending time demanding redshift measurements to deeper space. But after the first and second Reference 
Catalogues by de Vaucouleurs and collaborators were published in 1964 and 1976, with their compilations of hundreds 
and thousands of rcdshifts, many people started experimenting with the redshift in order to see the galaxy distribution 
in the radial direction. 

De Vaucouleurs (1953,1958) presented evidence, from the Shapley-Ames catalogue, for a local supercluster of 
galaxies centered at the Virgo cluster, and having an overall diameter of 30 Mpc. This system basically causes the 
well-known asymmetry in the number counts in the two hemispheres. It should be noted that there were also views 
regarding this enhancement of galaxy number density as a chance fluctuation (Bahcall & Joss 1976). 

A large increase in the number of known galaxy clusters came with George Abell's (1958) catalogue of 2712 rich 
clusters of galaxies. According to Abell' s selection criteria, a richness class is based on the counted number N of 
galaxies that are not more than 2 mag fainter than the third brightest galaxy. Abell introduced six richness classes 
(0,1,.. .,5), so that for class 30 < N < 49 and for class 5 N > 300. This catalogue covers the sky north of declination 
—27. The rich clusters of the southern sky were catalogued by Abell et al. (1989). Together these compilations contain 
4704 clusters. Such clusters are rare, but can be seen from large distances in space. 

Abell's collection was one outcome of the photographic survey of the entire northern sky, made by the large 
48 inch Schmidt telescope at Palomar Observatory - an incredibly important observational programme which gave 
astronomers huge amounts of data about stars and galaxies. The nine hundred 60 x 60 cm copies of the Palomar 
Sky Atlas photographs were a basic tool of observational astronomy at observatories all around the world for decades. 
Now the question was: Do Abell's rich clusters form superclusters? 



4.2.2 Santa Barbara 1961 conference. 

The international conference "On the stability of systems of galaxies" held in Santa Barbara, California in 1961 was 
devoted to the Ambartsumian hypothesis of instability of stellar and galaxy systems. Leading specialists on extra- 
galactic astronomy presented their works on galaxy clustering on different scales from binary galaxies to superclusters. 
We wish to mention that here de Vaucouleurs (1961) already presented his study of the density-radius relation for 
galaxy systems of different scales (see below). He emphasized that there was no indication of constant density level 
being reached in the observed range, so a constant mean density might be reached only for scales 10 times larger 
than the Local Supercluster (i.e. 300 h^Q Mpc). At the same conference Abell (1961) described his work on the 
rich clusters of galaxies and their superclustering. Typical second order clusters contain 10 clusters and have sizes 
of SO/ifg 1 Mpc. In their conference summary Neyman, Page, and Scott wrote: Both Abell and de Vaucouleurs feel that 
superclustering is established beyond doubt, and that the dimensions of second-order clusters (clusters of clusters of galaxies) are 
30 to 60 Mpc. This means that clusters cannot be treated as isolated systems embedded in an isotropic, homogeneous medium of 
field galaxies... 

Here the statisticians Jerzey Neyman and Elizabeth Scott recognize that their classical 2-level hierarchy model of 
the galaxy clustering (Neyman & Scott 1952) is not enough to explain the real galaxy distribution, and some radically 
new concept of clustering is needed for describing the observations. 

Further evidence for the large scale clustering among galaxies were obtained from the Lick Observatory galaxy 
survey by Shane & Wirtanen (1967). Its results were reviewed by Shane (1975) in an important paper, close to the 
end of the 2-dimensional period of galaxy catalogues. He summarized: Clustering seems to be a general, if not a universal, 
property among the galaxies. We find larger aggregations comprising numbers of clusters that extend over linear distances up to 
30 Mpc. There is suggestive evidence of still larger assemblages of galaxies on a scale of 100 Mpc or more. 

4.2.3 The cosmological de Vaucouleurs law. 

A most interesting discovery in the 1930's was that the non-uniformities of the galaxy distribution possess an intriguing 
regularity Edwin Francis Carpenter, an American astronomer studied clusters of galaxies and found that their galaxy 
number density depends on the cluster size so that the density is smaller in larger clusters. He calculated that the 
number of galaxies N in a cluster grows with the size r as 

N(r) cx r 1 ' 5 . (117) 

Carpenter (1938) regarded this relation as a cosmic restriction so that a cluster of a given extent may have no 
more than a limited number of members. This relation extends from pairs of galaxies to large systems of hundreds 
of members. This showed for him that small groups and large clusters do not essentially differ: the objects commonly 
recognized as physical clusterings are merely the extremes of a nonuniform though not random distribution which is limited by 
density . . . 

A next step was made by Kiang (1967) from the Dunsink Observatory in Ireland. After analyzing the distribution 
of Abell's clusters together with computer generated artificial distributions, he concluded that: the model of simple 
clustering by uniform clusters fails to represent the world of Abell's objects, in the same way as it has failed in the world of 
galaxies. Here he refers to the model introduced by Neyman & Scott (1952), where galaxies occur in clusters which 
are uniformly distributed throughout space. Kiang put forward a hypothesis that clustering of galaxies occurs on all 



scales - there are no clear-cut hierarchic levels. He wished to visualize the arrangement of galaxies so that the various 
clusters interpenetrate each other, which may guarantee that the average density does not depend on the volume. He 
came close to the modern view of fractal clustering, but did not make the crucial step to the fractal concept, and 
regarded the average cosmic density the same everywhere. Later Kiang & Saslaw (1969) found that the clustering 
extends at least to scales of 100 Mpc. 

Karachentsev (1968) added an important new aspect to Carpenter's result. He estimated average characteristics 
of 143 systems from binary galaxies to supcrclustcrs. He found evidence that both luminous and total (virial) mass 
densities are decreasing with increasing size of a system. This showed for the first time that the mass-radius behaviour 
of the hidden mass is also a power law, but with the exponent different than for the luminous matter. 

De Vaucouleurs made the decisive step in recognizing the cosmological significance of the clustering of galaxies. 
Based on his previous work and the studies of Carpenter, Kiang, and Karachentsev, he calculated from new data 
the density of matter inside galaxy clusters of different sizes. The results of his thinking he published in 1970 in the 
article "The Case for a Hierarchical Cosmology". Following his earlier paper (de Vaucouleurs 1961), he suggested the 
existence of a universal density-size law in the galaxy universe: 

p{r) = p (r )(r/r )-^ (118) 

where p(r) is the mass density in the sphere of radius r, p and r are the density and radius at the lower cutoff of 
the structure, and 7 w 1.7 is the power-law exponent derived on the basis of the available galaxy data. 

De Vaucouleurs (1970) summarized all at the time known properties of the galaxy clustering concluding that: In the 
1930s astronomers stated, and cosmologists believed, that, except perhaps for a few clusters, galaxies were randomly distributed 
throughout space; in the 1950s the same property was assigned to cluster centres; now the hope is that, if superclusters are here 
to stay (and apparently they are), at least they represent the last scale of clustering we need to worry about. . . 

He considered two extreme cases for the behaviour of the density-radius relation as a test for the nature of galaxy 
clustering. First, the density may decrease smoothly and monotonically, which would happen if there is no preferred 
sizes of galaxy clusters. In fact, this is what now is called stochastic fractal distribution. Second, the density curve 
may have a series of peaks at several preferred scales for which he suggested the values of 10 kpc (galaxies), 100 kpc 
(pairs and multiplcts), 1 Mpc (groups and clusters), 10 Mpc (superclusters) and 100 Mpc (third order clusters). In 
retrospect we may note that later the correlation analysis of observations pointed at the case without preferred sizes. 

4.3 Early arguments for the homogeneity of the galaxy distribution. 

It is intriguing that parallel to the piling data on galaxy clustering, there was a growing number of arguments favouring 
the uniformity of the galaxy distribution. In this view clusters were considered as exceptional objects in the sky. We 
summarize below the evidence for homogeneity, which were presented before the 1970s: 

• Hubble's galaxy counts. 

• Fluctuations in the number of galaxies due to variable Galactic dust extinction. 

• Sandage- Tammann- Hardy argument from the local linear Hubble law. 

• Isotropic distribution of distant objects 




Fi gure 7; Hubble's (1926) early result on the counts of bright galaxies. The single point at 16 mag was based on old data by Fath. The 
straight thick line has the slope of 0.6, which corresponds to the homogeneous cosmological model. The dashed line presents modern data 
from the LEDA database showing a slope of 0.44 in the interval 10 m -r 14 m , corresponding to the fractal dimension D = 2.2 (Teerikorpi 
2004). 

4.3.1 Hubble's counts of bright galaxies. 

Hubble (1926), from his bright galaxy counts, concluded that these correspond to the expectation from homogeneity 
up to m — 16. This was widely regarded as evidence for a homogeneous cosmological model. Thus he used the 
number-magnitude relation for galaxies as a cosmological test. 

The fundamental equation of stellar statistics, when applied to objects having a power law radial number density distribution 
n oc r~"' or a fractal number-radius relation N oc r D , implies that the number of objects N(m) having the magnitude less than 
m follows the relation 

N(m) = ^-m + const . (119) 
o 

with the fractal dimension D = 3 — 7. This result does not depend on the luminosity function of the objects. Hence for a 
homogeneous distribution (7 = 0; D = 3) one obtains the classical Seeliger law 

N{m) = 0.6m + const . (120) 



Hubble (1926) found that for photographic magnitudes in the interval 8™ 1 -j- 12 m plus at the point 16.7 the "0.6m- 
law" was valid, though he noticed a small systematic deflection which he ascribed to either observational errors or "a 
clustering of nebulae in the vicinity of the galactic system" (Fig. [7J. In the interval 10 m -j- 13 m this result was also 
confirmed by Shapley & Ames (1932) in their catalogue, for the sum of the counts in the northern and southern skies. 

For the bright galaxies (m < 16), the major shortcoming was the lack of data for the magnitude interval 12 m 16 m 
(see Fig. [7J. We will discuss modern data in this interval in sect. 5.4.3. 

Hubble was convinced that his data on bright galaxy counts already shows homogeneity. He considered that 
clusters of galaxies contain only a small fraction of all galaxies, and the true spatial distribution is quite homogeneous. 
Hubble's conclusion had a strong impact on the theoretical cosmology, as expressed by Einstein in 1933: Hubble's 
research has, furthermore, shown that these objects are distributed in space in a statistically uniform fashion, by which the 
schematic assumption of the theory of a uniform mean density receives experimental confirmation (cited by Peebles 1980). 



From this time on the picture of a uniform galaxy field with clusters as rare fluctuations became a paradigm of the 
homogeneous galaxy universe. 

Later Hubble's bright galaxy counts were used by Sandage, Tammann & Hardy (1972) to show the incompatibility 
with de Vaucouleurs's hierarchical model with 7 = 1.7. For this 7 Eg II 191 implies that N{m) = 0.26m + const, and 
this clearly differs from the value 0.6m from the counts by Hubble and Zwicky, available at that time. Hence they 
demonstrated that the hierarchical model with the fractal dimension D = 1.3 cannot explain the galaxy number 
counts. Here we should emphasize that modern counts of bright galaxies in the range 10™ -j- 14™ display a slope 0.44 
which corresponds to D = 2.2. This will be discussed in sect. 5.4.3. 

4.3.2 Hubble's deep galaxy counts. 

Hubble (1934, 1936) extended the number counts up to the magnitude 21 in his massive galaxy count program on 
1184 photographs on random positions in the sky. Each plate covered 0.25 deg 2 . Now he did not confirm the 0.6m 
prediction at faint magnitudes, but instead the counts followed a 0.5m law. He did not abandon the hypothesis of 
homogeneity, but tried to explain this deflection as due to a redshift effect on galaxy magnitudes. He made a kind of 
K-correction and concluded that the 0.6m law may be obtained only in a non-expanding universe. In fact, to the end 
of his life (1953), Hubble regarded that the cosmological redshift may be caused by some other effect than expansion. 

Sandage (1995) has analyzed Hubble's counting programme and found three kinds of systematic errors which 
influenced Hubble's calculations: 1) there were systematic errors in Hubble's photographic magnitudes, making them 
increasingly too bright for fainter objects; 2) there was an error in the K-term as applied by Hubble; 3) when calculating 
the prediction of spatial volumes in the Fricdmann model corresponding to different redshifts Hubble used an incorrect 
kind of distance. In the light of these problems, one may now see that Hubble's old result on faint galaxies cannot be 
used as evidence neither for homogeneity nor non-expansion. 

4.3.3 Variable dust extinction. 

One argument against the reality of superclusters referred to our "dusty window" to extragalactic space. Because the 
early evidence for superclustering of galaxies came from the sky distribution, a natural objection was that one should 
take into account the variable light extinction in different directions of the sky, due to the cloudy dust distribution in 
the Milky Way. 

Victor Ambartsumian developed the theory of a fluctuating Galactic extinction and applied it to the counts of 
extragalactic nebulae (Ambartsumian 1940; 1951). The paper by Neyman & Scott (1952) was devoted to their well- 
known two level clustering model. With references to Charlier (1922) and Ambartsumian (1951) they emphasize that 
there are in principle two approaches to galaxy clustering: 1) galaxies are really clustered in space; 2) the apparent 
clustering in the sky is caused by variable extinction due to interstellar dust clouds. Which factor dominates for 
clusters and superclusters? Later Holmberg (1975) and Fesenko (1975) further studied the role of inhomogeneous dust 
in the apparent galaxy clustering and concluded that the observed clustering is essentially modified by dust. 

In the 1950s Zwicky (1955, 1957) proposed that also intergalactic dust, concentrated in clusters of galaxies, could 
cause extinction of background clusters. He studied the distribution of 921 clusters and found a deficit of clusters in 
regions around the Virgo, Coma and Ursa Major clusters of galaxies. This he interpreted as obscuration by intergalactic 
dust. 



Also Karachentsev & Lipovetskii (1969) derived from counts of background clusters a positive mean value for the 
light extinction in clusters of galaxies. Their result was 0.2 magnitudes in the B band and they pointed out that there 
might have been selection effects influencing Zwicky's own counts: e.g., seen behind nearby clusters it is more difficult 
for more distant clusters to fulfil the identification criteria. We note that Mattila (1978) measured the diffuse light in 
the Coma cluster and concluded that a part of it could be scattered light from intergalactic dust inside the cluster. 
Later he measured together with Stickel et al. (1998) the far-infrared emission of the dust in Coma. The emission was 
detected and it was calculated that it could cause only about 0.2 magnitudes or probably less of extinction. 

Incidentally, Teerikorpi (2002) showed from the reddenings of quasars seen through galaxy halos (optical spectra of 
those quasars contain narrow absorption lines at much lower redshifts) that there is about 0.2 magnitudes of extinction 
per halo. Inside compact galaxy clusters and around the cores of rich clusters such as Coma the galaxy halos may 
almost overlap in projection and one might on this basis alone expect an extinction of the order of 0.1 mag. 

It is interesting to note that Zwicky & Rudnicki (1963) concluded, taking into account the effects of interstellar and 
intergalactic absorption . . . that the results obtained confirm the assumption that the distribution of clusters of galaxies is uniform 
within a space whose indicative radius is of the order of 10 9 pc. They inferred that clusters have different sizes up to a 
maximum 40 Mpc and there is no superclustering. 

However, the argument from extinction loosed its power after massive measurements of galaxy redshifts led to the 
discovery of the very lumpy 3-d galaxy distribution. We note that the study of the galactic extinction continues to be 
relevant for many extragalactic subjects, including the fluctuations in the cosmic background radiation. 

4.3.4 The classical "linearity" argument by Sandage, Tammann &: Hardy. 

De Vaucouleurs (1970) expected that mass density fluctuations cause deviations from the Hubble law and in particular 
in his hierarchical model the Hubble expansion rate should be reduced by gravitation inside the Local Supercluster. 
His PhD student Wcrtz (1971) developed a Newtonian expanding hierarchical model and calculated the predicted 
deflection from the linear Hubble law: the Hubble "constant" would increase with increasing distance within 20 Mpc 
from us. 

By this time Sandage had collected redshifts for 82 first-ranked cluster galaxies, which splendidly allowed one to 
test the prediction of the hierarchical model. This test was performed by Sandage, Tammann & Hardy (1972) who 
confronted the observations with the calculations by Wertz (1971) and Haggerty & Wertz (1972) for the predicted 
deflection from the linear redshift-distancc law. The result of the test was striking: the linear Hubble law with Hq = 
constant was established at all tested scales, while the hierarchical model predicted so strong a deflection that one 
should not see any cosmological expansion at distances closer than 20 Mpc. Later observations only strengthened the 
argument from linearity: a smooth linear Hubble flow starts in the vicinity of the Local Group, already at 1.5 Mpc 
(Sandage & Tammann 1975; Sandage 1986, 1987). 

Actually this test had a deeper meaning than being just a probe of the hierarchical distribution of luminous matter. 
It demonstrated a paradox: both empirical facts were true, i.e. the strongly inhomogeneous galaxy distribution in the 
local universe and the unperturbed linear Hubble law at the same scales. Sandage, Tammann & Hardy suggested two 
possible solutions for this paradox: 1) the mass density parameter could be very small, f^o <C 1, or 2) there could be an 
invisible uniform medium of high density. In both cases, the perturbations of the Hubble law would be tiny. In fact, 
from this post-classical cosmological test a whole new approach was developed, studying the properties of the cold 



very local Hubble flow (Chernin 2001; Baryshev, Chernin & Teerikorpi 2001; Maccio, Governato & Horcllou 2004). 

De Vaucouleurs (1972) did present evidence for a curved Hubble law in the local galaxy universe (the Hubble 
"constant" inferred from brightest group members strongly increased from small to large distances). Teerikorpi 
(1975a) studied the reality of such a behaviour using spiral galaxies with known van den Bergh's luminosity classes 
L c . These also showed a similar increase in H, but with a systematic shift depending on L c . This gave the crucial 
hint that the phenomenon was not real, but was caused by a selection bias influencing distance measurements. In 
Teerikorpi (1975b, 1982) the bias was satisfactorily modelled both for luminosity classes with Gaussian luminosity 
functions and for the brightest group member criterion. 

In Teerikorpi (1975a) the selection effect and its influence on the V/R ("Hubble constant") versus V (radial velocity) diagram 
of luminosity classified was explained in a simple manner, which serves as an example how a strong bias may be caused by 
observational limits. 

Assume that there is a limiting magnitude mi for the galaxy sample. Then in the M logR (R = distance in Mpc) 

diagram only galaxies below the line M — mi — 255 log R are available for mapping the kinematics of the local universe. 

Now consider a class of galaxies with a true average absolute magnitude Mo, then the selection begins to affect distance 
determinations at least at the distance Ro where Mo = m; — 25 — logro. It is easy to show that if R is the erroneous distance 
calculated from Mo at the real distance R — V/H, then the lower envelope of the points in the V/R vs. V diagram is defined 
by V/R — V/Ro (when 7? > Ro). Note that the slope depends on Ro, i.e. on the Mo of the luminosity class considered, when 
the limiting magnitude is constant. 

4.3.5 Isotropy in a homogeneous universe. 

Classically, from isotropy around one observer, together with the Copernican cosmological principle ("all points are 
alike"), one may infer the global homogeneity of the cosmological fluid (Walker 1944; for a simple geometric argument, 
see Weinberg 1977, p. 24). Before the 1970's there were three major observational evidences for an isotropic matter 
distribution around us. 

First, Hubble's deep counts of faint galaxies (up to m w 21) did not show large differences in different directions of 
the sky, after correction for Galactic extinction. It is true that Shapley (1934) noted a considerable difference in the 
numbers of bright galaxies (up to m pg = 13) for the northern and southern galactic hemispheres. He also found that 
at m « 18 the galaxy numbers might vary across the sky by a factor of 2 on angular scales of 30°. However, these 
deviations were viewed as local fluctuations. 

Second, the thousands of faint radio sources in the early catalogues (4C, Parkcs, Molongolo, and others) were found 
to be uniformly distributed in the sky within statistical uncertainty (e.g. Holden 1966). 

Third, already the first measurements by Conklin & Bracewell (1967) and Penzias, Schraml & Wilson (1969) of 
the cosmic background radiation found a remarkable isotropy at the level AT/T < 10~ 3 on arc-minute scales. 

These isotropies have generally been interpreted as strong evidence for homogeneity at scales larger than 1000 
Mpc. Is it possible to have local isotropy within inhomogeneous distribution of matter? The answer is "yes" for 
self-similar statistically isotropic fractal structures. This question is also connected with a more general formulation 
of the Cosmological Principle, which we will discuss in the last section. 



4.4 Results from angular galaxy catalogues 

4.4.1 Main catalogues of galaxies and clusters. 

In the 1950's and 1960's the first deep wide area catalogues of galaxies and galaxy clusters appeared, based on 
inspection of the monumental Palomar Sky Atlas photographic survey. This was the golden age of angular galaxy 
catalogues, from which the distribution of galaxies was investigated without distance information from the redshift. 
This important 2-d epoch in studies of galaxy clustering came to its end in the 1970's when the data analysis of all 
available galaxy catalogues was completed. 

The main catalogues were the Shapley-Ames (1932) Bright Galaxy Catalogue, Abell's (1958) catalogue of rich 
galaxy clusters, the Reference Catalogue by de Vaucouleurs & de Vaucouleurs (1964), the Catalogue of Galaxies and 
Clusters of Galaxies by Zwicky et al. (1961-68), and the Lick counts in cells by Shane & Wirtanen (1967). 

4.4.2 The angular correlation function analysis. 

The main mathematical instrument for analyzing the above mentioned catalogues was the angular correlation function 
(CF) techniques, which is described in detail by Peebles (1980). The angular CF w(9) is defined analogously to £(r) 
as the excess probability relative to the Poisson expectation to find an object in the solid angle 5£l at the angular 
distance 6 from randomly chosen objects in the sample: 



Here n is the mean surface density of objects in the sky. 

A practical estimate of w(9) for a sample of N galaxies with known positions on the sky can be obtained by 
counting objects in rings of radius 9 and width 59. The normalization condition is 



where Q s is the solid angle of the survey. An important practical point is that if the 2-D projection of an inhomogeneous 
spatial distribution is close to homogeneous, then this method cannot detect this true inhomogeneity but regards it 
as a structureless Poisson distribution. We will see that this just happened with spatial structures having fractal 
dimension D > 2. 

4.4.3 The relation between angular and spatial correlation functions. 

The primary aim of the analysis of angular catalogues was to estimate the spatial correlation function £(r) from the 
directly measured angular distribution of galaxies. The relation between w(6) and £(8) was first derived by Limber 
(1953), already in the context of galaxies (a useful review is by Fall 1979). For small angles 9, Limber's equation is: 



Here <fi(x) is the sample selection function, defined as the fraction of galaxies per unit volume of space observable at 
a distance x from Earth. 

Ea ll23l mav be inverted analytically in order to obtain £(r). However, a negative side is that this procedure requires 
differentiation of observed data which always contain some noise amplifying errors (Fall & Tremaine 1977). 



5P = nSSl[l + w(6)]. 



(121) 




(122) 




(123) 



An important particular case, which is useful in practice, is a power law solution of Ea ll23l For the spatial 
correlation function in the form: 



f(r) = Br^ = Br-^-^ (124) 
the corresponding angular correlation function is simply 

w{6) = AO 1 -' 1 =A6 D ~ 2 , (125) 

where the constant A comes from observations and the constant B depends on 7 and the selection function <j> (Fall 
1979). 

Caution: The power law solution ('ea ll25|) exists only if 

7 > 1, D < 2 (126) 

As we already discussed in sec. 2. 4. 5 according to theorem on fractals projection, this means that the method does 
not work for fractal structures with the fractal dimension D > 2. However, below we will see that the real galaxy 
distribution is characterized by the value of D which is just within this critical range! 

4.4.4 Hierarchical D = 1.2 models for angular galaxy catalogues. 

In their pioneering work Totsuji & Kihara (1969) derived from the angular data of the Lick galaxy counts a power law 
angular correlation function 

w(0) = A9-°- 8 , (127) 



from which they finally obtained 



Mr) = , ^r^Mpcy (12g) 



r 

In 1973 Peebles started an extensive programme analysing all angular galaxy and cluster catalogues, and asserted 
that all catalogues are characterized by almost the same power law with 7 = 1.77 and tq — 5/i^ 00 Mpc in scale intervals 
0.1 10 Mpc. A review of these results is given in Peebles (1980; 2001). 

In fact, this was a rediscovery of de Vaucouleurs's cosmological density-radius relation at small scales (0.1 -7- 10 
Mpc). Totsuji & Kihara (1969) and Peebles (1974a) had for the first time found that the galaxy correlation function 
is a continuous power law with no peaks, as would be expected from preferred scales (case 2 in de Vaucouleurs 1970). 
Thus it naturally reflects the self-similarity of fractal structures. The exponent 7 = 1.8 corresponds to the fractal 
dimension D = 1.2. 

The new emerging picture of continuous hierarchy inspired the construction of protofractal models for the spatial 
galaxy distribution. Peebles (1974b) tried to find a general luminosity function of galaxies consistent with the power 
law hierarchy. Soneira & Peebles (1977, 1978) constructed a static hierarchical model of the galaxy universe. They 
used regular hierarchical model of Fournier's type with 12 levels of galaxy pairs, so that — 2 and k r = 1.76 and 
according to ea llll the fractal dimension D = 1.23 or correlation exponent 7 = 1.77. The positions of galaxies are 
assigned according to the pair hierarchy, the luminosity function and apparent limiting magnitude were used to create 
a magnitude limited sample, and finally the galaxies were projected on the sky of an observer. Soneira & Peebles (1978) 
compared this binary hierarchy with Lick galaxy counts and concluded that even such a simple model reproduced 



well the angular correlation function of the observed galaxy distribution. Soneira & Peebles (1977) for Zwicky et 
al. (1961-68) Catalogue of Galaxies and Clusters of Galaxies, inferred that there is no evidence for a spatially uniform 
population of field galaxies. This confirmed the general tendency for galaxies to appear only in clusters. 

4.4.5 Tallinn 1977 conference 

A landmark event in the study of the galaxy distribution was the IAU Symposium N79 The large scale structure of 
the Universe held in Tallinn, Estonia, September 1977. Leading astronomers presented observational and theoretical 
works on galaxy clustering and it was the first wide discussion of all existing arguments for and against homogeneity 
at cxtragalactic scales. 

Peebles (1978) presented a review of the angular correlation function analysis applied for all available angular 
galaxy catalogues. He emphasized that "a small systematic error in the angular distribution can be translated into a very 
large error in the estimate of the spatial clustering" and that "redshift data will allow us to avoid this problem". Peebles 
argued that at distances < 10 Mpc the galaxy clustering is described by power law, but at larger scales the galaxy 
distribution became homogeneous. He referred the first time to the book by Mandelbrot (1977), where the idea of the 
fractal was extended to extragalactic scales. 

Joeveer & Einasto (1978) presented an analysis of existing redshift data and concluded that the galaxy Universe 
has a cell structure with the mean diameter of voids about 50 h^ Mpc. The presence of large structures and holes 
of various sizes was also demonstrated at the symposium by de Vaucouleurs, Tully, Fisher, Abcll, Tifft, Gregory, 
Kalinkov and others. 

In his concluding remarks Longair (1978) noted that "Everyone seemed to agree about the existence of superclusters 
... systems on scales ~ 30 - 100 Mpc." However, when he refers to the results of angular correlation function analysis 
by Peebles' group who derived a homogeneity scale at 10 Mpc, Longair asserted "I am still a firm believer in the basic 
correctness of the covariance analysis". Though the results presented by Kalinkov's team about the third-order clustering 
point to the tendency of continuous galaxy clustering. Longair emphasized that "One wonders whether their existence is 
consistent with the isotropy of the distribution of extragalactic radio sources and of the microwave background radiation". 

4.4.6 First evidence for the D = 2 distribution 

Though from the angular correlation function analysis the value of the fractal dimension D w 1.2 (7 = 1.8) was derived, 
there were other evidences pointing at the fractal dimension D w 2. Historically, it is interesting that Lundmark (1927) 
made in effect the first observational estimate of the fractal dimension on the basis of Charlier's model, noting that 
the second criterion is fulfilled, which corresponds to D = 2. 

Baryshev (1981) discussed some observational and theoretical arguments in favour of the fractal dimension D w 2. 
Using such observational cosmological data as: 1) the galaxy number counts N{m) in a wide magnitude interval from 
2 m up to 24 m ; 2) the virial mass density — radius relation p vir (R); 3) the peculiar velocity dispersion — radius relation 
cr v (R), he concluded that a hierarchical model of galaxy distribution with 7=1 (fractal dimension D = 2) is consistent 
with the observations. 

In that paper a new theoretical argument on a special property of fractals with the dimension D = 2 was also 
presented. It comes from Bondi's (1947) consideration of the global gravitational redshift part of the cosmological 



redshift. For a homogeneous matter distribution in the case of z « 1 the gravitational cosmological redshift is: 



5<f»(r) 1 GM(r) 1 / r 



were 5<fi(r) — <j>(r) — (f)(0) is the gravitational potential difference between the observer and the source, r# = c/H is 
the Hubble radius. 

It is very important that from the causality principle it follows that the source is in the centre of the matter ball 
with radius equal to the distance r between the source and observer. Note that Zeldovich & Novikov (1984, p. 97) and 
Peacock (1999, problem 3.4) put the observer to the centre of the ball and hence got the gravitational blueshift instead 
of Bondi's gravitational redshift. However such a choice of the reference frame violates the causality in the process 
considered. Indeed, the event of emission of a photon by the source (which marks the centre of the ball) must precede 
the event of detection of the photon by an observer. The latter event marks the spherical edge where all potential 
observers are situated. 

Generalizing ea ll29l to the case of a fractal distribution where M(r) oc r D one may derive the following relation 
for the gravitational part of the cosmological redshift within the fractal galaxy distribution: 

_ 4^?p rg fry' 1 

where pa, vq are the density and radius of the zero level of the fractal structure (galaxies in our case). 

For the fractal structure with D = 2 the cosmological gravitational redshift is a linear function of distance: 

Zgrav(r) = 27r ^2° r ° T = ~~~ T ( 131 ) 

where H g is the gravitational Hubble constant, which may be expressed as 

H g = 2np r - (132) 

For a structure with fractal dimension D — 2 the constant = po r o may be actually viewed as a new cosmological 
fundamental constant. If the value of the constant (3 — 1/2tt g/cm 2 ( e.g. po = 5.2 x 1CP 24 g/cm 3 and r Q — lOkpc), 
then H g = 68.6 (km/s)/Mpc. So the linear Hubble law within the fractal structure with D = 2 is possible, though it 
requires a very large amount of fractal-like distributed dark matter. 

4.4.7 Why did angular catalogues lose the D = 2 structure? 

An explanation of why it is difficult to study a fractal structure with dimension D > 2 from angular distribution 
of galaxies was given by Baryshev (1981). If we model a part of the fractal as a spherical cluster of particles inside 
the radius R, then one can derive the surface distribution F(a) of the particles, projected on the sky, using the Abel 
equation 

F(a)=2 f R p(r) ^L= (133) 

where a is the projected distance from the centre of the sphere. 

For a power-law representation of a spherically symmetric fractal structure p(r) = po(ro/r) 3 ~ D oc r~ 7 it is possible 
to obtain analytical solutions in closed form for fractal dimensions D 3, 2, and 1. For a homogeneous ball (D = 3, 
7 = 0) 

F(a) = 2 PQ r *Jl-^. (134) 



For a structure with D = 2 (7 = 1) we get 

F(a) = 2 Pa r a 

Finally, for D = 1 (7 = 2) the surface density will be 



1 -W<) +] "" 



(135) 



F(a) = 2poro — arccos — . (136) 

<7 li 

Hence for the values of a « R the surface density behaves approximately as 

F(a) ~ const , for D = 3 , (137) 

F(a) ~ In a, for D = 2 , (138) 

F(a) ~a-\ for D = 1 . (139) 

Both for D = 2 and D = 3 the surface density varies slightly, which means that the projected distribution appears on 
the sky with a homogeneous surface density. 

This derivation demonstrates that the angular correlation function analysis becomes inefficient for structures with 
the fractal dimension close to or larger than 2, because the information on the 3-d structures with D > 2 is lost as the 
projected on the sky 2-d distribution approaches homogeneity. As we already emphasized this result is a consequence 
of the general theorem on fractal projections (sect. 2.4.5) with its critical dimension D pr = 2. 



5 Debate on fractality: the epoch of spatial maps 

During the 1980s several galaxy redshift catalogues (see Table 1) became available for the 3-d analysis of spatial maps. 
This brought into light unexpected news about different behaviour of the ^-correlation function for different samples 
of galaxies and clusters of galaxies. It was then realized within some research teams that the fractal galaxy structure 
can naturally explain the peculiarities of the measured ^-function, and also that the appropriate mathematical tool 
for the analysis of fractals is the method of T-function. Since then the fractal approach has given rise to new fruitful 
directions for observational and theoretical studies. 



5.1 The fractal breakthrough in the 1980s. 

After the concept of the fractal entered the scene of the LSS of the Universe, the cosmological community as if divided 
into two parts. One continued to use the ^-function method, which led to the conclusion about the absence of a galaxy 
"fair sample". The second, smaller, group started to apply the fractal approach based on the T- function method, and 
they obtained concordant results from different galaxy samples. 

5.1.1 Davis & Peebles ^-function analysis of the CfA sample of galaxies 

The paper by Davis & Peebles (1983; hereafter DP83) was in many ways classical. It presented the first systematic 
analysis of the CfA data by the £ correlation function method. The CfA catalogue was the result of the first large 
redshift survey of the Harvard-Smithsonian Center for Astrophysics (CfA), which was complete to m B — 14.5 in the 
sky regions (S > , b > 40°) and (S > -2.5° , b < -30°)). It contained 2400 galaxies with redshifts. 



DP83 extracted a volume-limited subsample with Mb < —18.5 + 51ogft,i o which contained 1230 galaxies in the 
Northern zone and 273 galaxies in the Southern zone. 

In DP83 used the concepts of redshift-space and real-space correlation functions were utilized, developed by Peebles 
(1980). They also took into account the peculiar velocities, using the method, which we considered in sect. 3.2.5. 
They found the rms line of sight peculiar velocity distribution 

a v (r) = 340 ± 40( !_)0.i3±o.o 4 km/g (140) 

Mpch wo 

From the 3-d rms peculiar velocity the "cosmic energy equation" gives (Peebles 1980; Eq.74.9): 

< Wpcc > 1/2 ~ SSOSl^kms- 1 (141) 

which may be used to estimate the density parameter Oo- For CfA DP83 obtained Oo ~ 0.2 for the component of 
matter clustered with the galaxy distribution on scales r < l/ij~ Q Mpc. 

DP83 suggested the ^-function estimator which later came to be called standard. The calculation of data-random 
pairs was introduced for a correction of the edge effect and to reduce the shot noise on small scales. 

The main conclusion of DP83 was that the real-space two-point correlation function after the projection has a 
power-law form £ oc r~ 7 

&p(r)=(^^) 1 " M , (142) 

in the surprisingly wide interval of scales 

lO/i^kpc < r < 10 h^ Mpc . (143) 

Together with estimated errors the ^-function parameters were r = 5.4 ± 0.3 /i^qq Mpc and 7 = 1.74 ± 0.04. It is 
important that for scales r > 10 h^ 00 Mpc the ^-function drops, changes sign, and starts to oscillate near the zero-level. 
From what we considered in sect. 3. 2. 3., this is exactly what is expected for the ^-function estimator. 

After this pioneering work, the values of the unit scale r ~ 5hi 00 Mpc (defined as £(r ) = 1) and the correla- 
tion exponent 7 w 1.8 have been generally considered as standard cosmological numbers (see Peebles (2001) in the 
conference "Historical Developments of Modern Cosmology" ) . 

As we discussed in sec. 3 the ^-function method actually gives a distorted value for the intrinsic power-law exponent 
Itrue of a fractal structure, hence these results contain systematic errors. However at that time it seemed that the 
first 3-d map gave results which are consistent with analysis of angular catalogues and that the galaxy distribution 
becomes homogeneous on scales larger than 20 h^ Q Mpc . 

5.1.2 The puzzling behaviour of the ^-function 

When the first rcdshift surveys were studied by means of the ^-function method, a new unexpected problem appeared. 
The characteristic length tq was found to be dependent on certain parameters of the samples, such as the depth of a 
survey, the type and luminosity of galaxies and clusters, and the mean separation between the objects in the sample. 

Contrary to the expectation that "the spatial correlation function of galaxies is quite small for separations greater 
than about 2Qh^ Mpc", it was found by Bahcall & Soneira 1983 and Klypin & Kopylov 1983 that the characteristic 
length r (and hence the amplitude of the ^-function) become quite large for clusters of galaxies. 



Bahcall & Soneira (1983) calculated the redshift-space ^-function for a complete sample of N = 104 Abcll clusters 
with the distance class (< 4), and obtained the following estimates of its parameters: 

rg' w 25 /if oMpc , 7 c/ w 1.8 (144) 

Klypin & Kopylov (1983) studied another sample of Abell clusters with z < 0.08 and > 30°. Their catalogue 
contained N — 158 rich clusters of galaxies including redshifts measured with the 6- meter telescope of the Special 
Astrophysical Observatory of USSR Academy of Sciences. The redshift-space ^-function for the sample had the 
following parameters 

^25/HooMpc, ~ 1-6 (145) 

These results revealed a significant discripancy between the unit scales r for galaxies (5 Mpc) and for clusters (25 
Mpc). 

Moreover, when the ^-function was calculated for superclusters of galaxies (Bahcall & Burgett 1986; Lebedev & 
Lebedeva 1988), even larger scale were found: 

rf w 60 ftfooMpc , j cl ~ 1-8 (146) 

According to these data the correlation length r increases from 5 to 60 h^ Mpc when one considers increasingly 
massive objects in the universe. 

An important property of this new effect was also found by Einasto, Klypin & Saar (1986) who studied the 
behaviour of r within galaxy and cluster samples having increasing volumes. They used a cubic geometry for galaxy 
samples with the edge size I — R s (which is the depth of a survey) and found an approximately linear relation beween 
the unit scale and the depth of the sample: 

r oc R s . (147) 

In the framework of Gaussian density fluctuations on a homogeneous background this behaviour of ^-functions is an 
enigmatic fact, and to explain it several possibilities were discussed (Kaiser 1984; Bardin et al. 1986; Davis et al. 
1988; Bahcall 1988). 

The most popular one is Kaiser's idea of biased galaxy formation, based on a possible relation between the corre- 
lation functions for galaxy clusters and for the underlying mass density field. Here clusters are regarded as rare high 
density spots in the density field, so that one might expect 

^clusters = 6£dcnsity (r) , (148) 

where b is the bias factor which is about 10 if ^density = ^galaxies- This means that galaxies in clusters are formed from 
rare peaks above some global threshold in the primordial density field. However, the validity of this explanation in the 
case of Gaussian density fields was recently critisized by Gabriclli, Sylos Labini & Durrer (2000). They demonstrated 
that the increasing sparseness of peaks over the threshold in Gaussian random fields does not explain the observed 
increase of the amplitude of the correlation function £(r). 

Other alternative possibilities, like local inhomogencities, corrections for Galactic extinction, and luminosity seg- 
regation (Davis et al. 1988), look surprising and demand careful future studies with much larger galaxy samples. 



5.1.3 Pietronero's solution of the mystery of ro 

A radically new interpretation of the observed ^-function behaviour was found by Pietronero (1987) within the fractal 
approach to galaxy distribution. In this classical paper he introduced the T-function method for the 3-d galaxy map 
analysis and derived the relation between the T and £ functions which we considered in sec. 3. 4. 

For a spherical (or cubic) galaxy sample with the depth R s and a fixed luminosity of galaxies Pietronero (1987) 
obtained for the characteristic scale ro, defined as £( r o) = 1, the relation 

Hence within a fractal model one expects a linear dependence of r on R s . An important consequence of cci ll49l is that 
the increasing amplitude of ^-function and the corresponding increase of r for samples with larger R s is not due to a 
larger correlation length, but is simply an artificial effect caused by the definition of the reduced correlation function. 
In Fig. 3 of Pietronero (1987) it was clearly demonstrated that the cause of the increasing amplitude of ^-function is 
the increasing depth of a sample within the fractal structure. 

If the picture of the universal fractal galaxy distribution is true, then the method of T-function, appropriate for 
fractal structures, will reveal a power-law behaviour for future still deeper galaxy samples. The requirement of spherical 
geometry for a sample is the most important restriction for deep galaxy surveys. 

The first analysis of the CfA galaxy redshift catalogue by means of the T-function method was performed by 
Coleman, Pietronero & Sanders (1988). They found a power-law T-function with 7 = 1.5 ± 0.2, for VL samples with 
Ng a i — 226 and 442. The fractal behaviour was detected on the interval of scales from 1 /ij~ 00 Mpc up to 20/iJ" 00 Mpc 
without any characteristic scale, contrary to the homogeneity scale ro = 5 /ij~ Q Mpc derived by the ^-function method. 

5.1.4 Cellular fractal structure of the Universe 

Just a few months after Pietronero's (1987) paper, the fractal interpretation of the "correlation length" versus depth 
was applied to redshift data by Calzetti, Einasto, Giavalisco, Ruffini, & Saar (1987). They confirmed the linear relation 
between ro and R s , when the depth changed from 5h^ Mpc up to 50h^ Mpc. 

Then in a series of papers of a team led by Ruffini a cellular model of the Friedmann universe was developed, 
where within cells with sizes of about 100 Mpc the distribution of galaxies has the fractal dimension D sa 1.2, and on 
larger scales the universe becomes homogeneous (Ruffini, Song, & Taraglio 1988; Calzetti, Giavalisco & Ruffini 1988; 
Calzetti, Giavalisco & Ruffini 1989) . Their main conclusion was that de Vaucouleurs's density law may be reconciled 
with the homogeneous Friedmann model if there is a maximum scale of fractality. 

The Ruffini at al. model is based on the assumption that there are massive dark matter particles, called "inos", 
which obey Fermi statistics and are responsible for the initial density fluctuations of the cellular structure formation. 
The characteristic value of the "ino" rest-mass-energy is 0.4 10 eV, and the corresponding density parameter is 
Qinos — (0.4 -r- 1). An expected value of the size of the fractal cell is about 100 Mpc, which is determined by the Jeans 
length at the epoch when "inos" decoupled from matter. The angular scale of the corresponding CMBR fluctuations 
is about 1 degree. 

However, they used the lower value of the fractal dimension D m 1.2 derived from the ^-function analysis, which 
is actually a distorted value of the true fractal dimension D ss 2 obtained by means of the appropriate T-function 
method. It would be interesting to reconsider their model for D w 2. 



5.1.5 The multifractal confusion from the D 2 = 1.2 estimation using the ^-function method 

Almost at the same time another group of astronomers started to apply the fractal approach for a description of the 
galaxy distribution. However, now we may see that their work was affected by a distorted estimate of the correlation 
fractal dimension from the ^-function method. 

Jones, Martinez, Saar & Einasto (1988) and Martinez & Jones (1990) noted that when they applied the methods 
of box-counting and minimal spanning tree for determining the Hausdorf dimension Dh of the CfA redshift survey, 
they obtained the value Dh = 2.1 ± 0.1. However, from the ^-function for the same galaxy catalogue they concluded 
that the correlation dimension differs from this value, D 2 = 1.2, hence "the Universe is not a simple fractal. It is a 
more complex structure, a multifractal" (Martinez & Jones 1990). An analogous conclusion was made by Klypin et 
al. (1989) and Balian & Schaeffer (1989). 

Here we have an example of how inappropriate methods of data analysis may lead to erroneous theoretical con- 
clusions. Indeed, as we discussed in sec. 3. 4 the estimation of the fractal dimension from the £ correlation function as 
D 2 = 3 — 7£ at scales close to ro gives a distorted value for the true codimension 7 because there 75 (r = r ) ~ 27. 
Hence to calculate the true value of the correlation dimension from the slope of the ^-function near r = r one should 
take into account this distorsion, so D 2 = 3 — 7 = 3— (f^(r = r )/2) = 2.1 for the observed slope -f^(r = r ) = 1.8. 

Therefore one may conclude that actually D 2 « Dh, i.e. the correlation dimension is consistent with the Hausdorf 
dimension for the CfA catalogue and there is no need for multifractality based on a difference between these dimensions. 
Modern results of T-function analysis for CfA, 2dF,SDSS and other galaxy redshift surveys confirm the value for the 
correlation dimension D 2 ~ 2, and hence eliminate the above confusion on multifractality (see sec. 5. 3). 

5.1.6 Balatonfured 1987 conference and observational evidence for very large structures 

Parallel with developments of statistical methods for the analysis of fractal galaxy distribution in the 1980s, new obser- 
vational evidence appeared supporting the existence of galaxy structures with sizes much larger than the characteristic 
homogeneity scale r derived from the ^-function analysis. 

The results of a decade of intensive research after the Tallinn'77 conference were discussed at the 130th IAU 
Symposium "Large Scale Structures of the Universe", held in Balatonfured, Hungary 1987. It happened that this 
symposium became the last one on cosmology for Yakov Zeldovich (1914-1987) and for one of his talented pupils 
Victorij Shvartsman (1945-1987) who had just started to investigate the large scale structure of the Universe at the 
Special Astrophysical Observatory (SAO) of the USSR Academy of Sciences. 

At the Balatonfured'87 conference new observational data on the reality of galaxy structures with sizes of about 
100 Mpc were presented. Huchra, Geller, de Lapparent & Burg (1988) discussed an extension of the CfA redshift 
survey, which for several years was the main test bench for different statistical methods of the galaxy distribution 
analysis. They concluded that empty regions (voids) and filaments are common in the observed galaxy distribution, 
and the sizes of voids achieve 50 /ij~ o Mpc, which is much larger than the characteristic ^-function scale ro « 5 /ij" Q 
Mpc. 

Karachentsev & Kopylov (1988) reported the results of a spectral survey of 245 galaxies with raj < 17.5 in a 
narrow strip which passed through the Coma cluster. They confirmed a bubble- like type of structure within the Coma 
supercluster and estimated the parameters of the £ correlation function: r = 22 h^ Mpc and 7 = 1.5. The average 
size of 14 voids was estimated to be about 25 h^ Mpc. 



The most prominent structures were discovered by studies of rich galaxy clusters. Tully (1986, 1987) analyzed 
the distribution of 47 Abell clusters within a region up to z = 0.1c and found a flat structure having a size of about 
300 

^ioo Mpc. This is called the Pisces-Cetus Supercluster Complex. 

Similar results for an even deeper survey of galaxy clusters were obtained by Kopylov, Kuznetsov, Fetisova & 
Shvartsman (1988). They presented the first result of the program "The Northern Cone of Metagalaxy" which included 
measurements of redshifts up to z = 0.28 for 58 rich compact clusters of galaxies inside the cone with b 11 > 60°. From 
these data they made a preliminary statement on the existence of inhomogeneities in the distribution of galaxies on 
scales up to 500 /if 00 Mpc. 

In the Summary of the conference made by Peebles (1988) one may read: There is considerable evidence of 
structure on scales > 50 hT x Mpc, but I think it is fair to repeat the old questions: could this be an artifact of errors 
in the catalogues? Could the eye be picking patterns out of noise? If the answers were definitely "no" it would be very 
damaging for scale invariant cold dark matter. We all will be following the debate with great interest. 

It took one more decade of hard observational work to disclose the reality of such super-large structures in galaxy 
distribution, but the debate is still going on how damaging it actually is for the CDM models. 

5.2 Further steps in the debate 

After the new idea of fractality entered the studies of large-scale galaxy distribution there was a period of strong 
opposition from those who used conventional methods of analysis. Fortunately in science the collision of ideas is 
actually needed for a deeper understanding of the universe. This was also the case with the fractal galaxy distribution, 
which gave an alternative to homogeneous matter distribution in the Universe. 

5.2.1 Princeton "Dialogues'96" : Davis's evidence for homogeneity at scales larger 20 /i^Mpc 

In 1996 an international astronomical meeting under the intriguing title "Critical Dialogues in Cosmology" was held 
in Princeton. Remarkably the first subject which opened the conference was the dialogue between Marc Davis and 
Luciano Pietronero on the homogeneity of the galaxy distribution. 

Davis (1997) presented the position "that there is overwhelming evidence for large scale homogeneity on scales in 
excess of approximately 50 h^ Mpc, with a fractal distribution of matter on smaller scales" . He emphasized that the 
observed correlation function £(r) is well characterized by a power law, £(r) w (r/r )~ 7 , with r w bh^Mpc and 
7 = 1.8. Hence the fractal dimension at scales r < r is D = 1.2. 

Davis' arguments for homogeneity were: 

• Dl. Isotropy of the CMBR, X-ray and radio source counts. 

• D2. Observed counts of galaxies for magnitude range 14 < m < 18 have slope 0.6m. 

• D3. The observed angular correlation function w{9) is reliable for recovery of spatial correlation function £(r). 

• D4- Analysis of four VL samples from 1.2 Jy IRAS redshift survey give for £(r) standard values of 7 and ro 
when the volume limiting radius is increased from 60 h^ Mpc up to 120 /i^~ 00 Mpc. 

• D5. "The end of greatness" seen from the LCRS redshift survey. 

• D6. Ly-ct clouds detected in qso absorption spectra appear to be very nearly uniformly distributed in space. 



After discussing his arguments Davis concluded: "The measured two-point galaxy correlation function £(r) is a 
power law over three decades of scale and approximates fractal behavior from scales of 0.01 hi 00 Mpc < r < 10 /ij~ 00 Mpc, 
but on scales larger than w 20 hi 00 Mpc, the fractal structure terminates, the rms fluctuation amplitude falls below 
unity, and the Universe approaches homogeneity, as necessary to make sense of a FRW universe." 

5.2.2 Princeton "Dialogues'96" : Pietronero's arguments for fractality 

Pietronero presented the statistical method (T-function analysis) which is relevant for study the fractal structures, 
compared it with the ^-function method, and demonstrated the first results of application of both methods to the 
available redshift catalogues (Pietronero, Montuori, Sylos Labini 1997). 

Pietronero's main arguments for the fractality of the galaxy distribution were: 

• PI. The projection effect of a spatial fractal structure may lead to the observed isotropy on the sky for the angular 
distribution of astrophysical sources. 

• P2. A small number effect in the counts of bright galaxies may lead to observed 0.6m-law even for fractal 
distribution. 

• P3. The method of £ correlation function gives artificially distorted values both for the fractal dimension D and 
for the homogeneity scale rhom- 

• P4- The method of T-function (conditional density) is appropriate for estimation of the true value of the fractal 
dimension D and for detection of the crossover to homogeneity. 

• P5. The T-function analysis of available 3-d galaxy catalogues, CfA, PP, IRAS, LEDA, LCRS, ESP, gives the 
value D = 2.0 ± 0.2 for the fractal dimension at scales up to the radius R%%£ x of the largest sphere that can be 
contained in the sample. 

• P6. The homogeneity scale is not yet reached in existing galaxy catalogues and may be as large as 150 Mpc 
(LEDA result) and even as 1000 /j^o Mpc (number counts of ESP redshift galaxy survey). 

Both sides of the discussion agreed that the galaxy distribution is a fractal structure at least within scales 0.1 -j- 
10 /i^qq Mpc. But they disagreed about the values of the fractal dimension (Davis for D w 1.2 and Pietronero for D w 2) 
and about the value of the homogeneity scale (Davis for Rhom ~ 20 /iT/qq Mpc and Pietronero for Rhom > 150 /ij~ 00 Mpc). 
Pietronero also emphasized that possible existence of uniformly distributed dark matter may reconcile Friedmann 
homogeneous model with observed visible fractal structure. 

5.2.3 The problem of sky projection of fractals 

One of the strongest argument for the value of fractal dimension D — 1.2 and the homogeneity of galaxy distribution 
in space at scales larger than ro — 5 h^ 00 Mpc was the claim that the analysis of galaxy angular catalogues led to just 
such values for the parameters of the angular correlation function. 

Indeed, starting with Totsuji & Kihara (1969) in all angular galaxy gatalogues the analysts found the universal 
behaviour of the angular correlation function £, a ng(0) oc 0~ a with a w 0.8. For scales r < r^ this would correspond to 
the fractal dimension D in 3-d space D = 2 — a w 1.2. As this is less than 2, then according to the fractal projection 
theorem (sect. 2. 4. 5) one can estimate the true fractal dimension from the galaxy distribution projected on the sky. 



Unfortunately, this logic has a flaw. Indeed, if the real spatial galaxy distribution is a fractal structure with D > 2 
we would not be able to detect it from observations of the angular distribution of galaxies. This is because according 
to the fractal projection theorem the projected distribution imitates a homogeneous surface density. Consequently, 
3-d maps are necessarily required to discover fractal structures with D > 2. 

It is an intriguing fact that the fractal analysis of modern extensive 3-d galaxy maps has revealed D m 2.2, just 
putting the fractal dimension into that critical interval. 

A decade after the first warning on the possible "conspiracy" of structures with D > 2 by Baryshcv (1981), detailed 
studies started to appear on the complex problem of angular projections (Dogterom $ Pietronero 1991; Coleman & 
Pietronero 1992; Durer et al. 1997; Montuori & Sylos Labini 1997; Eckmann ct al. 2003). 

The method of angular T -function. The best demonstration of the consistency of the observed angular and spatial 
fractal structures with D ss 2 was given by Montuori & Sylos Labini (1997). They studied 3-d maps together with 
the corresponding angular distributions from several redshift catalogues: CfAl, SSRS1, Perseus-Pisces, APM bright 
galaxies, and Zwicky galaxies. For an undistorted estimation of the correlation exponent in angular data they used 
the conditional surface density r ang 

1 dN(9)_B ang D 

W*) - W) -^-9 (150) 

where S{6)d6 is the solid angle element (S(9) rts 2ir9 approximated for small angles S « 1), N(6) — B ang 9 D is the 
number of galaxies in the "polar cap" with the radius 9, D is the fractal dimension of the 3-d structure, which with 
the condition < D < 2 coincides with the fractal dimension D pr of the projected structure according to the fractal 
projection theorem (see sect. 2. 4. 5). a is the angular correlation exponent related to the fractal dimension D as 

a = 2-L> = 7-l (151) 

The last equality follows from the relation between the spatial and angular T-functions. The first one is the usual 3-d 
T-function T(r) oc r~ 7 with 7 = 3 — D, and the second one is the angular T-function r ang cx 9~ a with a = 2 — D. 

The result of this angular and spatial T-function analysis demonstrated very clearly that the fractal dimension of 
the observed 3-d structure is 

Z) = 1.9 ±0.1 (152) 

which was derived from independent analyses of the angular catalogues and 3-d maps. The angular correlation 
exponent was a = 0.1 ± 0.1 and the spatial correlation exponent was 7 = 1.1 ± 0.1 for all above mentioned catalogues. 

The most amazing fact is that the previously derived "universal" value of the angular correlation exponent a w = 0.8 
is an artificial effect caused by the method of the angular ^-correlation function w(9) (sect. 4.4). This reduced 
correlation function gives systematically distorted values of the true correlation exponent due to the normalization 
condition. 

We see here again that the whole story of how to find the true correlation exponent revolves around the difference 
between the power-law behaviour of the complete correlation function and the corresponding non-power-law reduced 
correlation function. This happens both for angular and spatial distributions. 



5.2.4 Modern research topics related to large scale fractality 

The debate on the nature of the large-scale structure of the visible matter in the Universe inspired many astronomers 
and physicists to study different aspects of the fractality galaxy distribution. The spectrum of subjects is very wide 



and shows that the principal questions, already raised by Einstein(1917, 1922) and Selety(1922, 1923) on the properties 
of hierarchical cosmological models, now are under careful investigation and actually they generate new branches of 
cosmological physics. 

In Table 3 we present a list of main reseach topics together with corresponding references which demonstrate the 
variety of cosmological aspects touched by the fractality of the large scale structure of the Universe. We shall discuss 
some results of the studies in the sec. 6 of the review. 

5.3 Recent results from the £ and T functions analyses 

As we discussed above the crucial parameters for correlation analysis of a galaxy distribution are: 1) the average sep- 
aration distance between nearest neighbour galaxies R sep , 2) the radius of the maximum sphere completely contained 
in a sample R^axi 3) the absolute magnitude interval AM, of selected galaxies, and 3) the number of galaxies in a 
volume limited sample N ga i. To extract reliable information on the correlation exponent and a homogeneity scale 
one shoud be aware of the restrictions of the method used for the estimation. We shall see that within the common 
interval of applicability both £ and T functions analyses give compatable results. 

5.3.1 The redshift space £- and T- functions 

£- and T -function analysis of the 2dF data. The final release of 2dF galaxy redshift survey (Colless et al. 2001; 2003) 
opens a new possibility for performing different kinds of statistical analysis of large samples of galaxies. 

The 2dF galaxy redshift survey contains about 220 000 galaxies in two (NGP and SGP) narrow slices of about 
90° x 15° (SGP) and 75° x 10° (NGP) complete up to bj = 19.5, with the effective reshift z s « 0.15, and the effective 
absolute magnitude M s — 5loghioo —20.0, corresponding to the luminosity L s s» 1AL* (Norberg et al. 2002). 

These data were analyzed using the ^-function method of the reduced correlation function by Hawkins et al. (2003). 
The redshift-space correlation function was approximated by two different power-law forms (their figs. 5, 6, 7), first, as 



at scales 3 < s < 20 Mpc. For larger scales 30 -j- 60 Mpc, £ z (s) becomes negative. 

Such a behaviour of the £ correlation function is consistent with the eg 1851 It means that 2dF galaxy distribution 
within the interval of scales 0.1 < s < 3 h^Q Mpc may be considered as a fractal structure with the fractal dimension 
D = 3 — 7 = 2.25. For scales s > 3 Mpc the ^-function continuously changing its slope and at scale tq « 5h^Q Mpc 
the exponent becomes j ro = 27 = 1.5 in perfect accordance with eqE3 Hawkins et al. (2003, fig. 7) found that 2dF 
sample has £ z (s) which is similar to Las Campanas and SDSS large slice-like surveys. 

The T-function method was applied to the 2dF VL samples by Vasiliev (2004) and Vasiliev et al.(2005), where 
they obtained that the conditional density for the 2dF data have power-law with the fractal dimension D = 2.2 ± 0.2 
at scales 0.5 < s < 40 Mpc. 

The value D — 2.2 ± 0.2 of the fractal dimension is consistent with results obtained by Sylos Labini, Montuori 
& Pietronero (1998) for all at the end of the 1990s avalable galaxy redshift catalogues: CfA, Perseus-Pisces, SSRS, 




(153) 



for the interval of scales 0.1 < s < 3 7i 100 Mpc, and second, as 




(154) 



Table 3: Main research topics related to large scale fractality. 



Subject 


References 


conditional density 
new methods of data analysis 
2-point conditional column density 


Pictroncro 1987; Coleman, Pietronero 1992; Gabrielli ct al. 2004 
Bharadwaj et al. 1999; Best 2000; Martinez, Saar 2002; 
Baryshev, Bukhmastova 2004; Vasiliev 2004 


fractal dimension 
of galaxy distribution 
scales of fractality 
mass-radius relation for all scales 


Coleman et al. 1988; Klypin et al. 1989; Lemson, Sanders 1991; 
Jones et al. 1988; Martinez, Jones 1990; Jones et al. 1992; 
Sylos Labini, Montuori, Pietronero 1998 
Sidharth 2000; Teerikorpi 2001; Rost 2004 


local radial galaxy distribution 
number counts, normalization 
luminosity function, multifractals 


Sandage 1995; Teerikorpi et al. 1998; Teerikorpi 2004 
Baryshev 1981; Joyce, Sylos Labini 2001; Courtois et al. 2004 
Sylos Labini, Pietronero 1996 


dependence of correlation function 
on depth, luminosity 
type of object 
peculiar velocities 


Einasto ct al. 1986; Calzetti et al. 1987; Davis et al. 1988 

Norberg et al. 2001, 2002 
Klypin, Kopylov 1983; Bahcall 1988; Bahcall et al. 2003 
Zehavi et al. 2002, 2004; Hawkins ct al. 2003 


local fractal dimension 
linearity and coldness of 
the local Hubble flow 


Tikhonov, Makarov, Kopylov 2001; Tikhonov, Makarov 2003 
Sandage et al.1972; Sandage 1986, 1987; Karachentscv et al.1996 
Ekholm et al. 2001; Karachentscv et al. 2003a,b; Whiting 2002 


protofractal (hierarchical) models 

relativistic fractal models 
cosmological tests, fractal universe 


Wertz 1971; Wesson 1975; Soncira, Peebles 1977, 1978 
Bonnor 1972; Ruffini et al. 1988; Ribeiro 1993; Gromov et al 2001 
Fang et al 1991; Baryshev et al. 1994; Joyce et al. 2000 


cosmological gravitational redshift 


de Sitter 1917; Bondi 1947; Baryshev 1981, 1994 


local tests of cosmological vacuum 
and dark energy within fractals 


Chernin 2001; Baryshev et al. 2001; Axenides et al. 2002 
Chernin et al. 2004; Maccio et al. 2004 


statistical mechanics of 
self-gravitating fractal gas, D = 2 


Perdang 1990; de Vega, Sanchez, Combes 1996, 1998 
Combes 1998; Huber & Pfcnniger 2001; 


N-body simulation, 
local structure 
initial conditions and discreteness 
stability, velocity, force 


Govcrnato et al. 1997; Moore et al. 2001 
Klypin ct al. 2003; Maccio et al. 2004 
Bottaccio et al. 2002; Bacrtschiger et al 2002; 2004 
Gabrielli, Sylos Labini, Joyce, Pietronero 2004 


origin and evolution of 
large scale fractals 


Haggerty 1971; Peebles 1974a; Alfvcn 1982; Lerner 1986 
Ostriker et al.1981; Pietronero et al.1986; Schulman ct al.1986 
Szalay, Schramm 1985; Maddox 1987; Luo, Schramm 1992 


cosmological principle 
and fractality, 
isotropy and homogeneity 


Mandelbrot 1975, 1977; Pietronero, Sylos Labini 1995; 
Rudnicki 1995; Wu, Lahav, Rees 1999; Baryshev,Teerikorpi 2002 
Mandelbrot 1989; Sylos Labini 1994 



IRAS, APM-Stromlo, LEDA, Las Campanas, ESP. They used T-function analysis so the probed scales were limited 
by the radius of maximum sphere r^jf completely embedded in the geometry of a catalogue, i.e. about 20 /ij~ 00 Mpc 
for the existed slice-like surveys and 100 /ij~ 00 Mpc for LEDA sample. 

SDSS: results from the £ and F* analyses. The Sloan galaxy redshift survey with its million galaxy redshifts and 
the wide sky coverage of about tt steradians is the ideal catalogue to settle definitely the on-going fractal debate. Up 
to now, however, only narrow (2.5° -j- 5°) slices - like in the 2dFGRS - have been completed and published as EDR, 
DR1, and DR2 catalogues (see the web site of SDSS). 

Zehavi et al.(2002) performed the ^-correlation function analysis of a sample of 29300 SDSS galaxies with radial 
velocities 5700 < cz < 39000km/s and absolute magnitude interval of —22 + 51og/iioo < M r < —19 + 51og/iioo- 
The redshift space £ correlation function of this SDSS sample definitely has a non-power law (see their fig. 5). The 
authors took the scale interval 2 < s < 8/i^~ 00 Mpc where they approximated the £- function by the power-law £ z (s) = 
(s/S.O/i^qqMpc) -12 . However it is clear from fig. 5 in Zehavi et al.(2002) that the £ z (s) has three characteristic intervals 
of scales: 1) for the interval 0.1 < s < 0.5ft.7/ oMpc the exponent 7 w 1.8, 2) for the interval 0.5 < s < 5/i^~ 00 Mpc 
the exponent 7 w 1, and 3) for the interval 5 < s < 30 /i^Mpc the exponent 7 w 1.8, As we discussed in sec. 3. 4 
such behaviour of the £ correlation function is just as expected for the fractal structure if one takes into account the 
characteristic scales R sep , ro, Rmax- 

The first T* analysis of the SDSS Luminous Red Galaxy sample was recently presented by Hogg et al. (2004) and 
also discussed by Joyce et al. (2005). The LRG sample so deep (average z about 0.3) that the radius of the maximum 
sphere reached the value R^x ~ 100 h± 00 Mpc. Hogg et al. (2004) for a sample with N ga i = 3658 found that the 
r*(r) has power-law corresponding to the fractal dimension D « 2 for the interval of scales 1-^25 h^ Mpc. For scales 
25 -j- 70 h± 00 Mpc there is a deflection from the power law and at scales 70 4- 100 h^ Mpc the T* achieves the constant 
value. This was interpreted as a detection of the homogeneity scale Rhom ~ 70 h^Q Mpc for the LRG SDSS galaxies. 

It should be noted that in a slice- like survey, such as the above considred sample of LRG, an artificial homogenization 
is possible starting from scales of about 0.25 R^ax' when the independent spheres in transversal direction cannot be 
completely embedded in the sample volume. Hence the above finding of a homogeneity should be in future reconsidered 
for larger spherical volumes. 

5.3.2 The problem of the peculiar velocity field 

Results from the correlation analysis of 2dF and SDSS data show that the galaxy distribution at least within interval 
of scales 0.1 -j- 20 h^ 00 Mpc is compatible with a fractal distribution having the fractal dimension D w 2.2. 

As we discussed in sec. 3 the restoration of the real-space ^-correlation function from the observed redshift-space 
£ z (s) involves a procedure of projection which is defined only for structures with fractal dimension D < 2 (i.e. 7 > 1). 
However the observed value of the fractal dimension e.g. for 2dF galaxies was D — 2.25 (7^ = 0.75), which violates 
the necessary condition for the restoration of the real-space correlation. Therefore the results of papers based on the 
procedure of projection should be reconsidered by using appropriate methods of restoration. 

Ignoring the problem of projection and using the standard procedure of the restoration of the real-space ^-function 
for 2dF data ( Hawkins et al. 2003) and for SDSS data (Zehavi et al. 2002) it was found that the deprojected spatial 
correlation ^-function has the canonical slope 7 w 1.7 and the unit scale ro ~ (5-j-6) h^[ 00 Mpc for the distance interval 
0.1 < r < 15/i^oq Mpc. The derived characteristic pairwise velocity dispersion a v = 500 -j- 600 km/s which is slowly 



decreasing with increasing distance in the interval 0.1 ~ 10/ij^gMpc. 

Though the effect of peculiar velocities leads to certain distorsion of the correlation function, which should be 
studied by appropriate methods, there is observational evidence that the value of the distorsion at scales 14-10 Mpc 
should be small because of the small observed velocity dispersion a v < 100 km/s in the local Hubble flow (Sandage 
1999; Ekholm et al. 2001; Karachentsev et al. 2003; Whiting 2003). This result is also in contradiction with the much 
higher value of the velocity dispersion derived from the ^-function analysis. 

Comparison of these results with N-body simulations in the ACDM Hubble Volume is restricted by a priori unknown 
biasing factor between the simulated cold dark matter density (and velocity) field and the real baryonic distribution 
of luminous galaxies. The arbitrary bias as a function of scale, b(r), is still not predicted by the theory of the large 
scale structure formation - so there is freedom to choose this function so that simulations can fit any observations. 

5.3.3 The problem of the dependence of ro on R s , L, d and galaxy type 

First of all as we already discussed in sec. 3 the unite scale r is not a proper characteristic of galaxy clustering but 
contains artificial distorsions due to individual properties of the sample geometry and the total number of galaxies. 
Hence a study of the relation between r and other physical parameters of a galaxy sample is not a correct approach 
to galaxy clustering. 

Secondly, the method of the projected ^-correlation function used by many groups (Norberg et al. 2001, 2002; 
Hawkins 2003; Madgwick et al. 2003; Zchavi et al.2004a;b) for derivation of the value of the real-space "correlation 
length" r utilized the procedure which eliminates the fractal structures with D > 2. 

Therefore the observed relations between r and i? s , L, d, and galaxy type (Norberg et al. 2001, 2002; Hawkins 
2003; Madgwick et al. 2003; Zehavi et al.2004a;b) contain a mixture of artificial and real effects which are difficult to 
separate. 

An appropriate method of study of the correlation properties of the fractal galaxy distribution is the dependence 
of the fractal dimension on luminosity or galaxy type, which is expected for multifractal structures (Pietronero 1987; 
Sylos Labini & Pietronero 1996; Gabrielli et al. 2004). The main problems for such studies are the small solid angle 
in the slice-like galaxy catalogues and the small number of galaxies in narrow absolute magnitude intervals. It is 
important to consider intervals of distances defined by the characteristic scales of samples (sec. 3. 4. 3). 

5.3.4 Power spectrum and intersection of fractals 

In sec. 3. 7 we considered the method of power spectrum analysis of galaxy distribution. As examples of its application 
to real data we consider CfA and SDSS galaxy samples. 

CfA redshift survey. Power spectrum (PS) of the CfA redshift survey was discussed by Park et al.(1994). They 
extracted four volume limited samples with about one thousand galaxies per sample and with the depths 60, 78, 101 
and 130 h^ Mpc. 

The PS is well described by two power laws: 1) at scales 5 4 30 Mpc the spectrum is P{k) cx k~ 2A , 2) at scales 
30 -T- 120 hiQ Mpc the spectrum is P(k) cx fc -1 - 1 . Because the radius of the maximum sphere is about 30 ft7/ Q Mpc, this 
means that the observed behaviour of the P(k) is consistent with the fractal structure with fractal dimension D = 2.1 
up to scales r — 120 h^ Q Mpc. At scales r > R^ x the survey effectively becomes 2-dimensional and in accordance 
with the theorem of intersection of fractal structure (sec. 2. 4. 5) the expected fractal dimension of the intersection should 



be Di n t = D — 1 = 1.1, which is just observed. 

SDSS galaxy redshift survey. The three dimensional power spectrum of several VL samples of the early data 
release SDSS galaxy redshift survey was analyzed by Tegmark et al.(2004). At their fig. 22 decorrelated real-space 
galaxy-galaxy power spectrum is presented, and again two power-law presentation of the PS is possible: 1) at scales 
10 -j- 60 /if 00 Mpc the spectrum is P(k) oc fc~ 2 , 2) at scales 60 -j- 200 hi 00 Mpc the spectrum is P(k) oc Again as in 
the CfA case, this is consistent with the fractal structure having the fractal dimension D ss 2 at all considered scales. 
The same behaviour of the power spectrum P(k) was found also for the 2dFGRS sample (Tegmark et al. 2003). 

5.4 Other results of the fractal approach 

5.4.1 The two-point conditional column density 

In section 3.6 we considered a new method of the fractal approach which is based on the calculation of the probability 
density to find a particle along the line of sight under the condition that the line is ended by two points of the structure. 
This method allows one to extend the fractal analysis up to the depth R s of a survey with the slice-like geometry. 

The analysis of LEDA and SDSS galaxies (Baryshev & Bukhmastova 2004) and 2dF VL samples (Vasiliev 2004; 
Vasiliev et al.2005) with the two-point conditional column density method, essentially extends the interval of probed 
scales from 20 ft-^g Mpc (reached by the T-function method) up to about 100 h^ Q Mpc. The result of the analysis is 
that the fractal dimension D 2.2 for the whole interval of probed scales. 

5.4.2 The local number counts of all-sky bright LEDA galaxies 

The local galaxy number counts as a function of apparent magnitude have two important applications in cosmology. 
Especially in the epoch of angular galaxy catalogues, the counts were regarded as a strong argument for a homogeneous 
galaxy distribution. Also, the local counts serve as a reference value for deeper galaxy counts. Here we give modern 
results on the the counts of bright galaxies on the basis of the all-sky LEDA data base. 

The Lyon extragalactic data base (LEDA) created by Georges Paturel in 1983 is a continuation of de Vaucouleurs's 
Reference Catalogue and its later editions. In fact, The Third Reference Catalogue was already based on the LEDA 
data. The LEDA extragalactic database offers currently a catalogue of homogeneous parameters of galaxies for the 
largest available whole sky sample. Among its over million galaxies there are about 50 000 galaxies with a measured 
B magnitude brighter than 16 mag. 

Made from the amalgamation of all available catalogues and continually being completed with the flow of new 
data, the completeness of the LEDA sample has been studied over the years inspecting the counts (Paturel et al. 1994; 
Paturel et al.1997; Gabrielli et al. 2004). Simultaneously with completeness, the counts give information about the 
slope of the bright end of the galaxy counts, hence on the spatial distribution law. 

Recently, Teerikorpi (2004), in connection with a study of the influence of the Eddington bias on galaxy counts, 
investigated the LEDA galaxy counts in the B magnitude range 10 -j- 16. All such galaxies from LEDA were taken 
having its total B magnitude and its a given, and having galactic latitudes | b |> 25 deg. The analysis of the counts 
indicated a slope of 0.44 in the B range 10 -j- 14, which corresponds to the fractal dimension D = 2.2 up to scales 
about lOO/i^Mpc. 

Rather similar results were obtained by Courtois et al. (2004a), who calculated the regression line up to B = 16 
for a somewhat different LEDA sample, and derived the slope of the counts of about 0.5. 



5.4.3 Radial number counts from the KLUN sample 

The usual studies of the large-scale structure utilize redshift as distance indicator. However, it is possible to use 
other distance measures, such as the Tully-Fisher and Faber-Jackson relations. Teerikorpi et al. (1998) derived 
the radial galaxy spatial distribution around our Galaxy, utilizing over 5000 Tully-Fisher distance moduli from the 
KLUN program. First results give evidence for a decrease in the average density consistent with the fractal dimension 
D = 2.2 ± 0.2 in the distance range 10 -j- 100 h^Mpc. This may be regarded as a new, independent argument in 
favour of the results obtained by conditional density methods, based on redshift distances, which also lead to the 
fractal dimension D m 2 at the same scale interval(Sylos-Labini et al. 1998). 

The method of distance moduli, differs from usual galaxy counts in the sense that the TF distance moduli probe 
with a better spatial resolution the distribution of galaxies. A magnitude measurement provides a very poor distance 
estimate, while a TF distance modulus has an error of about 0.5 mag. Furthermore, the method developed by 
Teerikorpi et al. (1998) takes into account the incompleteness of the sample, when it is applied. It will be interesting 
to apply the method to the larger KLUN+ sample with its 20000 galaxies in the future. 

Conclusions from modern local number counts. The above results imply several important conclusions: 

• there is no 0.6m-law in the bright galaxy counts, hence no homogeneity up to IOO/i^qq Mpc: 

• the observed QAAm-law of bright LED A galaxies and radial counts of KLUN galaxies are consistent with a fractal 
structure having D = 2.2; 

• it is necessary to revise the normalization of deeper galaxy counts 

First, there is no slope of 0.6 for galaxies in the B magnitude interval brighter than 14, i.e. there is no uniformity 
in the galaxy distribution up to 5000 km/s (100/i^qq Mpc). 

Second, the observed number counts N(m) = 0.44m + const in LEDA, the radial distribution N(r) cx r 2,2 for 
KLUN sample, and conditional density for main redshift catalogues T(r) cx r _1 are consistent with a fractal galaxy 
structure having D — 2.2 ± 0.2. 

Third, the absence of homogeneity on small scales influences the estimation of average number and luminosity 
densities (Joyce & Sylos Labini 2001). In case of homogeneity, these densities are constant, while for a fractal structure 
they depend on the radius of the volume in which they are calculated. This means that a usual normalization of deep 
galaxy counts based on the local homogeneity should be revised to take into account the local radial inhomogencity. 

6 Why fractality is important for cosmology 

The discovery of the fractal structure in 3-d redshift galaxy surveys opens new perspectives for understanding the 
origin and evolution of the large-scale structure of the Universe. One of the most important unsolved questions of 
theoretical cosmology is how to describe non-analytical fractal sources of the cosmological gravity field. Also much 
more studies should be done on the fractal velocity field and its evolution. This new situation in cosmology requires 
a careful reanalysis of the logic and structure of modern world models. 

6.1 Basic elements of cosmological models 

Modern cosmology is based on the following "corner-stones" : 



• cosmological principles, 

• fundamental physical theories, 

• cosmological observational data. 

In all these parts the fractal approach plays an essential role. Einstein's Cosmological Principle is extended up to 
Mandelbrot's CP. Fundamental physical theories include the new fractal mathematics. Cosmological observations 
uncover the scale-invariant properties of galaxy distribution. 

6.1.1 Three major empirical laws in cosmology 

The 20th century witnessed three major steps on the ladder of key discoveries, which unveiled three cosmological 
empirical laws: 

• the cosmological redshift- distance law 

• the thermal law of cosmic microwave background radiation 

• the fractal law of galaxy distribution 

Advances in astronomical instrumentation and spectroscopy were necessary for the discovery of the galaxy universe 
and then the Hubble law of redshifts in 1929. The development of radio astronomical devices led to the discovery 
of the thermal ocean of 3K photons in 1965. Finally, as we have discussed in this review, gathering of thousands of 
galaxy spectra with dedicated telescopes and applying appropriate methods of analysis, revealed the fractality of the 
large-scale galaxy distribution in the last two decades. 

6.1.2 The theoretical basis of modern cosmology 

All four fundamental physical interactions - the strong, the weak, the electromagnetic, and the gravitational - are 
used for modelling and understanding cosmological phenomena and for predicting observed astrophysical effects. 

But a special role in cosmology, among the fundamental interactions, is played by gravitation, the true astronomical 
force. It appears as the dominating force in the universe, starting from planetary and stellar scales up to cosmological 
distances. Its is natural that the theory of gravitation is at the heart of cosmology. It determines the large scale 
evolution of matter in the universe. 

Moreover, a cosmological model itself is the particular solution of the gravity field equations and this is why the 
gravity theory lies in the foundation of the whole building of modern cosmology (see discussion in Feynman et al. 
1995; Turner 2002a, b; Peebles 2002, 2003). This also explaines why the astrophysical tests of gravity theories should 
be considered as crucial cosmological observations (Baryshev 2003). 

6.1.3 The standard cosmological model 

At the start of the 21st century the standard cosmological model continues to be the Friedmann (1922, 1924) model of 
expanding space. It is based on Einstein's general relativity (GR) and the cosmological principle of homogeneity. The 
hot big bang picture also includes the process of growth of local inhomogcncities, caused by the gravity of non-baryonic 
cold dark matter (CDM). The global dynamics of the universe is at the present epoch determined by the antigravity of 



dark energy. The "ordinary" matter (stars, gas, dust), with which cosmology was concerned most of the last century, 
now is thought to make mere 0.5 percent of the mass of the universe and therefore cannot influence its expansion. 

Equations for gravity field. Einstein's equations of general relativity gives the relation between the geometry of 
space-time and the energy-momentum contents of matter (we use the notations from Landau & Lifshitz(1971) ): 

$ti-~gl$l=^Ti, (155) 

in which -ft^ is the Ricci tensor, gn~ is the metric tensor, and 

T k = diag(e, -p, -p, -p) . (156) 

is the total energy-momentum tensor (EMT) of the cosmological fluid in comoving coordinates. EMT includes two 
components: 1) ordinary matter with positive pressure, and 2) exotic substance, called dark energy or quintessence, 
with negative pressure. 

Einstein's Cosmological Principle. The cosmological principle of uniformity implies that the density and pressure 
of the cosmic fluid are functions of cosmic time only: 

g(f,t) = g(t), (157) 

p(f,t)=p(t). (158) 
The total energy density e = gc 2 and pressure p are given by sums of the above mentioned components: 

£ = Em + £de , P = Pm + Pde ■ (159) 

Here index "m" relates to different kinds of ordinary matter (dark and luminous) with positive pressure and with an 
equation of state 

p m = (3e m , 0</3<l, (160) 

e.g. e m — ecdm + £baryons + £rad + e v . Index "de" (dark energy/quintessence) means the more exotic substance with 
negative pressure and with an equation of state 

Pde=W£ d e, -1<W <0. (161) 

As a particular case (w = — 1) this includes Einstein's cosmological constant A, which may be interpreted as cosmo- 
logical vacuum. Quintessence as a new cosmic perfect fluid with the equation of state (|161f) was proposed by Caldwell, 
Dave & Steinhardt (1998) and Zlatev, Wang & Steinhardt (1998) as a solution of the cosmological coincidence problem 
(Peebles & Ratra 2003). Recently, values w < — 1 have also been considered in the equation of state for the so-called 
"phantom energy" (Caldwell, Kamionkowski & Weinberg 2003). 

Friedmann equations. The Friedmann model is an exact solution of the Einstein's equations feci ll55fl . Under the 
assumption of uniformity feas ll57H58|l the Robertson- Walker line element for homogeneous Riemannian spaces may 
be written in the form: 

ds 2 = c 2 dt 2 - S(t) 2 d X 2 - S(t) 2 I k ( X ) 2 duj 2 (162) 

where dui 2 = d9 2 + sin 2 6 'dip 2 , Xi^i^ are the "spherical" comoving space coordinates, t is the synchronous time 
coordinate, Ik(x) — sm (x)) X> smn (x) corresponding to curvature constant values k = +1, 0,-1 respectively and S(t) 
is the scale factor. 



Einstein's equations fea 11551) in the case of homogeneity are directly reduced to the Fricdmann's equation, which 
may be presented in the following form: 

d 2 S AttG ( 3p\ 
Te Bianchi identity implies the continuity equation 

£=- 3 ( e+ J)| ( 164 ) 

which must be added to ea 11681 Because the Lagrangian comoving coordinate x does not depend on time, one may 
rewrite ea ll63l using the definition of the proper metric distance r = S(t)x, as another form of the exact Friedmann 
equation: 

#r_ GM g (r) ^ 



dt 2 

where the gravitating mass M g (r) is given by 



M g = M m + M r + M v (166) 



and contributions from matter, radiation and vacuum are 



M m (r) = % (f? m + % )r i .107) 



M r (r) = —2Q r r 3 (168) 
Air 

M v {r) = -—2g v r z (169) 

Solving the Friedmann's equation fEa 1165(1 one finds the dependence on time for the metric distance r(t) or the scale 
factor S(t). A classification of two-fluid matter-dark energy Friedmann models is given by Gromov et al. (2004). 

Cosmological parameters. The FLRW model has two main parameters. The Hubble parameter H = S/S and the 
deceleration parameter q = —SS/S 2 which for the present time to & re H(to) = Hq and q(to) — go respectively. 

One also frequently uses the density parameter Q = g/ g cr where the critical density is 

to - gj- (170) 



Ea ll65l mav be written also in the form 

q = ln(i + 

gc 



^->f l • ^r) 11711 



where Q,p, g are the total quantities, and 

Q = 1 + fl k (172) 

with n k = kc 2 /S 2 H 2 . 

The old standard model has the following parameters 

n = n {rn)0 = i, n v = o, qo = o.5 (173) 

The new standard model which is currently accepted has 

fi fa 1, fi m Pa 0.3, Q v « 0.7, q ~ -0.6 . (174) 

This means that in the Friedmann model the expansion of the present universe is accelerated and that the dominant 
force in the universe is cosmological antigravity of the vacuum. 



6.1.4 Fractal sources for gravity field 

The results of proper statistical analysis of spatial galaxy distribution have revealed that the galaxy distribution is 
essantially inhomogeneous at least up to scales about 100 Mpc. It means that in modern cosmology one should consider 
more general models than homogeneous Friedmann ones. 

Mandelbrot's Cosmological Principle. The global mass-radius relation M(r) is the main characteristics of the 
fractal matter distribution which determines the fractal density field p(r) as the sources of gravitational field. 

In the next section we shall consider Mandelbrot's cosmological principle , which is a generalization of the Einstein's 
CP for the case of inhomogeneous cosmological models with isotropic fractal structures , and may be written as 



It is important that here the variable r is the radius of a ball around each point of a structure (see sec. 2. 3). 

Fractal cosmological models are based on solutions of gravity field equations with the sources described by the 
fractal density law. 

Lemaitre-Tolman-Bondi model in general relativity. Lcmaitre-Tolman-Bondi (LTB) models are exact solutions 
of Einstein's equations for 1) spherical symmetry, 2) pressureless matter (dust) and 3) motion with no particle layers 
intersecting. Originally studied by Lemaitre, Tolman and Bondi, these models are the simplest generalization of the 
Friedmann- Robertson- Walker (FRW) models with a non-zero density gradient (Bondi 1947). 

The LTB model has been used for understanding the kinematics and dynamics of galaxies around individual 
mass concentrations. For example, Teerikorpi et al. (1992), and Ekholm et al. (1999) could put in evidence the 
expected behaviour in the Virgo supcrcluster: 1) Hubble law at large distances, 2) retardation at smaller distances, 
3) zero-velocity surface, and 4) collapsing galaxies at still smaller distances. 

Bonner (1972) was the first to apply the LTB model to the hierarchical cosmology. He used de Vaucouleurs' density 
law p <~ r~ 7 with 7 = 1.7. Ribeiro (1992, 1993) has developed a numerical method for the solution of the problem. 
A new approach to the solution of LTB equations for a fractal galaxy distribution with large scale asymptotic FRW 
behaviour was presented by Gromov et al. (2001). 

However, the application of LTB models to a fractal distribution leads to a conceptual problem, because the original 
LTB formulation contained a central point of the universe, around which the density distribution is isotropic. In a 
fractal distribution there is no unique centre, but every object of the structure may be treated as a local centre which 
accommodates the LTB centre. Every structure point is surrounded by a spherically symmetric (in average) matter 
distribution. 

In this sense, the application of the LTB model to fractals means that there is an infinity of LTB exemplars with 
centres on every structure point. Their initial conditions are slightly different, because for any fixed scale the average 
density is approximately constant. For different scales the density is a power law. This excludes geocentrism and 
makes possible the use of LTB models as an exact general relativistic cosmological model where expansion of space 
becomes scale dependent. A similar problem of an apparent centre exists even for Friedmann models where the rate 
of space expansion within distance r from a fixed galaxy is determined by the total mass of the sphere around this 
galaxy. For LTB fractal models the space expansion at distance r from a fixed point of the fractal structure is also 
determined by the average mass of the sphere around this point. 



Q(r, t) = g(r, t) 



(175) 



p(r, t) = p(r, t) . 



(176) 
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Fi gure 8; The Hubble law within the fractal galaxy distribution. From Baryshev et al. (1998). 

6.2 The origin and evolution of large scale fractals: challenges for theoretical models 

The observed fractality of large-scale galaxy distribution opens new aspects in the process of the structure formation. 

6.2.1 Hubble law within fractal galaxy distribution 

Two fundamental empirical laws which have been established from extragalactic data are in apparent conflict with the 
SCM. First, there is the power law density-radius relation (de Vaucouleurs law) which corresponds to fractal struture 
with fractal dimension D ss 2 up to the scales of about 100 Mpc. Second, Cepheids, TF-distance indicator and Type 
la supernovae confirm the linearity of Hubble's rcdshift-distance law within the distance scales 1 -j- 100 h^ 00 Mpc, just 
where the fractality exists. 

As we noted in sec. 4. 3 Sandage, Tammann & Hardy (1972) first recognized the contradiction between the observed 
linearity of the Hubble law and a possible hierarchical galaxy distribution. They used the existence of the very local 
Hubble law as an argument against the de Vaucouleurs hierarchical model, which was based on calculations of the 
redshift-magnitude relation performed by Haggerty & Wertz (1972) for hierarchical cosmologies. 

A strong deflection from linearity expected within the fractal inhomogeneity cell was confirmed by Fang et al. 
(1991). The fractal aspect of this problem was discussed by Baryshev (1992, 1994) and Baryshev et al. (1998), who 
emphasized that the observed linear redshift-distance relation inside the power-law fractal density distribution creates 



the problem which they called the Hubblc-dcVaucoulcur's paradox, i.e. the observed coexistence of both laws on the 
same scales (see Fig. 8). There are several possible solutions of the problem. 

Asymptotically homogeneous LTB models. Using the LTB model Gromov et al. (2001) found the necessary 
conditions for the linear Hubble law existing within the fractal structure with fractal dimension D = 2. The larger 
the scale of homogeneity Rhom the smaller should be the density parameter f2 m , e.g. for Rhom = 100 Mpc the linear 
Hubble law exists at distances r > 1 Mpc if the density parameter f2 m < 0.01. 

Open Friedmann model. Joyce et al. (2000) suggested an open Friedmann model where the decreasing density 
of the fractal matter component at distance Rhom becomes less than the density of the cosmic background radiation, 
which plays the role of the leading homogeneous component of the FRW universe. Hence the density parameter at 
the present epoch is fio = Q. ra d and space expands within the fractal structure. They also considered qualitatively the 
modifications to the physics of anisotropy of the CBR, nucleosynthesis and structure formation. 

Influence of dark energy. Chernin et al. (2000), Chernin (2001), Baryshev et al. (2001) suggested that the cold 
local Hubble flow is a signature of the dominance of the cosmological vacuum or dark energy. Recent cosmological 
N-body simulations by Maccio et al.(2004) confirmed that inclusion of cosmic vacuum in the calculations leads to 
significantly lower velocity dispersion in Local Volume-like regions than what happens without vacuum. 

6.2.2 Problem of the origin of the fractal structure 

Since the 1970s Haggerty (1971), who worked in the University of Texas' Center for Statistical Physics, has developed 
a Newtonian model for the creation of hierarchical self-gravitating structures. This was based on Prigogine's and 
Severne's study of a non-Markofnan kinetic theory of binary gravitational interactions with irreversible growth of 
correlational energy. 

A scenario of explosive origing of large-scale structures was suggested by Ostriker & Cowie (1981) where large 
cavities and the superclusters of galaxies at scales up to 100 Mpc. This picture of galaxy formation is close to the 
way star formation is viewed. Within this approach Schulman & Seiden (1986) applied the percolation theory and 
numerical simulations for the derivation of the correlation exponent of the resulting structures. They got 7 ~ 1 which 
corresponds to the fractal dimension D s=a 2. To obtain the commonly accepted value 7 = 1.8 they tried to find a 
mechanism for changing the initial slope. 

Models for origins of the fractal large-scale structures were considered by Szalay & Schramm (1985), Pietronero & 
Kupcrs (1986), Maddox (1987), Luo & Schramm (1992). 

Statistical mehanics and thermodynamics of self-gravitating gas was considered by Perdang (1990), de Vega et 
al. (1996, 1998), Combes (1998). The intriguing result is that Newtonian self-gravitating N-body systems have a 
quasi-equilibrium fractal state with a dimension of ~ 2, though several problems are still open and need more studies. 

N-body simulations in order to study the gravitational growth and possible origin of fractal structures from initial 
small density fluctuations were made by Bottaccio et al. (2002), Baertschiger et al. (2002, 2004). They discovered 
that the discreteness has a strong influence on the results of standard cosmological simulations and noted that the 
Hubble time scale is too short for developing observed fractal structures within the whole range of scales. 

Large density fluctuations even at high redshifts lead to a conflict with small temperature fluctuation of the 
CMBR and to the problem of formation time for largest structures. A possible solution is that the standard scenario 
of gravitational growth of large-scale structure from small initial density seeds should be revised. Also the observed 



anisotropy of CMBR can be essentially distorted by the interving matter, so the fractal initial conditions may not 
relate to the observed AT/T (Schwarz et al. 2004). 

6.3 The Cosmological Principle 

There is much confusion in the literature on how Mandelbrot's Cosmological Principle of Fractality and its relation to 
the Einstein's Cosmological Principle of Homogeneity should be understood. Below we compare these principles and 
show that the principle of fractality is a natural generalization of the principle of homogeneity. 

6.3.1 Einstein's Cosmological Principle 

Einstein (1917) applied his general relativity to the cosmological problem, i.e. for constructing a model of the universe 
as a whole. Not yet knowing about galaxies, he imagined a world filled with stars and argued that the stars have a 
natural spatial distribution, which is uniform: matter concentrations around any preferred centre should with time 
evaporate and disperse uniformly all over the universe. Later, in correspondence with Selety, Einstein rejected the 
hierarchical distribution of stars (see sec.4.1). 

Einstein postulated a uniform matter distribution, and put relativistic gravity into cosmology. The result was 
a world with uniform geometry. The Copernican cosmological principle has many faces, as has been interestingly 
discussed by Rudnicki (1995) in his book. Its one formulation "all places in the Universe are alike" is naturally 
fulfilled in a homogeneous world. Besides the absence of a centre, another plus-side of uniformity was a simplification 
of Einstein's equations, which permitted him to derive the static spherical world model. Finally, Fricdmann (1922) 
liberated the universe from this stiff state, allowing the uniformly distributed matter and space to expand. 

The name Einstein's Cosmological Principle for the hypothesis of the homogeneity of large scale matter distribution 
in the universe was coined by Edward Milne who analyzed the foundations of cosmology in the 1930s. Milne formulated 
also an observer-oriented version of the Cosmological Principle as: the whole world-picture as seen by one observer 
(attached to a fundamental particle or galaxy) is similar to the world-picture seen by any other observer. 

In these early years of modern cosmology there was no direct observational evidence for the uniformity of the 
universe and it was theoretical reasoning which guided the cosmologist. 

Modern view on the Cosmological Principle. Every cosmological model has its beginnings in cosmological prin- 
ciples, special kind of hypotheses which are regarded as valid for the whole, perhaps infinite universe, even though 
observations can be made only from a finite part of it. Some principles touch questions of epistemology and are tacitly 
supposed to be true, for example that the knowledge about the whole universe is accessible for us, and that the laws 
of physics are the same everywhere and produce similar things. 

In the current cosmological thinking the contents of the Cosmological Principle are rather deep and complicated. 
We point out three distinct statements which are usually regarded as a starting point for constructing cosmological 
models. 

CP1 Physical laws the same in all space and time 
CP2 Fundamental physical constants are true constants 

CP3 Particular physical properties, including measuring standards, are the same in all space and time 



When cosmology advances, these statements may require adjustments or even radical changes in the light of new 
observations and theoretical ideas. For example, multi-dimensional theories predict variations of physical constants. 
However, all modern observational tests strongly restrict any such variations (Uzan 2002). The requirement of the 
same physical laws (CP1) also includes the possibility of discovering new laws of cosmological physics which appear 
only on very large spatial, temporal and mass scales. An important principle "More is Different" was introduced by 
Anderson (1972), who emphasized that each new level of complexity of material systems introduces new laws for their 
behaviour. 

In modern cosmology the Cosmological Principle is understood in a narrow sense as a hypothesis on the large-scale 
distribution of matter in the universe. For example, the statement that matter has a homogeneous and isotropic 
distribution is at the basis of the standard cosmological model. Nowadays observations show that Mandelbrot's 
cosmological principle of fractality has become an important alternative for the spatial distribution of luminous matter. 

6.3.2 Derivation of uniformity from local isotropy 

Walker (1944), a British mathematician who worked closely with Milne, proved that uniformity follows from his 
hypothesis of "Local Spherical Symmetry" which supposes that isotropy exists locally about each point of a Riemannian 
manifold. 

A simple reasoning leading to homogeneity when there is isotropy around each point, may be found in the First 
Three Minutes by Weinberg (1977) . He shows how one can go from any one point to another arbitrary place 
along circle arcs on which the density remains the same. Hence the density is the same on every point. However, 
strictly speaking this conclusion is based on a hidden mathematical assumption of regularity, i.e. the existence of a 
smooth density around each point, only then from the left-hand side does the conclusion follow: Local Isotropy + No 
Centre + Regularity => Uniformity. Here "local isotropy plus no centre" means that all points are equivalent and 
around each point the density law does not depend on the direction (though it might depend on the distance from 
this point). The "regular" matter distribution is described by continuous, smooth mathematical functions. 

In the chapter "Simplifying assumptions of cosmology" in his Introduction to Cosmology Narlikar (1993) explicitly 
introduces the assumption of the smooth fluid approximation, which essentially means going over from a discrete 
distribution of particles to a continuum density distribution. This is what we call "regularity" . It means that one may 
use the concept of the mass density at each point of space, like in a fluid. It is thus the union of local isotropy, no 
centre, and smoothness which gives homogeneity. In this case the Cosmological Principle reduces to the statement of 
the uniformity of the matter distribution. 

6.3.3 Mandelbrot's cosmological principle 

In his book Fractals: form, chance, and dimension, Mandelbrot (1977) foresaw that galaxies are fractal-like distributed 
and gave the first mathematical description of the fractal properties of such a distribution. He recalls how around 
1965, his ambition was to implement the law of decreasing density with a model where there is "no centre of the 
universe" or "the center is everywhere" . 

Mandelbrot views the fractal galaxy distribution as a major conceptual step in the description of the cosmological 
matter distribution. It is a kind of synthesis of hierarchical structures ("thesis") and homogeneity ("antithesis"), 
essentially based on randomness. Indeed, there is a fundamental difference between true random fractals and stiff 



hierarchical protofractals. Into protofractals the hierarchy is injected "ex-nihilo", by defining explicitly its levels. But 
fractals internally contain a scale invariance (self-similarity) and the impression of a hierarchy follows as an unavoidable 
consequence. A useful example is the Levy dust which is created by a random walk process in which the direction of 
each step is chosen isotropically and the length of a step follows a certain probability distribution. 

Fractality carries within itself also a trace of uniformity. Within a fixed radius, i.e. for a fixed scale, every observer 
counts the same number of elements, on average. But upon changing the radius, a "new uniformity" is found with 
a new mean number density. Furthermore, there is no centre for random fractals - this is another "relic" from 
homogeneity. 

Thus Mandelbrot made the first step for genuine fractals in cosmology, generalizing Einstein's cosmological princi- 
ple corresponding to D = 3 by fractality, which allows a non-uniform galaxy distribution with D < 3. His "Conditional 
Cosmographic Principle" states that all observers see similar cosmic landscapes around them, but only under the con- 
dition that they make observations from a structure element (galaxy) . "Conditional" in the spatial context emphasizes 
that each observer occupies a material element of the structure (c.f. "conditional density"). 

Mandelbrot's cosmological principle of fractality - that the observers attached to the material structure elements 
are equivalent - is close to what Milne presented. Thus the fractality of the universe perfectly satisfies Milne's 
Cosmological Principle. It also automatically makes what Karachentsev (1975) has called "the ecological correction 
to the Copernican principle" - the real observer can live only on or close to a material celestial body. This is usually 
called the Weak Anthropic Principle. In this sense the Copernican cosmological crinciple is contained by Mandelbrot's 
principle of fractality and hence the assertions about an "unprincipled" fractal universe (see e.g. Coles 1998; Wu et 
al. 1999) are not true. 

Does isotropy always imply uniformity? The proof of uniformity is based on the density being smooth around all 
points, which is valid for regular distributions, but not for fractals. It is smoothness which wipes out fractality and 
makes uniformity. Thus strictly speaking from local isotropy and the principle of no centre one can infer a fractality 
of the structure, of which homogeneity is only a special case with D = 3 (see Sylos Labini, 1994). 

Of course, there is never an exact local isotropy around every observer, not even in a uniform world, and still less 
inside a fractal distribution. Instead one may speak of a statistical isotropy, so that the sky observed from any galaxy 
"looks much the same" . In particular the counts of galaxies as a function of magnitude are similar around each galaxy, 
because radial distributions of galaxies are similar. It is natural to conjecture that for distributions made of discrete 
points, there is a generalization of the above chain of reasoning: Statistical Isotropy + No Centre => Fractality. 
Statistical isotropy and no centre, without the restrictive assumption of smooth mathematics, would thus lead to an 
isotropic fractal. 

As considered in sec. 2.4. 5, the most important reason for an apparent isotropic celestial distribution of galaxies is 
the projection of a fractal structure with fractal dimension D > 2 . According to the theorem on fractal projections 
the resulting distribution will have the fractal dimension D pr = 2, which means homogeneity on a 2-d plane or isotropy 
on the celestial sphere. 

Two other factors which lead to apparent isotropy on the sky are lacunarity and the luminosity function. It is now 
known that the patchiness on the sky depends not only on the fractal dimension, but also on the lacunarity, which 
is a measure of how frequent large voids arc. Numerical simulations have shown that fractals with a small lacunarity 
can have rather smooth projections on the sky. The second factor which smooths out the patchiness, is the large 



differences in the luminosities of celestial bodies. As a result two objects with equal apparent brightness actually may 
have widely different distances. The mixing of nearby and distant objects hides clusters and fills in holes. For example, 
this decreases the celestial anisotropy for very distant radio sources. 

6.3.4 Towards Einstein Mandelbrot concordance 

We do know of genuinely uniform components of the universe: the observed photon gas of the cosmic background 
radiation, the ocean of possible low-mass neutrinos, and maybe more importantly, the suggested physical vacuum 
or dark energy. As the average density of the fractal matter decreases with increasing scale, there will eventually 
be a scale beyond which the density of the uniform component is larger than the density of the fractal component. 
Hence one may regard, after all, the universe as homogeneous on such scales. However, this is not due to the galaxy 
distribution, but because of the uniformity of the relativistic matter component! 

As to the fractal galaxy distribution, there are two alternatives - a finite or infinite range of fractality. True, 
there are no scale limits to a pure mathematical fractal. The name 'fractal universe' is often linked with an infinite 
fractal. Such a universe would have zero average density. But real physical objects usually have lower and upper 
cutoffs between which the fractal properties are observed. Thus it is also expected that the fractal galaxy distribution 
appears only within a finite interval of scales. One possibility is that it becomes homogeneous on some maximum scale 
flhom- Thus one may, as Mandelbrot did, allow for the possibility that the matter distribution may become uniform 
on large scales, while being fractal on smaller scales. With any uniform matter component, such as the photon-gas or 
the cosmological vacuum, the universe becomes homogeneous on a sufficiently large scale. Such a universe would have 
a non-zero average density. Thus the intuitions of both Einstein and Mandelbrot appear to have grasped fundamental 
features of the universe. 

7 Concluding remarks 

The historical milestones of the path of the study of the large-scale structure of the universe are outlined in Table 
2. There one may see that during the period up to the 1970s many pioneering observational results and theoretical 
considerations appeared. The idea of a hierarchical matter distribution originated in a mathematical form thanks to 
Fournier, Charlier and Selety. It is intriguing to see that the problems discussed by Selety and Einstein in the 1920s, 
presently define whole directions in large-scale structure physics (see Table 3). 

Observers discovered on the sky a very lumpy distribution of galaxies and prepared angular catalogues of galaxies 
and clusters, from which the first superclusters were detected up to sizes of 100 Mpc. These gathered clouds over 
homogeneous cosmological models and there was a sharp debate about the significance of such inhomogeneities. At 
the time when there was no extensive distance information, it seemed that convincing arguments for a homogeneous 
galaxy distribution were found. First, the number counts of bright galaxies appeared to agree with the 0.6m law of 
homogeneity. Second, fluctuations in Milky Way's dust extinction was a factor producing apparent clustering. Third, 
the linearity of the Hubble law at small scales was seen as evidence for the galaxy homogeneity at similar small scales. 

In the meanwhile, a conceptual breakthrough was the suggestion by Mandelbrot that the galaxy clustering may 
be described as a stochastic fractal. After the nature of spiral nebulae was settled, the discussion around the fractal 
dimension and the maximum scale of galaxy clustering forms the core of the new Great Debate. 



During those early years the correlation function method, used in statistical physics for the description of homo- 
geneous systems with small density fluctuations, was extensively applied to the angular distribution of galaxies. In 
an important step, Totsuji & Kihara (1969) discovered the power law dependence of the angular correlation function, 
with the exponent 7 = 1.8 and a small homogeneity scale r$ = Mpc. For a comprehensive description of the 

statistical results from the angular galaxy catalogue epoch see the review by Fall (1979) and the text-book by Peebles 
(1980). 

Since then the correlation function has been applied for many angular and spatial galaxy catalogues, and in the 
words by Kerscher, Szapudi & Szalay (2000), " the two-point correlation function became one of the most popular sta- 
tistical tools in astronomy and cosmology" . These analyses created the feeling of confidence that the galaxy distribution 
is uniform on scales larger than 10/i^qq. 

However, when the study of the galaxy clustering entered the new epoch of 3-d maps, a series of remarkable 
discoveries followed, which were unexpected in view of the 5 Mpc "correlation length" . The first findings of real 
spatial large scale structures up to 50 Mpc were reported at the Tallinn 1977 conference. After that larger and larger 
structures have been discovered up to our days, proving that the earlier detections of superclusters were not rare 
chance events, but actually probed essential properties of the galaxy clustering. 

A decisive step was the introduction by Pietronero (1987) of the statistical method which is proper for the analysis 
of a fractal distribution of galaxies. This fractal-inspired method of conditional density has revealed the scale invariant 
galaxy clustering with the fractal dimension D m 2 up to scales where the method is applicable (Sylos Labini et al. 
1998). 

Also, it has become clear that the classical 2-point correlation function suffers from two major drawbacks. First, 
due to the projection effect the angular correlation function cannot detect spatial structures with the fractal dimension 
D > 2. Second, due to the normalization condition, both angular and spatial correlation functions yield systematically 
distorted values for the homogeneity scale (6>o, ro) and the correlation exponent (jang, 7), as compared with the true 
values of the complete correlation function. The importance of these effects is seen from the derivation of the correlation 
dimension D 2 = 1.2 from the ^-method, when the T-method gives D 2 ~ 2, for the main 3-d catalogues. The effects 
explain why the results from the £ and T analyses have led to different conclusions about the homogeneity scale and 
the correlation exponent. These conflicting results on the value of the fractal dimension were also a cause for the 
controversy around the multifractal analysis of the galaxy data. 

New extensive data, 200 000 galaxies for the 2dF survey (instead of 2000 for the CfA) , together with the T function 
analysis show that undistorted values of the correlation exponent is 7 w 0.8 which corresponds to the fractal dimension 
D = 2.2 ± 0.2. The scale of homogeneity i? max has a lower limit corresponding to the maximal sphere completely 
embedded into the volume of the galaxy survey, which is about 40/iJ~ 00 Mpc for 2dF and 70/i u1 1 q Mpc for the first release 
SDSS. In order to obtain an estimate of -Rhom on larger scales, one must enlarge the solid angle of the survey. Also, 
indirect evidence exist, such as the 450/i^qq Mpc Sloan Great Wall and 300/ij" O o Mpc inhomogencitics in the 2dF QSO 
distribution, that i?hom could reach such large values. An application of the new method of the two-point conditional 
column density to the 2dF and SDSS data confirmed that Rhom > 100/iJ~ 00 Mpc. 

In order to obtain a reliable mass-radius relation M(r) from observational data it is essential to use appropriate 
statistical methods of analysis (like conditional density T(r)) and control the distance scales R sep and R^ax which 
determine the region of applicability of the used methods. 



Three large numbers of galaxy redshift surveys needed for further progress in the large-scale sctructure analysis: 
large solid angle on the sky, large number of galaxies and large depth/separation ratio. 

A special study of dark matter distribution is requiered, which now is possible by using gravitational lensing 
technics ( Mellier 1999; Wittman ct al. 2000; Gray et al. 2002; Refregier 2003). The first results of the weak lensing 
analysis have shown that the dark and luminous matter are spatially similarly distributed. 

We summarize the following main challenges that the modern galaxy clustering research faces: 

• due to projection on the sky the angular distribution of galaxies loses the information on structures with fractal 
dimension D > 2; 

• due to the normalization condition the £ function estimators yield incorrect values of the homogeneity scale i?hom 
and the fractal dimension D; 

• the crucial role in statistical analysis of a galaxy sample plays the three quantities: the average distance between 
galaxies R sep , the radius °f completely embedded spheres in the sample geometry, and the total number of 
galaxies N ga i in VL subsamples; 

• essential fluctuations in the derived value of D are caused by the "cosmic variance" (only one realization of the 
stochastic process) and the "sampling variance" (finite available volume); 

• the method of projected redshift- space correlation function cannot be applied to fractal structures with D > 2, 
hence new methods for extracting real-space correlation properties of the galaxy distribution are required; 

• the pairwise velocity dispersion 500 — 600 km/s on scales 1 — 10 Mpc, obtained by the ^-function method, is 
surprisingly high in comparison with the velocity dispersion a v < 100 km/s in the Local Volume with r < 10 
Mpc; 

• the differences in clustering of galaxies with different types and luminosities may require multifractal analysis; 

• understanding the origin and evolution of the observed fractal distribution of galaxies may require a revision of 
the standard picture of gravitational growth of structure, by considering the consequences of primordial fractal 
density fluctuations. 

We may conclude that the fractality of galaxy clustering has become a fundamental empirical phenomenon of 
observational cosmology, which should be explained by theoretical models of the Universe. 
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